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ABSTRACT 


V' 

Multivariable  design  techniques  developed  by  Professor 
Brian  Porter  of  the  University  of  Salford,  England  are  applied 
to  the  X-29A  forward  swept-wing  demonstrator.  The  aircraft 
model  is  in  a  linearized  state-space  form  at  a  number  of 
different  flight  conditions.  The  computer  aided  design  pro¬ 
gram  MULTI  is  used  to  develop  the  control  laws. 

Separate  longitudinal  and  lateral  controllers  are  devel¬ 
oped  for  each  of  eight  maneuvers  of  five  different  flight 
conditions*  0.4  Mach,  Ok  ft;  0.9  Mach,  50k  ft;  1.2  Mach,  15k 
ft;  0.7  Mach,  15k  ft;  and  0.9  Mach,  Ok  ft.  A  universal  later¬ 
al  controller  that  is  effective  for  all  maneuvers  at  all 
flight  conditions  is  developed.  Finally,  a  universal  lateral 
controller  requiring  only  three  distinct  gains  is  developed. 
This  minimum  gain  universal  lateral  controller  produces  mini¬ 
mal  degradation  of  response,  while  permitting  the  elimination 
of  gain  scheduling  usually  required  of  digital  flight  control¬ 
lers.  In  addition,  maximum  possible  maneuvers  are  performed 
without  exceeding  control  input  limitations. 

Suggestions  for  improving  MULTI  to  allow  artificial 
intelligence  algorithms  are  discussed.  These  would  enable  the 
designer  to  determine  the  number  of  point  designs  required  to 
satisfactorily  cover  the  flight  envelope,  and  at  the  same  time 
provide  the  designer  with  the  capability  of  designing  for  many 

points  with  minimal  effort. 

A. 


CHAPTER  I 


INTRODUCTION 


Background 

Future  aircraft  will  probably  use  digital  computers  for  a 
number  of  functions.  Integrated  fire/flight  control,  naviga¬ 
tion,  threat  assessment/management,  and  weapons  delivery  are 
areas  which  show  great  promise.  Digital  computers  offer  re¬ 
duced  weight  and  cost  with  improved  performance.  Future  air¬ 
craft  will  use  more  control  inputs  in  order  to  achieve  greater 
maneuverability.  An  increase  in  the  number  of  control  sur¬ 
faces  means  an  increase  in  pilot  workload  that  can  be  eased 
with  digital  flight  control. 

Future  fighters  will  achieve  this  greater  maneuverability 
only  with  a  corresponding  decline  in  static  stability.  This 
means  that,  without  proper  control,  future  fighters  will  be 
unstable.  Since  no  pilot  can  continuously  provide 
compensating  inputs  to  an  inherently  unstable  aircraft,  digi¬ 
tal  computers  will  have  to  control  the  aircraft  for  him.  The 
ability  of  computers  to  solve  complicated  algorithms  in  a 
short  time  also  offers  the  opportunity  for  new  maneuvers, 
called  con tro 1-conf igured  vehicle  (CCV)  maneuvers,  that  are 
now  possible  with  multiple  control  surfaces. 

Since  computers,  or  controllers,  will  stabilize  and  con¬ 
trol  the  aircraft,  the  question  arises  as  to  the  best  method 
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for  designing  these  complex  flight-control  algorithms.  Clas¬ 
sical  control  theory  designs,  which  examine  s ing le - inpu t/ - 
single-output  systems  (SISO),  are  ill-suited  for  multiple- 
input/multiple-output  systems.  As  more  and  more  inputs  are 
added  more  trial  and  error  is  required  in  the  design.  Linear 
quadratic  gaussian  (LQG)  techniques  can  be  applied  to  MIMO 
systems,  but  the  states  must  be  known  or  estimated.  Unfortu¬ 
nately,  no  known  algorithms  exist  for  choosing  the  proper 
weighting  elements,  and  much  trial-and-error  and  insight  ob¬ 
tained  from  experience  must  be  employed. 

Professor  Brian  Porter  of  the  University  of  Salford,  U.  K. 
has  developed  a  design  method  that  can  synthesize  digital 
flight  controllers.  Porter's  method  decouples  the  outputs 
while  minimizing  the  trial-and-error  or  tuning  required.  This 
thesis  is  a  third  generation  attempt  to  apply  Porter's  method 
to  aircraft.  Many  improvements  and  refinements  are  added  to 
the  computer-aided-design  (CAD)  program  MULTI.  This  thesis 
suggests  a  number  of  future  possible  improvements  to  MULTI 
that  would  eliminate  most  of  the  trial-and-error  presently  re¬ 
quired.  This  would  change  the  design  engineer's  role  to  that 
of  selecting  the  best  computer-generated  designs,  rather  than 
trying  to  generate  them  himself.  These  algorithms  offer  the 
ability  for  a  single  design  engineer  to  synthesize  a  longitu¬ 
dinal/lateral  controller  with  sensors,  actuators,  and  lags  for 
the  entire  flight  envelope  within  a  short  time  period. 

Grumman  Aerospace  Corporation  is  developing  a  demonstrator 
aircraft  under  the  auspices  of  the  Defense  Advanced  Research 
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Projects  Agency  (DARPA) .  The  demonstrator,  the  X-29A,  incor¬ 
porates  a  number  of  advanced  technology  features.  These  fea¬ 
tures  include  composite  materials,  variable  camber  wings, 
forward  swept  wings,  fly-by-wire,  and  close-coupled  canards. 
Problem 

This  thesis  first  develops  a  controller  using  Porter's 
method  for  each  of  six  longitudinal  maneuvers  at  five  flight 
conditions.  Next,  a  controller  for  each  of  two  lateral  maneu¬ 
vers  at  five  flight  conditions  is  developed.  A  robust  con¬ 
troller,  capable  of  producing  good  responses  for  both  lateral 
maneuvers  at  all  flight  conditions,  is  developed.  Finally,  a 
robust  lateral  controller  requiring  only  three  distinct  gains 
is  developed.  The  outputs,  control  surface  deflections,  con¬ 
trol  surface  rates,  and  states  of  interest  are  examined  and 
cosipared  for  different  flight  conditions  and  controllers.  The 
flight  conditions  chosen  are  0.4  Mach,  Ok  feet;  0.9  Mach,  50k 
feet;  1.2  Mach,  13k  feet;  0.7  Mach,  15k  feet;  and  0.9  Mach,  0k 
feet.  The  longitudinal  maneuvers  performed  are  direct  climb, 
two  direct  lifts,  pitch-pointing,  and  two  vertical  transla¬ 
tions.  The  lateral  maneuvers  performed  are  beta-pointing  and 
coordinated  turn.  In  each  case  the  maximum  possible  maneuver 
without  exceed!"  -ontrol  input  limits  is  performed. 

Scope 

Lastly,  *  I i s  sis  incorporates  improvements  made  to 

MULTI.  Because  o t  the  large  number  of  maneuvers  and  flight 
conditions,  the  models  used  do  not  include  sensors,  actuators, 
lags,  or  control  input  rate  limitations.  However,  these  im¬ 
provements  can  easily  be  added  by  using  MULTI  to  adjust  the 
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controller  to  accomodate  these  refinements.  Since  the  ability 
to  determine  easily  the  input  rates  was  not  added  to  MULTI 
until  after  the  design  work  was  completed,  the  controller  and 
maximum  inputs  were  not  adjusted.  The  rates  are  presented, 
though,  and  where  they  exceed  allowable  maximums,  can  be 
decreased  by  reducing  the  maximum  steady-state  input  or 
lengthening  the  time  the  aircraft  has  to  perform  the  maneuver. 
This  is  easily  done  and  does  not  require  retuning  the  control¬ 
ler.  Because  of  limited  time,  no  attempt  was  made  to  develop 
a  robust  longitudinal  controller,  but  the  data  necessary  to  do 
so  are  included. 

Assumptions 

Grumman  developed  a  nonlinear  simulation  program  for  the 
X-29A.  Linearized  data  in  state-space  form  at  five  flight 
conditions  is  used.  Specifically,  "data  set  #3",  first 
promulgated  in  August  1983  is  used. 

The  controllers  synthesized  in  this  thesis  are  developed 
using  the  latest  MULTI  version  released  in  August  1983.  MULTI 
is  written  in  FORTRAN  and  is  run  on  the  CDC-6600  computer 
system  at  AFIT.  Calcomp  plotters  are  used  to  produce  the 
plots  in  this  thesis,  although  plots  at  the  terminal  are  also 
available.  Note  that  MULTI  is  an  interactive  design  tool, 
allowing  the  engineer  many  design  iterations  in  a  short  span 
of  time. 

Sequence  of  Presentation 

Chapter  II  briefly  describes  the  theory  used  to  develop 
MULTI  and  provides  a  list  of  references  for  the  serious  inte- 
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rested  industrious  reader.  Chapter  III  describes  the  X-29A. 
Chapter  IV  details  the  design  process  for  the  longitudinal 
controllers  and  presents  results.  Chapter  V  details  the 
design  process  for  the  lateral  controllers  and  presents 
results.  In  addition  a  single  robust  lateral  controller 
requiring  only  three  gains  is  obtained  for  all  aaneuvers  at 
all  five  flight  conditions.  Chapter  VI  summarises  the  results 
achieved  in  the  thesis,  gives  recommendations  for  future  work, 
and  details  a  class  of  future  possible  additions  to  MULTI  to 
improve  results  and  greatly  improve  productivity. 
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CHAPTER  II 


MULTIVARIABLE  CONTROL  LAWS  USING 


SINGULAR  PERTURBATION  METHODS 


Iotroduc  tton 

The  designs  presented  in  Chapters  IV  and  V  are  obtained 
using  the  interactive  computer-aided-design  (CAD)  tool 
MULTI.  MULTI  implements  the  multivariable  control  laws 
developed  by  Professor  Brian  Porter  of  the  University  of 
Salford,  England.  This  chapter  gives  a  brief  introduction 
to  two  of  the  design  methods  used  to  obtain  high  performance 
tracking  systems  for  the  X-29A. 

It  should  be  noted  that  the  continuous- time  methods  have 
a  corresponding  discrete-time  method  with  no  loss  of 
generality.  It  is  advantageous  to  examine  the  designs  in 
the  continuous- time  method  because  it  is  conceptually  easier 
to  consider  pole  and  transmission  zero  locations  in  the  s- 
plane  than  the  z-plane.  Before  accepting  the  discrete 
design,  sampling  effects  and  time  delays  must  be  considered 
(Ref  1:58-77). 

The  state  and  output  equations  are  described  by  first- 
order  linear  differential  equations  written  in  state-space 
format : 


Ax  +  Bu 


(2-1) 

(2-2) 


where 


■  the  continuous- time  plant  matrix  (n  x  n) 


B  •  the  continuous- time  control  matrix  (n  x  o) 

C  »  the  continuous- time  output  matrix  (p  x  n) 
x  ■  state  vector  (n  x  1) 
y  ■  output  vector  (p  x  1) 
u  ■  control  input  vector  (m  x  1) 

Note  that  Porter's  method  does  not  allow  a  feedforward 
matrix  (D).  For  systems  containing  a  D  matrix,  a  satisfac¬ 
tory  equivalent  representation  may  be  obtained  by  redefining 
the  control  inputs  as  states.  This  imbeds  the  D  matrix  into 
the  C  matrix  and  causes  the  actuator  inputs  to  become  the 
new  control  Inputs.  This  transformation  is  often  necessary 
for  systems  that  have  accelerations  as  outputs. 

The  proportional  plus  integral  (PI)  cascade  compensator 
performs  both  tracking  and  disturbance  rejection.  The 
results  in  this  thesis  are  obtained  using  a  discrete- time  PI 
implementation  in  MULTI.  The  equivalent  digital  equations 
are : 

;T 

exp{AT}  *  Bdt  *  u^(kT)  (2-3) 
*  0 

jr(kT]  -  C  *  x(kt)  (2-4) 

where 

expfAT}  "  sampled  data  plant  matrix 

C° 

J  exp-fAT^.  *  Bdt  •  sampled  data  control  matrix 


C  ■  sampled  data  output  matrix 
T*  sampling  period 

Since  u(kT)  is  piecewise  constant  over  the  sample  period,  it 
may  be  moved  outside  of  the  integral  (Ref  2:7-24).  Figure 


a 


2-1  shows  the  PI  controller.  MULTI  Implements  two  samplers, 


one  In  the  feedback  loop  from  the  plant,  and  the  other  in 


the  commanded  input.  The  single  sampler  shown  in  Figure  2-1 


is  equivalent  to  the  other  two  implemented  in  MULTI. 


The  A,  B,  and  C  matrices  are  partitioned  in  the 


following  manner: 


A  '  A 
111  12 


1  1 
+ 


(2-5) 


A  1  A 
21  I  22 


[  C  J  C 
1  2 


(2-6) 


where  A  is  a  matrix  of  order  ((n-m)  x  (n-m)) 
11 


where  A  ,  C  ,  and  B  are  matrices  of  order  (m  x  m) 
22  2  2 


(Ref  3:39). 


A  restriction  of  Porter's  method  is  that  the  number  of 


inputs  must  be  equal  to  the  number  of  outputs.  Control 


inputs  for  the  X-29A  aircraft  model  refer  to  canards, 


symmetrical  flaperons,  s trake - f laperons ,  thrust,  rudder,  and 


differential  flaperons.  Outputs  for  the  X-29A  are  selected 


s  ta  tes  . 


Two  design  methods,  Regular  and  Irregular,  are 


presented.  The  Regular  design  method  is  applied  to  the 


longitudinal  state-space  equations  while  the  Irregular  de- 


»**  w* 


v  • 


-*  v;  ?  7. 7»;  i 


>  ii  _ 


sign  method  is  applied  to  the  lateral  state-space  equations. 

A  major  advantage  of  Porter's  method  is  that  output 
feedback  is  used.  Consequently,  not  all  the  states  must  be 
known,  just  those  selected  in  an  Irregular  design. 

Known  Plant-  Regular  Design 

A  Regular  system  has  the  product  [CB]  equal  to  the 
number  of  outputs.  Thus  the  system  has  a  first  Markov 
parameter  of  maximal  rank  and  minimum  phase.  [CB]  must  have 
full  rank  (be  nonsingular)  so  that  it  is  be  invertible  (Ref 
3:7)  . 

A  requirement  of  Porter's  method  is  that  the  system  be 
controllable  (A,  B  a  controllable  pair)  and  observable  (A,  C 
an  observable  pair).  In  addition,  the  matrix 


(2-7) 


must  have  full  rank  to  assure  controllability  using  the  PI 
controller. 

The  control  law  for  the  digital  controller  is  given  by: 

u(kT)  -  (1/T)[K  i(kT)  +  K  £(kT)  ]  (2-8) 

0  1 

where 


K  is  the  proportional  gain  matrix  for  the  error  signal 
0 

je(kT) 

K  is  the  gain  matrix  for  the  backward  difference  of 
1 

the  error  signal,  which  is  designated  as  £(kT) . 
Equivalently,  in  the  continuous  case: 


10 


'  v’  OV  O  ■  *  * '  O  .  •  .  * 
lV.iV.lV jlVaVlV  fcV  •VVU.eV’e. 


where 


g  is  the  forward  path  gain 

K  is  the  proportional  gain  for  the  error  signal 
0 

1C  is  the  gain  for  the  integral  for  the  error  signal. 

1 

For  the  discrete  case: 

e(kT)  -  v(kT)  -  ^(kT)  (2-10) 

£[(k+l)T]  -  £( kT)  +  T£(kT)  (2-11) 

The  controller  matrices  are  as  follows: 

1 

K  -C^( C  )  Y  (2-12) 

0  2  2  == 

1 

K  -  (C  B  )  Y  (2-13) 

1  2  2=* 

where 

y  -  diagonal  matrix  {O’ »  O’ . O’}1 

12  p 

e 

Note  that  the  cascade  vector  integrator  drives  the 
steady  state  error  vector  to  zero  for  a  constant  command 
input  vector.  The  state  derivative  is  driven  to  zero  for 
all  cases  except  where  an  output  is  the  derivative  of  one  of 
the  states.  For  example,  in  the  longitudinal  mode  if  pitch 
rate  q  is  commanded  as  a  step  then  the  state  derivative  of 
pitch  angle  will  reach  the  commanded  value.  This  causes  a 
transmission  zero  at  the  origin.  For  this  case,  if  an 
Irregular  design  is  used,  the  state  derivative  must  not  be 
chosen  to  be  fed  back. 

When  Porter's  method  is  applied,  it  may  be  desirable  but 


not  necessary  for  the  B  matrix  in  equation  (2-4)  to  be  in 


the  form  |  0 

1  B 


MULTI  accepts  either  form. 


For  the  continuous  case,  as  the  gain  factor  g  increases 
to  infinity  in  the  limit,  the  system  transfer  matrix  given 
as  follows: 


G ( s )  -  C(sl  -  A)  B 


(2-14) 


assumes  an  asymptotic  form 


r<\>  -  r<\>  ♦  r<\> 


(2-15) 


where 


3T*  (\)  is  the  "slow"  transfer  function  matrix 
A 

T1  (\)  is  the  "fast"  transfer  function  matrix 

The  roots  of  these  transfer  functions  are  grouped  into 

three  categories  as  Z  ,  Z  ,  and  Z  .  The  slow  modes  are 

12  3 

the  roots  in  Z  and  Z  .  As  the  gain  factor  increases  in  the 
1  2 

limit,  the  poles  of  Z  become  uncontrollable  while  the  poles 

1 

of  Z  become  unobservable.  As  this  occurs,  the  fast  modes 
2 

corresponding  to  Z  increasingly  dominate  the  system  re- 

3 

sponse.  Expressions  for  the  asymptotic  modes  are  given  in 
Table  2-1  (Ref  1:69-70).  Therefore,  as  the  sampling  period 
T  approaches  zero,  the  closed-loop  transfer  function 
becomes : 


<r 


£ 


r<X>  "  p  <X>  “  <Xi  "I  +  c  B  K  )  C  £  K 

p  p  220  220 


(2-16) 


Thus  if  jC  is  chosen  as  follows: 
0 


_C  ji  K  «  diagonal  ma  tr  ix  -T  (J  ,  (J  ,....,  \  »Y 

220  12  n  J  — 


then 


F(X)“  diagonal  ma  tr  ix  /  (\-l+^- )  ,  .  .  ,  j*/ (^-1+ QT)}  (2-17) 

1  1  p  p 


The  fast  mode  roots  and  the  transmission  zeros  must  lie 

within  the  unit  circle  in  the  z-plane  (discrete)  for  the 

system  to  be  stable.  The  transmission  zeros  are  a  subset  of 

the  slow  modes  and  can  be  calculated  from  the  equation  in 

Table  2-1  since  Z  *  Z  . 

t  2 

From  equation  (2-17)  it  can  be  seen  that  decoupling  is 
achieved  as  the  sampling  period  goes  to  zero.  If  the 
transmission  zeros  are  not  in  the  stable  region,  then  the 
output  matrix  C  must  be  changed.  Alternately,  a  lower  limit 
on  the  sample  period  may  be  acceptable.  If  not,  the  differ¬ 
ent  set  of  outputs  requires  that  a  new  design  be  attempted. 
Known  Plant-  Irregular  Design 

When  the  product  [ GB ]  is  rank  deficient,  then  an 
Irregular  design  must  be  used.  Extra  output  measurements 
are  required.  The  new  feedback  vector  is  given  as  follows: 


w(  t) 


[*■* 


MA  C  +  1U  1  x  (O 

11  2  12  1 

x  (t) 
2 


(2-18) 


where 


f!  is  of  order  (p  x  (n-p)) 
is  of  order  (p  x  (n-p)) 

1 

F^  is  of  order  (p  x  p) 

2 

The  system  is  shown  in  Figures  '2-2  and  2-3.  The  ele¬ 
ments  of  the  measurement  matrix  ]1  are  selected  so  that  the 
matrix  [rB]  has  full  rank.  Thus  the  F  matrix  replaces  the  C 
matrix  in  the  design  method.  The  control  law  equation  (2-8) 
is  unchanged,  but  the  error  vector  is  redefined  as: 

^(kT)  «  v(kT)  -  w(kT)  (2-20) 

Since  the  command  input  v(t)  is  a  constant  vector 

•  • 
lim  x  (t)  “  lim  [  A  x  (t)  +  A  x  (t)]  -  0 
t-vOO  1  t-*oo  11  1  12  2 

then 

lim  e(kT)»lim  {v(kT)  -^(kT) -M  [  A  x  (kt)  +  A  x  (kt)]} 
k-voo  k  -e»oe  11  1  12  2 

-lim  £v(kT)  -  ^(kT)} 

k-^ao 

-  0  (2-21) 
Thus  the  error  vector  is  driven  to  zero  in  the  steady  state 
for  a  constant  command  input  vector. 


For  an  Irregular  design,  the  closed-loop  transfer  func¬ 
tion  matrix  of  equation  (2-15)  approaches  an  asymptotic 


value  as  the  sampling  period  goes  to  zero  (Ref  3:121-122): 


$ 


.  _!  _1  _1  _1 
P<\)-(C  -C  F  F  )(\l  -I  -TA  +  TA  F  F  )  (TA  F  ) 
1221  n-p  n-p  11  12  2  1  12  2 

1  1 

r<\>  -  £  l  <Ai  - 1  *  i  B  K  )  F  1  K  (2-22) 

22  p  p  220  220 

The  transmission  zeros  are  determined  as  for  a  Regular 
system,  except  that  the  F  matrix  is  substituted  for  the  C 
matrix.  The  resulting  transmission  zeros  must  be  in  the 
open  left  half  s-plane  (continuous  case)  or  within  the  unit 
circle  in  the  z-plane  (discrete  case). 

For  an  Irregular  design,  as  the  gain  factor  goes  to 
infinity,  the  slow  modes  corresponding  to  the  poles  Z 

1 

become  uncontrollable.  The  fast  modes  corresponding  to  the 

poles  of  Z  remain  controllable  and  observable  as  the  gain 
3 

1/T  goes  to  infinity.  However,  unlike  a  Regular  design,  the 

slow  modes  corresponding  to  the  poles  of  Z  remain  both 

2 

observable  and  controllable.  Thus  the  slow  transfer  func¬ 
tion  matrix  does  not  approach  zero  as  the  sampling  period  is 
reduced  to  zero.  Therefore  the  slow  and  fast  transfer 
matrices  given  in  (2-22)  must  be  diagonal  for  decoupling  to 
take  place  in  the  limit. 

The  measurement  matrix  M  is  chosen  such  that: 

1.  M  is  as  sparse  as  possible. 
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2.  P  B  has  full  rank. 

2  2 

3.  The  resulting  transmission  zeros  must  be  in  the 
stable  region.  Otherwise,  there  is  a  lower  limit  on  the 
sampling  period. 

4.  Decoupling  of  outputs  in  the  limit  as  the  sampling 
frequency  goes  to  infinity  is  guaranteed  only  if  the  product 

_1 

C  F  is  a  diagonal  matrix. 

2  2 

5.  The  system  has  closed-loop  roots  in  the  stable 
region . 

The  controller  matrices  must  be  chosen  such  that: 

1 

K  m  01  K  -01  <F  B  )“  diag{0*  ,g*  /r  > 

0  1  2  2  12  p 


(2-23) 


CHAPTER  III 
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THE  X-29A  DEMONSTRATOR 


Aerodynamic  s 

The  X-29A  demonstrator  aircraft  integrates  several  differ¬ 
ent  advanced  technologies  of  aircraft  design  (see  Figure  3-1). 
This  chapter  discusses  f orward- swep t  wings,  close-coupled 
canards,  strake-f laps ,  variable  camber,  fly-by-wire,  relaxed 
stability,  and  the  interactions  between  them.  Although  no 
knowledge  of  the  actual  aircraft  design  is  needed  to  apply 
Porter's  method,  it  is  often  helpful  to  gain  some  insight  into 
how  the  aircraft  "should"  behave. 

In  general,  swept  wings,  fore  and  aft,  help  to  move  the 
aircraft  through  the  air  more  efficiently  by  reducing  compres¬ 
sibility  drag.  Drag  slows  the  aircraft  by  transferring  energy 
from  the  aircraft  to  surrounding  air.  As  a  wing  travels 
through  air,  it  compresses  the  air  in  front  of  it,  creating 
compressibility  drag  (Ref  4).  As  the  aircraft  increases 
speed,  compressibility  drag  increases.  By  sweeping  the  wings 
at  an  angle  nonperpendicular  to  the  forward  velocity  of  the 
aircraft,  the  velocity  component  of  the  airstream  perpendicu¬ 
lar  to  the  leading  edge  of  the  aircraft  is  less  than  that  of 
the  aircraft  itself.  This  causes  the  air  to  move  along  the 
length  of  wing,  instead  of  just  from  leading  edge  to  trailing 
edge . 

As  the  aircraft  enters  supersonic  speeds,  wave  drag  also 
occurs.  In  the  subsonic  region,  the  pressure  wave  created  by 


the  wing  pushes  air  out  of  the  way 


But  in  the  supersonic 
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region  the  pressure  wave,  which  only  travels  at  the  speed  of 
sound,  is  overtaken  by  the  aircraft  which  is  traveling  faster 
than  the  speed  of  sound.  These  differences  in  speed,  and 
hence  pressure,  cause  a  sonic  boom.  This  sonic  boom  is  caused 
by  a  transfer  of  energy  from  the  aircraft  to  the  surrounding 
airstream.  This  wave  drag  can  be  minimized  by  presenting  a 
smooth  minimized  cross-sectional  area  of  the  aircraft  to  the 
oncoming  airstream.  Forward-swept  wings  have  a  much  smoother 
cross-sectional  area  than  aft-swept  wings,  and  hence  have  less 
wave  drag. 

Forward-swept  wings  are  placed  further  back  along  the 
fuselage,  further  away  from  the  center  of  gravity,  than  aft- 
swept  wings.  In  order  to  minimize  pitching  moment,  the  pay- 
load  of  the  aircraft  is  placed  at  the  center  of  gravity  of  the 
aircraft.  Aft-swept  wings,  being  close  to  the  center  of 
gravity  of  the  aircraft,  present  a  dilemma  for  the  designer. 
If  the  fuselage  is  "pinched- in"  to  present  a  smooth  cross- 
sectional  area  ,  then  payload  is  reduced.  If  payload  capabil¬ 
ity  is  not  reduced,  then  the  resulting  uneven  cross-sectional 
area  increases  the  wave  drag  at  supersonic  speeds. 

Forward-swept  wings  have  an  important  advantage  over  aft- 
swept  wings  in  " s ta 1 1 -progress  ion" .  In  order  to  generate  more 
lift  for  the  aircraft  the  angle  of  attack  (between  the  wing 
and  the  airstream)  is  increased.  As  this  happens,  the  air 
passing  across  the  top  of  the  wing  speeds  up.  By  Bernoulli's 
principle,  this  difference  in  speeds  causes  a  pressure  differ¬ 
ential  that  causes  lift.  If,  however,  the  angle  of  attack 


becomes  too  great  then  the  airstream  can  not  follow  the  con¬ 
tour  of  the  wing,  causing  stall.  Stall  progression  occurs  in 
the  opposite  direction  of  airflow.  Thus  stall  progression 
moves  from  trailing  edge  to  leading  edge  of  the  wing  and  along 
the  wing  from  the  aft-most  part  of  the  wing  to  the  fore-most 
part  of  the  wing.  In  aft-swept  wings,  the  aft-most  part  of 
the  wing  is  at  its  tips,  where  the  ailerons  are  located.  This 
means  that  the  ailerons  are  no  longer  effective  and  control  is 
lost.  In  f orward- swep t  wings  the  aft-most  part  of  the  wing  is 
at  its  root,  or  the  part  of  the  wing  that  is  attached  to  the 
fuselage.  Thus  stall  first  occurs  away  from  the  ailerons. 
Moving  the  ailerons  to  the  root  for  aft-swept  wings  is  not 
feasible  since  this  would  reduce  the  roll  moment  by  decreasing 
the  lever  arm.  These  advantages  allow  the  X-29A  to  fly  at 
angles  of  attack  of  up  to  70  degrees.  Present  aircraft  with¬ 
out  vectored  thrust  cannot  fly  at  angles  of  attack  greater 
than  30  degrees  or  stall  occurs. 

Close-Coupled  Canards 

Close-coupled  canards,  located  in  front  of  the  forward- 
swept  wings,  can  be  used  to  channel  air  to  the  root  of  the 
wing  in  order  to  resist  stall.  This  further  increases  the 
capability  of  the  X-29A  to  fly  at  greater  angles  of  attack 
before  stall  is  induced. 

In  statically  stable  aircraft,  the  lift  vector  is  located 
aft  of  the  center  of  gravity.  This  means  that  the  aircraft 
tends  to  have  a  nose-down  pitching  moment.  In  conventional 
aft-swept  wing  aircraft,  this  is  countered  by  placing  a  rear 
horizontal  stabilizer  to  provide  a  downward  force  to  counter 
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wing  lift.  While  this  cancels  out  the  pitching  moments  about 
the  center  of  gravity,  it  also  reduces  the  total  lift  of  the 
aircraft.  In  forward- swep t  wing  aircraft  the  horizontal  sta¬ 
bilizers,  or  canards,  are  placed  in  front  of  the  center  of 
gravity.  Thus  while  it  still  cancels  out  the  pitching  moment 
about  the  center  of  gravity,  it  now  adds  to  total  lift.  With 
the  lift  vectors  from  the  wings  and  canards  on  opposite  sides 
of  the  center  of  gravity,  the  center  of  gravity  can  be  moved 
around,  permitting  much  greater  flexibility  in  payload  distri¬ 
bution  . 

Strake-Flaps 

Strakes,  horizontal  surfaces  attached  to  the  fuselage  from 
the  wings  to  the  tail,  have  flaps  attached  at  the  tail  to 
provide  a  nose-down  pitching  moment  at  high  angles  of  attack. 
(Ref  5).  The  strakes  also  aid  in  increasing  pitch  damping  by 
resisting  movement  as  the  rear  of  the  fuselage  moves  up  and 
down  (Ref  4).  The  strakes  can  be  moved  in  concert  with  the 
close-coupled  canards  and  symmetrical  flaps  on  the  wings  to 
obtain  optimal  air  flow  as  the  airstream  passes  from  front  to 
back  of  the  aircraft. 

Wing  Composition  and  Shape 

Past  applications  of  forward- swept  wings  have  been  lim¬ 
ited  to  small  sweep  angles  because  of  structural  divergence  of 
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the  wings.  As  the  angle  of  attack  of  the  wing  increases,  the 
wing  is  twisted.  In  aft-swept  aircraft,  the  leading  edge  of 
the  wing  twists  downwards,  thus  relieving  the  load  by  decreas¬ 
ing  the  angle  of  attack.  In  forward - swep t  wing  aircraft, 
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however,  the  wing  is  twisted  upwards,  increasing  the  angle  of 


attack.  This  leads  to  even  greater  twisting  of  the  wing  until 
the  structure  cannot  bear  the  load  and  could  be  torn  away. 

Grumman  has  been  able  to  overcome  this  problem  by  using  a 
process  known  as  "aeroelas  tic  tailoring".  In  aeroelastic 
tailoring,  composite  materials,  such  as  boron  fiber  and  graph¬ 
ite  in  a  plastic  matrix,  are  used  to  produce  a  nonmetallic 
wing  that  is  307.  lighter  than  a  metal  wing  of  the  same 
strength.  In  addition,  the  wing  structure  is  optimized  by 
varying  the  direction  and  thickness  of  the  layers  of  bonded 
fiber  composite  material.  This  layering  effectively  distrib¬ 
utes  the  twisting  load,  which  is  greatest  at  the  wingtips, 
over  the  entire  wing. 

The  X-29A  uses  supercritical  wings  to  increase  its  criti¬ 
cal  Mach  number.  The  critical  Mach  number  is  the  speed  at 
which  the  aircraft  can  fly  before  the  air  hitting  the  front  of 
its  wings  increases  to  the  speed  of  sound.  As  the  air  passes 
over  the  top  of  the  wing  it  slows  down,  producing  a  shock 
wave.  This  phenomenon,  known  as  the  Bernoulli  effect,  in¬ 
creases  drag  and  decreases  lift.  Supercritical  wings  have  a 
greater  leading  edge  radius  of  curvature  and  a  flatter  surface 
top.  This  causes  the  air  to  reach  supersonic  speeds  earlier 
and  sustains  the  supersonic  speed  further  along  the  wing. 
This  pushes  the  shock  wave  to  the  trailing  edge  of  the  wing, 
where  drag  is  much  less  and  lift  is  not  decreased  as  much. 
Variable  Camber 

The  X-29A  uses  trailing  edge  surfaces  to  vary  the  camber 
of  the  wing  to  provide  the  optimal  wing  shape  over  differing 
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flight  conditions.  The  surfaces  are  composed  of  three  parts 
along  the  wing.  The  outermost  part,  called  the  outboard 
flaperon,  is  the  equivalent  of  an  aileron  which  is  used  to 
provide  roll  control.  The  inboard  flaperon,  composed  of  two 
parts,  is  similar  to  a  traditional  flap  to  provide  additional 
lift  during  landing.  The  trailing  edge  control  surfaces  are 
separated  into  three  parts  to  combat  aerodynamic  flutter. 
The  inboard  flaperon  uses  one  actuator  to  move  both  trailing 
edge  surfaces.  This  is  done  because  of  weight  and  space 
constraints.  Both  actuators  are  moved  together.  Independent 
movement  of  the  inboard  and  outboard  flaperons  is  not  possible 
because  of  nonlinearities  in  the  actuators.  Thus  all  three 
parts  move  together  as  symmetrical  flaperons  to  provide 
variable  camber  and  as  differential  flaperons  to  perform  la¬ 
teral  maneuvers.  This  results  in  increased  maneuverability 
and  reduced  drag  throughout  the  flight  envelope  (See  Figure  3- 
2).  An  integrated  flight  control  computer  is  used  to  automa¬ 
tically  control  optimum  wing  shape  during  steady  level  flight. 
To  provide  variable  camber,  the  flaperons  are  moved  as  a 
function  of  canard  displacement.  If  the  pilot  wishes  to  use 
the  symmetrical  flaperons  to  perform  a  maneuver,  he  must 
switch  to  manual  camber  control. 

Common  Components 

In  order  to  minimize  cost,  5  67.  by  weight  of  the  X-29A 
consists  of  components  which  are  common  to  existing  aircraft. 
The  forward  fuselage  and  cockpit  area  are  from  the  F-5A  fight¬ 
er;  the  main  landing  gear,  control  surface  actuators,  auxilia- 
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ry  power  system,  and  engine  accessory  drive  are  from  the  F-16; 
and  the  engine  is  a  GE  F404  turbofan  that  is  from  the  F-18. 


Fly-by-Wire 

Past  aircraft  have  used  mechanical  inputs  from  the  pilot 
to  power  hydraulic  actuators  which  move  the  control  surfaces. 
The  X-29A  uses  electric  signals  from  the  pilot,  called  "fly¬ 
by-wire",  to  activate  the  hydraulic  actuators  that  move  the 
control  surfaces.  This  advanced  technology  reduces  weight  and 
space  and  may  provide  increased  reliability  through  multiple 
redundancy  and  distributed  processing.  Future  aircraft  might 
possibly  be  " a  1 1 - e 1 ec t r ic"  ,  in  which  not  only  the  signals  from 
the  pilot  will  be  electric,  but  the  actuators  that  move  the 
control  surfaces  will  be  electrically  powered  as  well.  This 
has  the  additional  advantage  of  greatly  reducing  the  "dead- 
zone",  the  lag  in  time  between  when  the  inputs  are  commanded 
and  when  the  control  surfaces  actually  start  to  move.  Future 
technology  might  also  provide  fiber-optic  links  to  the  hy¬ 
draulic  system.  In  addition  to  further  reducing  weight  and 
space,  such  circuits  are  less  vulnerable  than  electric  cir¬ 
cuits  are  to  "EMP"  (electromagnetic  pulse). 

Relaxed  Stability 

The  X-29A  is  statically  unstable  at  subsonic  speeds  (see 
Figure  3-3).  Since  an  uncompensated  perturbation  in  the  lon¬ 
gitudinal  plane  would  cause  the  X-29A  to  go  unstable  in  a 
fraction  of  a  second,  digital  flight  control  computers  must 
provide  "artificial  stability"  for  the  aircraft.  This  has  the 


advantage  that  smaller  and  lighter  flight  control  surfaces  are 


needed.  But  this  also  has  the  disadvantage  that  if  the  compu¬ 
ter  hardware  should  fail,  or  if  the  software  of  the  controller 
is  faulty,  then  the  pilot  will  be  unable  to  recover  the  air¬ 
craft.  Hardware  failure  risk  is  minimized  by  providing  a 
backup  power  supply  and  redundant  digital  computers.  Software 
failure  risk  is  guarded  against  by  providing  a  backup  "analog 
reversion-mode"  system  that  enables  the  pilot  to  recover  the 
aircraft.  But  this  adds  weight  and  cost  to  the  aircraft. 

At  supersonic  speeds  the  center  of  pressure  (point  at 
which  the  lift  vector  of  the  aircraft  is  applied)  moves  back 
towards  the  center  of  gravity  of  the  aircraft.  The  aircraft 
is  then  neutrally  stable  and  smaller  control  surface  deflec¬ 
tions  are  needed  to  stabilize  and  maneuver  the  aircraft.  This 
decreases  drag  since  every  time  a  control  surface  is  moved,  a 
transfer  in  energy  from  the  aircraft  to  the  surrounding  air- 
s  tream  occurs . 


Summary 

Since  the  X-29A's  aerodynamics  enable  it  to  move  more 
efficiently  through  air,  less  fuel  is  required  to  fly  it  than 
a  conventional  aircraft  of  similar  size.  Smaller  weight  and 
size  also  permit  a  reduction  in  cost  of  building  the  X-29A. 
Thus  fixed  and  marginal  costs  are  kept  to  a  minimum.  The 
" 1 i f t- to- trim"  canards  offer  the  advantage  of  freedom  of  pay¬ 


load  distribution.  The  reduction  in  size  offers  increased 


survivability,  since  smaller  aircraft  are  harder  to  detect. 
Relaxed  stability  permits  greater  maneuverability,  and  hence 
greater  survivability  and  mission  effectiveness.  Digital 


flight  control  systems  may  also  offer  better  performance 
battle  damage  is  sustained  by  permitting  reconfiguration 
the  control  laws. 

Unfortunately,  this  increase  in  performance  must  be 
for  with  increased  complexity  in  design. 


CHAPTER  IV 


t  m  .  * 

LONGITUDINAL  CONTROLLER  DESIGNS 

Introduc tl on 

This  chapter  presents  the  results  of  the  longitudinal 
controller  designs  for  the  five  flight  conditions  for  each  of 
six  maneuvers.  A  sample  flight  condition,  1.2  Mach  at  15  000 
feet,  is  presented.  Although  each  flight  condition  is  discus¬ 
sed,  supporting  state-space  matrices,  design  parameters,  com¬ 
manded  inputs,  controller  matrices,  figures  of  merit,  and 
plots  for  the  other  four  flight  conditions  are  included  in  the 


appendices  . 

Grumman  Aerospace  Corporation  developed  a  nonlinear  simu¬ 
lation  computer  program  for  the  X-29A  demonstrator  aircraft  at 
different  flight  conditions.  From  this  program  linearized 
data  in  state-space  form  at  five  flight  conditions  is  ob¬ 
tained.  Many  design  methods,  including  Porter's  method, 
assume  linear  behavior  in  order  to  achieve  mathematically 
tractable  results.  While  the  simplification  of  linearized 
performance  of  the  X-29A  for  small  perturbations  about  the 
chosen  flight  condition  is  a  good  one,  it  becomes  less  valid 
the  greater  the  perturbations  become.  The  usual  remedy  is  to 
design  controllers  at  many  flight  conditions,  and  then  use 
gain-scheduling  between  the  discrete  design  points.  This 
approach  requires  that  the  design  method  lend  itself  readily 
to  development  of  many  controllers  with  minimum  effort. 

In  order  to  further  simplify  the  model,  the  longitudinal 
and  lateral  controllers  are  developed  separately.  This  is 
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possible  in  the  X-29A  because  no  cross-coupling  occurs.  If 
this  were  not  the  case,  though,  this  would  be  a  good  first 
approach.  After  separate  designs  are  obtained,  the  complete 
long! tudina 1 / la tera 1  controller  could  be  retuned  to  accomodate 
cross-coupling.  If  the  cross-coupling  is  too  great,  however, 
then  the  controller  would  have  to  be  designed  using  the  full 
model.  This  increases  complexity,  though,  and  does  not  allow 
the  advantage  of  separate  design  teams  for  the  longitudinal 
and  lateral  controllers  (Ref  6).  Decoupling  also  permits 
separate  computers  for  the  longitudinal  and  lateral  control¬ 
lers. 

In  order  to  further  simplify  the  model,  a  reduced  number 
of  states  is  used  in  the  Grumman  model.  The  states  used  are 
q,  pitch  rate;  theta,  pitch  angle;  alpha,  angle-of-attack;  and 
VTOT,  total  velocity.  This  controller  design  further  assumes 
that  VTOT  is  approximated  by  0,  forward  velocity.  Note  that 
since  the  model  used  assumes  steady  level  flight  as  the  ini¬ 
tial  condition,  U  can  be  replaced  with  u,  the  change  in 
forward  velocity. 

Symmetrical  flaperons,  thrust,  canards,  and  s trake- f laper- 
ons  (strakes)  are  the  available  longitudinal  inputs.  Since 
the  X-29A  uses  variable  camber,  this  initial  design  attempt 
uses  only  thrust,  canards,  and  strakes  as  control  inputs. 
This  allows  a  separate  controller  to  be  designed  to  adjust  the 
symmetrical  flaperons  to  provide  variable  camber. 

Porter's  method  requires  that  the  number  of  inputs  be 
equal  to  the  number  of  outputs.  It  is  difficult  to  command 
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four  separate  outputs  that  do  not  conflict  with  one  another. 
If  all  four  inputs  are  used,  then  all  four  states  can  be 
controlled  as  outputs.  Since  commanding  q  and  theta  is  redun¬ 
dant,  a  transmission  zero  at  the  origin  occurs.  This  signi¬ 
fies  that  redundant  inputs  are  present.  Even  though  the 
outputs  commanded  are  redundant,  the  required  maneuver  can  now 
be  performed  with  four  control  surfaces. 

Thus  all  of  the  longitudinal  maneuvers  performed  under¬ 
state  the  actual  longitudinal  performance  capabilities  of  the 
X-29A.  Future  designs  should  use  two  separate  longitudinal 
controllers  for  the  X-29A.  One  design  would  treat  the  ca¬ 
nards,  symmetrical  flaperons,  and  strakes  as  one  control  input 
in  order  to  minimize  drag.  This  design  would  more  properly 
model  aerodynamic  effects,  since  changes  in  the  canard  affect 
the  airflow  over  the  inboard  flaperons.  Changes  in  the  in¬ 
board  flaperons,  which  occur  with  active  variable  camber,  in 
turn  affect  the  airflow  passing  over  the  strakes.  Thus  move¬ 
ment  of  the  flaperons  would  be  a  function  of  canard  movement, 
while  movement  of  the  strakes  would  be  a  function  of  flaperons 
or  canards.  This  controller  would  be  used  all  of  the  time  and 
would  be  the  only  controller  used  when  the  X-29A  is  cruising. 
Since  the  number  of  independent  control  inputs  is  reduced  to 
only  two  (canards,  flaperons,  and  strakes;  thrust),  only  two 
outputs  could  be  controlled.  The  other  longitudinal  con¬ 
troller  would  treat  all  the  control  inputs  as  independent 
inputs,  permitting  maximum  performance.  This  controller  would 
enable  the  aircraft  to  achieve  maximum  maneuverability,  and 
would  be  activated  by  the  pilot  when  he  wished  to  perform  a 


maneuver,  or  when  an  automatic  set  of  pre-programmed  maneuvers 


was  initiated,  as  in  integrated  fire/flight  control.  The 
second  controller  would  take  precedence  over  the  first,  since 
this  controller  would  be  used  in  combat.  When  absolute  maxi¬ 
mum  maneuvers  were  not  being  used,  then  the  first  controller's 
inputs  would  be  added  to  the  second  controller's  inputs  in 
order  to  perform  the  desired  maneuvers  with  the  greatest 
efficiency . 

The  outputs  chosen  to  be  controlled  are  the  states  q,  u, 
and  alpha.  It  is  extremely  important  that  the  same  output  (C) 
matrix  be  used  for  all  of  the  maneuvers.  Otherwise,  the  pilot 
would  have  to  flip  a  switch  to  tell  the  aircraft  what  maneuver 
he  was  about  to  command.  In  addition,  different  output 
matrices  make  attainment  of  a  universal,  or  robust  controller, 
extremely  doubtful.  Lastly,  the  advent  of  digital  flight 
control  permits  automatic  maneuvering  while  in  combat.  This 
occurs  when  the  pilot  has  maneuvered  the  aircraft  within  close 
proximity  to  the  target,  and  gives  authority  to  the  controller 
to  take  over.  Different  output  matrices  can  be  accomodated  by 
a  digital  flight  control  system  that  uses  a  predetermined  set 
of  maneuvers,  but  now  the  problem  of  large  gain  transients  in 
switching  back  and  forth  to  different  controllers  must  be 
overcome . 

Description  o f  Maneuvers 

Controllers  are  designed  for  six  longitudinal  maneuvers  at 
each  of  the  five  flight  conditions.  The  initial  design  at¬ 
tempts  used  a  different  controller  for  each  maneuver,  for  a 


total  of  forty  controllers.  However,  it  is  noticed  that 
performance  is  only  slightly  degraded  by  using  the  same 
controller  for  three  of  the  maneuvers,  another  controller  for 
two  of  the  maneuvers,  and  a  separate  controller  for  the 
remaining  maneuver.  In  keeping  with  the  overall  goal  of 
obtaining  a  robust  controller,  this  reduces  the  total  number 
of  controllers  to  fifteen  longitudinal  controllers. 

Note  that  all  outputs  must  be  commanded  to  some  value;  no 
output  is  allowed  to  "float".  Since  the  X-29A  is  a  new  air¬ 
craft,  maximum  possible  maneuvers  are  commanded.  Note  that  in 
many  instances  transient  responses  are  degraded  in  order  to 
obtain  maximum  performance.  The  maximum  possible  maneuver  is 
determined  by  simply  increasing  the  commanded  inputs  until 
input  limits  are  met.  MULTI  allows  input  limits  to  be 
reached,  but  not  exceeded  during  a  simulation,  and  notifies 
the  designer.  This  frees  the  design  engineer  from  having  to 
determine  these  maximums  by  using  a  rigorous  a-priori  aerody¬ 
namic  analysis.  Since  the  option  to  easily  determine  control 
input  rates  was  not  added  to  MULTI  until  after  design  work  was 
completed,  control  input  rates  are  exceeded  for  some  maneu¬ 
vers.  Where  the  rates  are  exceeded,  the  maximum  commanded 
steady-state  input  need  only  be  reduced  or  the  maneuver  time 
lengthened,  but  the  controller  does  not  have  to  be  refined. 
In  order  to  facilitate  easy  comparison  of  aircraft  capabili¬ 
ties  at  different  flight  conditions,  each  longitudinal  maneu¬ 
ver  was  commanded  over  the  same  time  duration.  Different 


flight  conditions  performed  different  maneuvers  better  than 
others,  and  flight  condition  number  three,  at  1.2  Mach  at  15 


000  feet  is  a  good  representative.  One  of  the  additions  to 
MULTI  that  this  thesis  makes  is  to  incorporate  a  prefilter  of 


the  pilot's  commanded  inputs.  The  prefilter  allows  an  in¬ 
crease  of  approximately  three  percent  in  the  commanded  steady- 
state  value  of  the  maximum  maneuver.  This  is  traded  off 
against  a  degradation  in  settling  time  values,  however. 

The  six  longitudinal  maneuvers  designed  for  are  direct 
climb,  direct  liftl,  direct  lift2,  pitch-pointing,  vertical 
trans lationl ,  and  vertical  trans lation2 .  One  controller  is 
used  for  direct  climb,  direct  liftl,  and  direct  lift2  at  each 
flight  condition.  In  direct  climb,  q  is  commanded  as  a  pulse, 
u  is  ramped-up  and  held  to  a  steady-state  value,  and  alpha  is 
commanded  to  zero.  This  maneuver  allows  the  pilot  to  point  the 
nose  of  the  aircraft  at  a  target  and  increase  the  airspeed  in 
order  to  chase  it.  In  direct  liftl,  q  is  commanded  as  a  pulse 
of  0.7  seconds  duration,  while  u  and  alpha  are  commanded  to 
zero.  This  gives  an  indication  cf  how  well  the  aircraft  can 
change  its  pitch  in  a  short  time.  In  direct  lift2,  q  is 
commanded  as  a  pulse  of  3.2  seconds  duration,  while  again  u 
and  alpha  are  commanded  to  zero.  Note  that  since  alpha  is 
commanded  to  zero,  che  flight-path  angle  gamma  is  now  equal  to 
theta.  Theta  can  be  calculated  by  taking  the  integral  of  the 
commanded  q  input. 

A  single  controller  is  used  for  pitch-pointing  at  each 
flight  condition.  In  p i t ch - p o i n t ing ,  q  is  commanded  to  a 
value  that  produces  a  pitch  angle,  theta,  equal  to  the  com¬ 
manded  value  of  alpha.  U  is  commanded  to  zero.  This  maneuver 


is  of  greatest  benefit  to  aircraft  that  cannot  fly  at  large 
angles  of  attack.  Pitch-pointing,  also  called  fuselage¬ 
pointing,  allows  the  aircraft  to  increase  its  pitch  angle 
without  increasing  its  flight  path  angle,  gamma.  Gamma  is 
held  to  zero  since  theta  and  alpha  are  kept  equal. 

A  single  controller  is  used  for  vertical  translationl  and 
vertical  trans la t ion2 .  In  vertical  translationl,  q  and  u  are 
commanded  to  zero  while  alpha  is  commanded  as  a  pulse.  This 
maneuver  moves  the  aircraft  upwards  a  fixed  distance.  This 
maneuver  is  useful  for  refueling  and  positioning  the  aircraft 
within  a  formation.  In  vertical  trans la t ion2 ,  q  and  u  are 
again  commanded  to  zero,  but  this  time  alpha  is  commanded  to 
ramp-up  and  hold  a  steady-state  value.  This  maneuver  causes 
the  aircraft  to  rise  at  a  fixed  velocity  while  keeping  the 
pitch  angle  zero. 

Development  o  f  Aircraft  Model 

The  first  step  in  digital  flight  control  design  is  the 
development  of  the  aircraft  equations  of  motion  based  upon 
certain  assumptions  (Ref  7).  Such  assumptions  may  include: 

1.  The  X  and  Z  axes  lie  in  the  longitudinal  plane  with  the 
center  of  the  chosen  axis  system  at  the  center  of  gravity  of 
the  aircraft. 

2.  The  mass  of  the  aircraft  is  constant.  Although  much  fuel 
may  be  expended  (up  to  23.87.  by  weight  of  the  X-29A)  during  a 
mission,  the  amount  of  fuel  expended  during  a  16  second 
maneuver  can  be  safely  disregarded. 

3.  The  earth  is  an  inertial  reference.  For  an  aircraft  this 


is  a  good  assumption. 


4.  The  perturbations  from  equilibrium  (steady  level  flight) 
are  small. 

5.  Quasisteady  flow  occurs.  This  assumes  that  the  airstream 
surrounding  the  X-29A  changes  instantaneously  when  the 
aircraft  is  disturbed  from  equilibrium.  Although  at  higher 
Mach  numbers  (transonic  and  above)  compressibility  drag 
renders  this  a  less  satisfactory  assumption,  it  is  a  good 
first  approximation. 

6.  The  aircraft  is  a  rigid  body.  After  the  design  is 
completed,  a  bending  mode  analysis  can  be  completed  to  see  if 
this  is  a  valid  assumption. 

Next,  the  equations  of  motion  must  be  used  to  develop  a 
nonlinear  simulation  program.  Typically,  an  aerodynamist 
develops  the  equations  of  motion  and  aircraft  model  computa¬ 
tion  program  for  the  digital  flight  control  engineer.  Note 
that  a  limited  nonlinear  simulation  capability  has  been  added 
to  another  program  supporting  MULTI  (Ref  8).  In  this  thesis, 
the  linearized  data  received  from  NASA  are  assumed  to  be 
correct  since  the  equations  of  motion  are  not  made  available. 
After  transforming  the  linearized  data  containing  the  plant  (A 
matrix)  and  inputs  (B  matrix)  into  state-space  form,  the 
designer  must  decide  what  outputs  (C  matrix)  will  be  used. 
For  certain  outputs  a  feedforward  (D  matrix)  term  may  also 
have  to  be  added.  The  outputs  chosen  will  determine  whether  a 
Regular  or  Irregular  design  method  is  first  attempted.  A 
Regular  design  offers  three  advantages.  First,  a  Regular 
design  requires  fewer  gains  than  an  Irregular  design  does. 


Second,  a  Regular  design  offers  the  possibility  of  faster 
performance  since  some  slow  modes  become  uncontrollable  in  an 
Regular  design.  Lastly,  other  things  being  equal,  a  Regular 
design  is  easier  since  fewer  parameters  must  be  chosen. 

The  design  of  the  longitudinal  controllers  in  this  thesis 
uses  a  Regular  design  with  a  constant  C  matrix  with  q,  u,  and 
alpha  as  outputs.  After  obtaining  the  A,  B,  C,  and  D  matrices 
it  is  next  necessary  to  determine  the  transmission  zeros  and 
to  be  sure  that  there  are  no  decoupling  zeros  by  using  the 
computer  program  ZERO  (Ref  9).  Alternately,  decoupling  zeros 
may  be  determined  by  examining  the  system  matrix  P  (Ref  10). 
The  system  is  controllable  and  observable  if  P  has  full  rank, 
(see  equation  4-3).  A  third  method  of  checking  controllabili¬ 
ty  is  to  examine  the  M  and  M  matrices. 

c  o 

n-  1 

M-  [_B  A  B  ...  A  _B  ]  (4-1) 

c 

n-1  T 

M-  [CCA  . . .  GA  ] 

o  (4-2) 

must  have  full  rank  for  the  system  to  be  controllable, 
c 

fl  must  have  full  rank  for  the  system  to  be  observable, 
o 

After  entering  the  data  in  state-space  form,  the  program 
gives  the  decoupling  and  transmission  zeros.  There  must  be  no 
inp u t - d e co up  1 ing  zeros,  otherwise  the  system  is  uncontrolla¬ 
ble.  If  the  system  is  uncontrollable  because  B  does  not  have 
full  rank,  then  redundant  controls  are  present.  The  B  matrix 


iitl 


caa  be  made  to  have  full  rank  by  combining  the  linearly  depen¬ 


dent  columns  so  that  all  columns  are  linearly  independent. 


Ou tpu t- d ecoup 1 ing  zeros,  which  determine  observability,  possi¬ 


bly  can  be  removed  with  an  Irregular  design.  If  the  system  is 


unobservable  because  C  is  rank  deficient,  then  additional 


measurements  must  be  added  to  the  system.  This  is  done  by 


adding  an  M,  or  measurement  matrix  in  an  Irregular  design. 


Transmission  zeros  must  be  in  the  closed  left-half  plane.  If 


not,  they  possibly  can  be  moved  from  the  right-half  plane  to 


the  left-half  plane  by  using  an  Irregular  design. 


If  the  product  of  the  C  and  B  matrices  has  full  rank,  then 


a  Regular  design  may  be  attempted  if  no  transmission  zeros  are 


in  the  right-half  plane.  Although  a  design  may  be  possible, 


good  results  are  not  guaranteed.  If  it  appears  that  only 


designs  excessively  requiring  large  inputs  (e.  g.,  too  large  a 


canard  deflection)  resulting  from  large  controller  gains  are 


possible,  then  the  singular  decomposition  of  P  (Ref  11)  must 


be  examined.  The  singular  decomposition  is  a  measure  of  the 


singularity  of  a  matrix.  If  P  is  barely  nonsingular,  then 


large  gains  will  be  required  to  stabilize  the  aircraft.  This 


can  only  be  remedied  by  changing  the  system.  The  easiest  way 


to  do  this  is  to  change  the  C  and  D  matrices  by  choosing 


different  outputs.  Changing  the  B  matrix  requires  a  change  in 


control  inputs,  while  changing  the  A  matrix  requires  a 


modified  aircraft. 


P  -  I  -A  B 


■C  D 


(4-3) 


■■  .'-.V 


The  closed-loop  roots  can  be  determined  by  calculating  the 


closed  loop  transfer  matrix  G(s): 

-1 

G(s)»  C[  sl_  -  A]  B  (4-4) 

Calculation  of  the  final  transmission  zeros  (those  ob¬ 
tained  after  the  system  model,  design  parameters,  and  sampling 
period  are  chosen)  provides  a  theoretical  limit  to  the  perfor¬ 
mance  of  the  system  with  the  chosen  design.  As  the  sampling 
period  approaches  zero  in  the  limit,  the  sampling  frequency 
approaches  infinity.  As  the  sampling  frequency  approaches 
infinity,  some  of  the  closed-loop  eigenvalues  asymptotically 
approach  the  transmission  zeros. 

In  order  to  calculate  the  theoretical  limit  the  B  matrix 


is  transformed  into 
has  as 


£ 

B 


form.  The  £  matrix  is  square  and 
2 


many  rows  as  system  inputs.  If  enough  zero  rows  already 
exist,  then  the  order  of  the  state  equations  is  changed  to 
place  the  zero  rows  on  top.  If  not,  then  a  transformation 
matrix  T  must  be  found  to  change  B  into  the  required  form.  As 
the  sampling  frequency  approaches  infinity,  the  closed-loop 
roots  are  driven  to  the  final  transmission  zeros.  Note  that 
this  analysis  only  applies  to  the  limits  of  performance,  not 
to  the  actual  performance  itself. 

-  TB  (4-5) 


Transforming  the  B  matrix  has  three  advantages.  First, 


the  transmission  zeros  can  be  calculated  by  hand,  thus 


validating  the  results  obtained  from  the  computer  program 
ZERO.  Second,  in  the  limit  as  the  sampling  frequency  goes  to 
infinity,  decoupling  of  outputs  is  guaranteed.  Lastly,  for 
Irregular  designs,  algorithms  exist  for  picking  M  (measurement 
matrix)  elements.  Unfortunately,  these  algorithms  may  not 
select  M  elements  that  produce  good  results,  and  the  designer 
may  have  to  resort  to  guessing  which  elements  to  pick. 

Matrix  manipulation  cau  be  facilitated  with  the  use  of 
MATLAB,  a  computer  program  developed  by  Cleve  Moler  at  the 
University  of  New  Mexico. 

Next,  the  design  parameters  specifying  the  digital 
controller  are  selected. 

Finally,  steady-state  analysis  of  performance  can  be 
accomplished  by  using  G(s)  and  R(s)  (formed  from  the  commanded 
inputs).  Transient  analysis  is  evaluated  using  control  rate 
limits  . 

To  summarize,  preparation  for  selection  of  design 
parameters  (steps  1-5),  design  (step  6),  and  evaluation  (steps 
7-10)  require  the  following: 

1.  Validation  of  the  aircraft  model  (A,B,C,  and  D  matrices). 

2.  Selection  of  outputs  (C  matrix). 

3.  Determination  of  Regular  or  Irregular  design. 

4.  Check  of  the  controllability  (B  and  A  matrices)  and 

observability  (C  and  A  matrices)  of  the  system.  This  may  be 

accomplished  by  using  a  computer  program,  examining  P,  or 

examining  M  and  M  . 

c  o 

5.  Determination  and  evaluation  of  the  transmission  zeros 
using  a  computer  program ,  Note  that  Porter's  method  works 


equally  well  for  both  analog  and  digital  systems.  However, 
the  discrete  time  transmission  zeros  are  also  affected  by  the 


sampling  period.  As  the  sampling  period  approaches  zero, 
increasingly  tight  tracking  occurs. 

6.  Calculation  of  G(s).  (A,B,  and  C  matrices) 

7.  Transformation  of  the  B  matrix  using  a  T  matrix  to  obtain 
transmission  zeros,  guaranteed  decoupling  of  outputs,  and 
algorithms  for  picking  M  elements  (M  matrix  only  used  for 
Irregular  design).  (If  possible,  not  required). 

8.  Selection  of  design  parameters  (see  Application  of 
Porter's  Method)  .  Evaluation  of  singularity  using  singular 
decomposition  if  required. 

9.  Evaluation  of  steady-state  performance  using  G(s)  and 
commanded  inputs  R(s). 

10.  Evaluation  of  transient  performance  based  upon  control 
rate  limitations. 

App lication  of  Porter's  Method 

The  state  and  output  equations  for  flight  condition  number 
three  (1.2  Mach  at  IS  000  feet)  are: 
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Using  the  computer  program  ZERO  to  obtain  the  decoupling 
and  transmission  zeros  yields: 

1.  There  are  no  decoupling  zeros. 

2.  A  transmission  zero  exists  at  (0,0).  This  is  expected 
since  the  number  of  transmission  zeros  should  be  equal  to  the 
number  of  states  minus  the  number  of  inputs  (Ref  12). 

3.  The  above  values  are  obtained  for  the  other  four 
flight  conditions  as  well. 

Note  that  the  state-space  equations  have  !~0  ~ I  form. 

L5  J 

The  resulting  product  below  is  recognized  as  having  full  rank. 

CB-C  B  (4-8) 

2  2 

Thus  a  Regular  design  is  permitted.  Since  and  A^  have  full 

12  21 

rank,  it  is  seen  that  the  £  matrix  in  equation  (4-1)  is  non¬ 
singular.  Therefore,  the  system  is  controllable  and 
observable,  verifying  the  results  obtained  from  the  program 
ZERO. 


The  plant  matrix  can  be  partitioned  into  the  form: 


A  - 


where 


A 

11 

A 

21 

A  is 
11 


A 

12 

A 

22 


(4 


(n-m)  x  (n-m)  with  n«number  of  states  and 


manunber  of  inputs. 


Note  that  no  transmission  zeros  are  in  the  right-half 
plane.  Therefore,  a  regular  design  is  permitted. 

Since  the  B  matrix  could  not  be  easily  transformed,  a  T 
matrix  is  not  explicitly  determined.  MULTI  uses  any  form  of 
the  B  matrix  when  the  controller  matrices  are  computed. 

Next  it  is  necessary  to  choose  the  best  set  of  design 
parameters  to  obtain  a  controller  that  gives  the  best 
performance.  For  both  Regular  and  Irregular  designs  the 
following  parameters  must  be  chosen: 

1.  ALPHA-  Proportion  of  integral  to  direct  feedback  such 

that  K*ALPHA  *  jC.  A  special  incentive  exists  to  try  to 
0  1 

achieve  a  design  with  ALPHA  equal  to  one  since  this  reduces 

the  required  gains  by  half.  Using  "cache-memory",  controller 

gains  can  be  stored  in  registers  (Ref  1 H)  .  Since  K  and  K  are 

0  1 

equal  when  ALPHA  is  one,  the  number  of  required  registers  is 
reduced  by  a  factor  of  two.  Heuris tical ly ,  ALPHA  should  be 
greater  than  zero  and  less  than  two. 

2.  EPSILON-  Sigma  matrix  multiplier.  Heuris tica 1 ly , 
EPSILON  should  be  greater  than  zero  and  less  than  four. 

3.  SIGMA  weighting  matrix-  SIGMA  has  p  weighting 
diagonal  entries.  The  individual  SIGMA  elements  correspond  to 
the  respective  commanded  inputs  and  outputs.  Because  the 
SIGMA  matrix  is  diagonal,  the  outputs  can  be  decoupled. 

4.  SAMP l. E  PERIOD-  This  thesis  uses  0.025  seconds.  The 
SAMPLE  PERIOD  should  be  greater  than  0.01  seconds  and  less 
than  0.10  seconds.  As  the  SAMPLE  PERIOD  is  reduced,  the 


closed-loop  roots  of  the  model  migrate  to  the  transmission 
zeros.  If  transmission  zeros  should  exist  in  the  right-half 
plane,  then  there  is  a  lower  limit  on  the  allowable  value  of 
the  SAMPLE  PERIOD. 

5.  Commanded  inputs  are  initially  set  to  some  maneuver 

values  that  the  designer  is  sure  the  aircraft  can  perform. 

The  magnitude  of  the  steady-state  control  inputs  u  to  the 

s  s 

aircraft  for  a  desired  maneuver  y  can  be  computed  from  the 

s  s 

open-loop  aircraft  transfer  function  G(s)  by  using: 

-1 

u  -[£(0)]  y  ( A- 10 ) 

s  s  s  s 

To  start,  the  design  parameters  ALPHA,  EPSILON,  and  SIGMA 
are  set  equal  to  one.  Then  the  K  matrices,  or  controllers, 
are  examined.  In  general,  the  entries  of  the  K  matrices 
should  have  orders  of  magnitude  greater  than  or  equal  to 
negative  two  and  less  than  or  equal  to  two.  In  addition  to 
most  designs  being  in  this  general  area,  the  dynamic  range  of 
the  elements  is  minimized.  This  aids  in  producing  computa¬ 
tionally  correct  results. 

The  following  is  a  suggested  design  procedure  for  using 
the  computer-aided-design  tool  MULTI.  It  should  be  noted  that 
no  unique  solution  exists  to  the  design  process.  Tradeoffs  are 
made  between  transient  and  steady-state  response,  between 
maximum  commanded  inputs  and  control  input  responses,  and 
between  length  of  maneuver  time  and  control  input  rates. 
Additional  tradeoffs  must  be  made  in  the  search  for  a  robust 


controller.  The  design  parameters  affect  one  another,  with  the 
greatest  coupling  between  ALPHA  and  EPSILON,  and  among  the 
SIGMA  elements . 

1.  Aircraft  model: 

Enter  the  A,  B,  C,  and  D  matrices  in  state-space  form. 
For  an  initial  design  attempt,  no  sensors  or  actuators  are 
used . 

2.  Simulation  parameters: 

Set  the  state  and  integrator  initial  values,  X(0)  and 
Z(0).  Since  the  aircraft  data  provided  are  for  steady  level 
trimmed  flight,  these  vectors  are  set  to  zero.  Since  only 
perturbations  are  of  interest,  the  trim  values  of  alpha  and 
theta  do  not  affect  the  simulation.  However,  it  should  be 
noted  that  if  these  trim  values  are  large,  the  results  could 
show  that  additional  increments  are  possible  when  in  fact  they 
would  exceed  the  allowable  values. 

Set  the  commanded  inputs  that  correspond  to  the  maneuver 
to  be  performed.  Choose  a  sampling  period  (e.  g.,  .025  sec¬ 
onds).  Set  the  simulation  time  (for  this  thesis  16  seconds) 
to  be  at  least  as  long  as  the  commanded  input.  Next,  set  the 
calculation  step  size  in  time  equal  to  the  sampling  period. 

For  an  initial  design  attempt  no  lags  are  used.  The  lag 
capability  in  MULTI  models  the  synchronous  ana  log- to-d igi ta 1 
delay.  The  computation  delay,  which  is  asynchronous,  is  usu¬ 
ally  much  shorter  than  the  ana  log- to-d i gi ta 1  delay.  The  digi¬ 
ts  1- to-ana  log  delay  is  negligible.  No  capability  currently 
exists  in  MULTI  for  modeling  input  "dead-zones",  the  delay 
between  the  commanded  start  of  the  control  input  and  when  the 


control  inputs  actually  start  to  respond.  For  example,  with¬ 
out  afterburners,  the  dead-zone  for  thrust  is  approximately 
0 . 5  seconds . 

Lastly,  set  the  control  input  limits: 

Canards  (-60,30)  degrees 

Strakes  (-19.5,13)  degrees 

Thrust  (-14782,2000)  pounds 

Note  that  these  limits  are  set  from  trim  values.  Since  14  782 
pounds  are  required  to  achieve  steady  level  flight,  total 
thrust  may  be  decreased  by  the  same  amount. 

3.  Design  parameters: 

Initially,  ALPHA,  EPSILON,  and  SIGMA  are  set  to  one.  If 
all  of  the  entries  of  the  controller  matrices  exceed  the 
specified  orders  of  magnitude,  then  EPSILON,  the  multiplier  of 
these  entries,  is  adjusted  in  the  opposite  direction.  Next, 
the  entries  of  each  column  in  the  controller  matrix  are  exam¬ 
ined.  Again,  if  the  entries  in  a  particular  column  exceed  the 
heuristically  determined  controller  limits,  the  corresponding 
SIGMA  element  is  adjusted.  Once  the  controller  entries  are 
adjusted  to  a  best  initial  guess,  a  simulation  is  run.  MULTI 
alerts  the  designer  if  control  inputs  are  exceeded. 

If  control  inputs  are  exceeded,  the  designer  should  make 
sure  that  the  maneuver  commanded  is  one  the  aircraft  is 
capable  of  performing.  If  so,  the  design  parameters  should  be 
adjusted  by  small  increments,  and  the  simulation  rerun.  Once 
the  simulation  produces  inputs  that  are  within  limits,  the 
outputs  should  be  examined. 


The  outputs  are  quickly  examined  using  the  option  that 
provides  figures  of  merit.  Design  parameters  are  retuned 
until  the  best  figures  of  merit  are  obtained.  Next,  plots  of 
outputs,  states,  and  inputs  are  examined  at  the  terminal.  If 
satisfactory,  the  commanded  inputs  are  increased  until  control 
input  limits  are  reached  in  order  to  determine  the  largest 
allowable  command.  Lastly,  tradeoffs  are  made  between  design 
objectives.  Note  that  overshoots  can  sometimes  be  reduced  by 
increasing  ALPHA.  When  a  satisfactory  design  is  achieved, 
calcomp  plots  are  generated. 

Results 

The  design  parameters  and  maximum  commanded  inputs  selec* 
ted  are  presented  in  Table  4-1.  The  controller  matrices 
resulting  from  the  selected  design  parameters  are  presented  in 
Table  4-2.  Figures  of  merit  are  provided  in  Table  4-3. 
Control  rates  are  provided  in  Table  4-4.  Calcomp  plots  of 
outputs,  commanded  inputs,  control  inputs,  theta,  and  flight- 
path  angle  are  provided  in  Figures  4-1  to  4-19.  Note  that 
only  nonzero  commanded  inputs  are  plotted.  Similar  information 
is  provided  for  the  other  four  flight  conditions  in  appendices 
B,  D,  E,  F,  G,  and  H. 

The  first  longitudinal  maneuver,  direct  climb,  is  best 
performed  by  the  flight  condition  at  1.2  Mach,  15  000  feet.  A 
smoothed  pulse  pitch  rate  of  approximately  47  degrees/second 
is  commanded  and  closely  followed  (see  Figure  4-1).  The 


TABLE  4-1 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS  FOR  LONGITUDINAL  CONTROLLER 
Flight  Condition  -  1.2  Mach,  15  000  feet 


IT 

Maneuver 

IT 

Command  Vector  v 

5T 

ALPHA 

51  EPSILON 

5T 

SIGMA 

5T 

IT 

IT 

* 

5T 

5T 

5T 

(Diagonal  values) 

5T 

IT 

IT 

IT 

5T 

51 

51 

IT 

Direct 

NT 

1.2,. 8, 2, 3. 2 

IT 

5T 

5T 

5T 

IT 

Climb 

IT 

1.2,405,16,16 

5T 

1.1 

5T  0.9 

5T 

2,1,0 

5T 

IT 

IT 

0 , 0 , 0 , 0 

IT 

5T 

5T 

5T 

IT 

IT 

5T 

5T 

5T 

5T 

IT 

IT 

5T 

5T 

5T 

5T 

IT 

Direct 

IT 

.2, .45, .5, .7 

5T 

5T 

5T 

5T 

IT 

Lif  tl 

IT 

0 ,0,0,0 

IT 

5T 

5T 

5T 

IT 

IT 

0,0, 0,0 

5T 

5T 

5T 

5T 

IT 

IT 

5T 

5T 

5T 

5T 

5T 

IT 

5T 

5T 

5T 

5T 

IT 
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IT 

1.6,  .72,2.4,3.2 

5T 

5T 

5T 

5T 

IT 

Lif  t2 

IT 

0,0, 0,0 

IT 

5T 

5T 
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IT 
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IT 
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5T 
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Pitch- 

IT 

.8,  .0225,1.2,1.6 
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5T 

NT 

Pointing 

5T 

0 ,0 ,0 ,0 

5T 

1.1 

5T  0.9 

5T 

2,1,  .15 

5T 

5T 
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5T 

5T 
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IT 

5T 

5T 
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5T 
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IT 
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IT 

0,0, 0,0 

5T 

5T 

5T 

5T 

IT 

Trans  la  - 

IT 

0 ,0 ,0 ,0 

5T 

5T 

5T 

5T 

IT 
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IT 

1.8, .027,3.2,4.5 

IT 

5T 

5T 

5T 

5T 

IT 

5T 

.82 

5T  .  9 

5T 

1.5,1.  , .067 

5T 

IT 

IT 

-5T 

5T 

5T 

5T 

IT 
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IT 
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IT 

5T 

5T 

5T 

IT 
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IT 
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5T 

5T 

5T 
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5T 

5T 

5T 

5T 

IT 

5T 

IT 

5T 

5T 

5T 

IT 

IT 

IT 

5T 

5T 

5T 

♦Note:  Each  pulse  entry  in  has  four  parts: 

1)  The  time  (seconds)  the  input  reached  s teady- s ta te , 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


TABLE  4-2 


CONTROLLER  MATRICES 


Flight  Condition  -  1.2  Mach,  15  000  feet 


IT  Maneuver  IT  Command  Vector  v  IT  ALPHA  II 


irect  IT  1.2, .8, 2, 3. 2 

limb  IT  1.2,405,16,16 

IT  0 , 0 , 0 , 0 


Direct 


IT  Li  f  1 2 


IT  Pitch-  IT  .8  , .0225,1.2  , 
TT  Pointing  IT  0,0, 0,0 
IT  IT  .8,  .0225,16,16 


0 ,0 ,0 , 0 


1.6,  .72,2.4,3.2  IT 
0 ,0 ,0 , 0  IT 
0 , 0 , 0 , 0  IT 


.8  , .0225, 1.2, 1.6  IT 
0,0, 0,0  IT  1.1 

n  a  r  *  r  a  r  «v 


Vertical  IT  0,0, 0,0  IT 

Transla-  IT  0,0, 0,0  IT 

tionl  IT  1.8  , .027,3.2,4.5  IT 


Vertical  IT 
Transla-  IT 
tion2  IT 


,0,0,0 

,0,0,0 

.8, .024,16,16 


. 1902E+01 
- . 1006E+02 
. 31 64E- 02 


. 1630E-02 
-  .5156E-01 
. 9144E-02 


IT  .  1902E+1  .  1 6 3 0E -  2 

1T-.1006E2  -  .  5156E-1 

IT  .3164E-2  .9144E-2 


H 


IT  .  1064E1  .1215E-2 

1T-.  5624E1  -  .  3844E-  1 

1T.  1769E-2  .  6816E-  2 


u 


■  .  2  5  7  2E2  IT 
•  .  8040E2  IT 
-.2141  IT 


.  8563E1  IT 
.  26  7  7E2  IT 
.  7 1 3E-  1  IT 


♦Note:  Each  pulse  entry  in  v  has  four  parts: 

1)  The  time  (seconds)  the  input  reached  steady-state 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero , 

4)  and  the  time  the  input  returns  to  zero. 


**No  te :  K 


ALPHA  *  K 


TABLE  4-3 


LONGITUDINAL  FIGURES  OF  MERIT 
Flight  Condition  -  1.2  Mach,  15  000  feet 


51  Maneuver  IT  Command  Vector  v 


Direct 

IT 

1.2,  .8, 2, 3. 2 

IT 

15 

IT 

1.65 

IT 

.8113 

Climb 

IT 

1.2,405,16,16 

IT 

1.35 

IT 

1 . 5 

5T 

408.6827 

IT 

0,0, 0,0 

IT 

8 . 70 

IT 

2.1 

IT 

.3000 

IT  Direct 


IT  Lif  t2 


0 , 0 , 0 , 0 


IT  12.15 


5T  IT 

1.6, .72,2.4 
0 , 0 , 0 , 0 

0 , 0 , 0 , 0 

,3.2 

5T  15  IT  1.95 

IT  15  IT  2.55 

IT  15  IT  2.4 

IT  IT 

Pitch- 

IT 

.8  ,  .0225,1.2,1.6 

IT 

15 

IT 

1.05 

Pointing 

IT 

0 , 0 , 0 , 0 

IT 

15 

IT 

1 . 2 

IT 

.8  ,  .0225,16,16 

IT 

4.8 

IT 

15 

IT  4.00E-8 


.  7501 
-.5716 
.  2888 


255 

0196 


02250 


Vertical  IT 
Transla-  IT 


,0,0,0  IT  15  IT  1.95 

,0,0,0  IT  15  IT  1.95 

.8, .027,3. 2,4.5  5T  15  IT  3.3 


033 

0089 

230 


Vertical  IT 
Transla-  IT 
tion2  IT 


,0,0,0 

,0,0,0 

.8  , .024,16  ,16 


5T  15  5T  1.05 
5T  15  IT  1.2 
IT  6.15  5T  15 


037 

0094 

240 


♦Note:  settling  time  in  seconds 

♦♦Note:  peak  time  in  seconds 

***Note:  peak  value  (greatest  magnitude),  angles  in  radians 

Note:  commanded  inputs  are  smoothed 

Note:  commanded  inputs  in  order  are  -  pitch  rate,  change 

in  velocity,  and  ang 1 e - of  -  a t tack 
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forward  velocity  is  increased  by  405  feet/second  and  held. 
The  commanded  u  input  and  output  are  indistinguishable  (see 

Figure  4-2).  The  largest  maximum  commanded  inputs  (in  de¬ 
creasing  order)  are  at  1.2  Mach,  15  000  feet;  0.9  Mach,  sea 

level;  0.4  Mach,  sea  level;  0.7  Mach,  15  000  feet;  and  0.9 
Mach,  50  000  feet. 

Since  thrust  is  an  input,  velocity  control  is  excellent 
for  all  maneuvers  performed.  The  additional  thrust  required 
to  increase  velocity  ranges  from  320  pounds  at  0.4  Mach,  sea 
level,  to  1000  pounds  at  0.9  Mach,  sea  level.  At  1.2  Mach, 
15  000  feet,  340  additional  pounds  of  thrust  is  required  (see 
Figure  4-3).  The  limiting  input  for  direct  climb  is  the 

strakes  (see  Figure  4-3).  The  direct  climb  maneuver  commands 
alpha  to  zero,  so  f-’ght-path  angle  should  equal  theta  in  the 
steady-state  since: 

gamma  ■  theta  -  alpha  (4-11) 

An  examination  of  the  gamma  and  theta  plots  reveals  that 
after  an  initial  transient,  gamma  does  indeed  track  theta 
closely  (see  Figure  4-4).  The  best  steady-state  tracking 
(performance  is  presented  in  decreasing  order)  is  performed  by 
1.2  Mach,  15  000  feet;  0.9  Mach,  sea  level;  0.7  Mach,  15  000 
feet;  0.9  Mach,  50  000  feet;  and  0.4  Mach,  sea  level.  The 
best  transient  response  is  provided  by  0.9  Mach,  sea  level; 
1.2  Mach,  15  000  feet;  0.7  Mach,  15  000  feet;  0.4  Mach,  sea 
level;  and  0.9  Mach,  50  000  feet.  Steady-state  gamma  varies 
from  90  degrees  (more  could  have  been  commanded,  but  the 
assumption  of  small  perturbations  is  stretched  here  at  best) 


at  1.2  Mach,  15  000  feet  to  31  degrees  at  0.9  Mach,  50  000 

f  eet . 


Maximum  allowable  deflection  rates  for  the  canards  and 
strakes  are  100  degrees/second  and  30  degrees/second, 
respectively.  Actual  canard  deflection  rates  are  well  within 
the  allowable  maximum  for  all  flight  conditions.  Strake 
deflection  rates,  however,  greatly  exceed  the  maximum  at  all 
flight  conditions.  For  example,  the  maximum  negative  and 
positive  strake  deflection  rates  at  1.2  Mach,  15  000  feet,  are 
-173.6  and  131.9  degrees/second.  Whenever  control  input  rates 
are  exceeded,  the  maneuver  time  must  be  lengthened  and  the 
commanded  steady-state  value  must  be  decreased,  but  the  con¬ 
troller  does  not  have  to  be  retuned. 

Direct  liftl,  a  commanded  pitch-rate  pulse  of  only  0.7 
seconds  duration,  has  the  same  controller  at  each  flight 
condition  as  direct  climb.  For  this  maneuver,  the  largest 
commanded  inputs  are  at  1.2  Mach,  15  000  feet;  0.7  Mach,  15 
000  feet;  0.9  Mach,  sea  level;  0.6  Mach,  sea  level;  and  0.9 
Mach,  50  000  feet.  Transient  alpha  responses  are  best  mini¬ 
mized  at  0.9  Mach,  sea  level;  1.2  Mach,  15  000  feet;  0.7  Mach, 
15  000  feet;  0.4  Mach,  sea  level ;  and  0.9  Mach,  50  000  feet. 
At  1.2  Mach,  15  000  feet,  pitch-rate  is  commanded  to  a  27 
degree/second  pulse,  while  alpha  is  commanded  to  zero.  Pitch- 
rate  tracks  the  commanded  output  excellently,  while  after  a 
short  transient  alpha  goes  to  zero  (se”  Figure  4-5).  The 
change  in  forward  velocity  is  commanded  to  zero.  The  u 


transient  is  small,  with  the  greatest  magnitude  a  decrease  in 


u  of  -0.21  feet/second  occurring  at  1.2  Mach,  15  000  feet  (see 
Figure  4-6) . 

Again,  the  strake  is  the  limiting  control  input  to  a 
greater  commanded  input.  A  negligible  increase  in  thrust  of 
16  pounds  is  required  at  1.2  Mach,  15  000  feet  (see  Figure  4- 
7).  Again,  since  alpha  is  commanded  to  zero,  flight-path 
should  track  theta.  In  the  s teady- s ta te  ,  gamma,  the  flight- 
path  angle  does  indeed  equal  pitch  angle  at  all  flight  condi¬ 
tions.  The  best  transient  tracking  of  theta  by  gamma  is 
provided  at  0.9  Mach,  sea  level;  1.2  Mach,  15  000  feet;  0.7 
Mach,  15  000  feet;  0.4  Mach,  sea  level;  and  0.9  Mach,  50  000 
feet.  A  steady-state  flight-path  angle  of  13  degrees  is 
achieved  at  1.2  Mach,  15  000  feet  (see  Figure  4-7). 

Unfortunately,  allowable  control  input  rates  are  again 
exceeded.  For  example,  maximum  strake  rates  are  -166  and  167 
degrees/second  at  1.2  Mach,  15  000  fe-  t*  Maximum  canard  rates 
are  within  acceptable  limits  except  at  0.4  Mach,  sea  level 
where  the  maximums  are  -141  and  143  degrees/second. 

Direct  lift2,  which  is  the  same  maneuver  performed  with 
longer  maneuver  time,  permits  greater  commanded  steady-state 
values.  In  decreasing  order  of  capability,  the  greatest 
commanded  inputs  are  at  1.2  Mach,  15  000  feet;  0.4  Mach,  sea 
level;  0.9  Mach,  sea  level;  0.7  Mach,  15  000  feet;  and  0.9 
Mach,  50  000  feet.  At  1.2  Mach,  15  000  feet,  pitch-rate  is 
commanded  to  a  pulse  of  41  deg ree s / s econd  (see  Figure  4-8). 
Again,  excellent  tracking  is  performed  of  the  commanded  input 
(see  Figure  4-8).  A  minimal  undershoot  of  pitch-rate  output 
occurs  at  the  end  of  the  commanded  input,  but  it  decays  within 


Again,  the  strakes  are  the  limiting  control  input.  The 
shape  of  the  control  inputs  is  similar,  with  a  steady-state 
negative  strake  deflection  required  to  provide  a  nose-down 
pitching  moment  (see  Figure  4-9). 

Again,  in  the  steady-state  gamma  equals  theta.  The  best 
transient  tracking  of  theta  by  gamma  occurs  at  0.9  Mach,  sea 
level;  1.2  Mach,  15  000  feet;  0.7  Mach,  15  000  feet;  0.4  Mach, 
sea  level;  and  0.9  Mach,  50  000  feet.  At  1.2  Mach,  15  000 
feet  a  steady-state  flight-path  angle  of  83  degrees  is 
achieved  (see  Figure  4-10). 

Note  that  an  increase  in  maneuver  time,  although  coupled 
with  an  increase  in  the  commanded  steady-state  value,  causes  a 
marked  reduction  in  control  input  rates.  The  previously  ex¬ 
ceeded  maximum  canard  rates  decline  to  -35  and  59  degrees/sec¬ 
ond  at  the  worst-performing  flight  condition,  0.4  Mach,  sea 
level.  At  1.2  Mach,  15  000  feet,  maximum  canard  rates  are  -8 
and  16  degrees/second .  The  maximum  strake  rates  are  -86  and 
41  degrees/second. 

Pitch-pointing,  the  maneuver  that  commands  pitch-rate  and 
ang le - o f - a t tack  such  that  flight-path  is  zero,  has  a  separate 
controller.  The  greatest  commanded  inputs  are  at  0.4  Mach, 
sea  level;  0.7  Mach,  15  000  feet;  1.2  Mach,  15  000  feet;  0.9 
Mach,  sea  level;  and  0.9  Mach,  50  000  feet.  The  best 
transient  and  steady-state  tracking  of  the  commanded  inputs  by 
the  outputs  occurs  at  0.4  Mach,  sea  level;  0.9  Mach,  50  000 
feet;  0.7  Mach,  15  000  feet;  1.2  Mach,  15  000  feet;  and  0.9 


Mach,  sea  level. 


At  0.4  Mach,  sea  level,  2.1  degrees  of 


J  OUTPUT  TO  OIREC  LIFT2 


CRNflRD  i  S TRRKE  INPUTS  TO  DIRECT  LIFT2:I.2H.  ' 5K  FT 


Figure  4-9 
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pi tch- po in t ing  (the  value  of  alpha  and  theta)  is  possible.  At 
1.2  Mach,  15  000  feet,  the  response  is  degraded  to  1.2  degrees 
of  p i tch- po in t i ng  (see  Figure  4-11). 

The  input  responses  at  the  five  flight  conditions  are 
similar.  Again,  the  strake  is  the  limiting  control  input. 
Negative  steady-state  control  values  occur  in  which  the  canard 
deflection  i3  approximately  half  that  of  the  strake 
deflection  (see  Figure  4-12).  The  deflections  provide  cancel¬ 
ing  pitching  moments  about  the  center  of  gravity,  resulting  in 
a  net  rotation  of  the  nose  about  the  center  of  gravity.  Since 
the  symmetrical  flaperons,  and  not  the  strakes,  are  intended 
to  be  the  primary  lift-vector  surfaces,  the  strakes  are  a  poor 
match  against  the  canards. 

The  flight  conditions  minimizing  gamma  steady-state  error 
are  0.4  Mach,  sea  level;  0.9  Mach,  50  000  feet;  0.7  Mach,  15 
000  feet;  1.2  Mach,  15  000  feet;  and  0.9  Mach,  sea  level.  At 
1.2  Mach,  15  000  feet,  1.6  degrees  of  p i tch - po in t ing  with  a 
gamma  steady-state  error  of  .35  degrees  is  achieved. 

Maximum  canard  rates  are  well  within  allowable  limits,  but 
the  maximum  strake  rates  are  again  violated.  At  the  best¬ 
performing  flight  condition,  0.4  Mach,  sea  level,  the  maximum 
strake  rates  are  -39  and  27  d eg r e e s / s e c ond .  At  1.2  Mach,  15 
000  feet,  the  maximum  canard  rates  are  -18  and  6 
degrees / s econd .  The  maximum  strake  rates  at  this  flight 
condition  are  -58  and  22  d eg r e e s / s e c ond  . 

One  controller  is  used  for  both  vertical  translationl  and 
vertical  t r ans 1  a t i on 2  .  The  largest  commanded  alpha  inputs 
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Figure  4-11 
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occur  at  0.7  Mach,  15  000  feet;  0.4  Mach,  sea  level;  1.2  Mach, 
15  000  feet;  0.9  Mach,  sea  level;  and  0.9  Mach,  50  000  feet. 
Pitch-rate  is  commanded  to  zero,  but  the  resulting  positive 
transient  is  canceled  by  an  undershoot.  In  general,  the  alpha 
output  follows  the  general  shape  of  the  commanded  input, 
although  the  tracking  lags  by  about  1.6  seconds.  At  1.2  Mach, 
15  000,  alpha  is  commanded  to  1.3  degrees  (see  Figure  4-14). 

The  control  input  shapes  are  similar  at  all  flight 
conditions,  with  the  strake  hitting  the  control  limit.  Thrust 
control  provides  excellent  u  response  (see  Figure  4-15). 

Pitch  angle,  which  should  be  zero  since  pitch-rate  is 
commanded  to  zero,  does  have  a  positive  transient  that  even¬ 
tually  decays  to  zero.  At  1.2  Mach,  15  000  feet,  the  pitch 
angle  transient  is  .45  degrees,  while  the  resulting  flight- 
path  angle  is  -0.9  degrees  (see  Figure  4-16).  The  best  tran¬ 
sient  response  is  provided  by  0.4  Mach,  sea  level;  0.9  Mach, 
50  000  feet;  0.7  Mach,  15  000  feet;  1.2  Mach,  15  000  feet;  and 
0.9  Mach,  sea  level. 

Maximum  canard  control  input  rates  are  within  allowable 
limits  for  all  flight  conditions.  Maximum  strake  control 
input  rates  are  within  allowable  limits  for  all  flight 
conditions  but  0.9  Mach,  50  000  feet;  -40  and  41 
degrees/second.  At  1.2  Mach,  15  000  feet,  maximum  canard 
rates  are  -5  and  7  degrees/second.  Maximum  strake  rates  are 
-14  and  18  degrees/second. 

In  vertical  trans lation2 ,  alpha  is  commanded  to  a  fixed 

The  greatest  commanded  alpha  inputs  are 
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steady-state  value. 
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Figure  4-14 


at  0.7  Mach,  15  000  feet;  0.4  Mach,  sea  level;  1.2  Mach,  15 
000  feet;  0.9  Mach,  sea  level;  and  0.9  Mach,  50  000  feet.  The 
best  alpha  transient  response  occurs  at  0.9  Mach,  50  000  feet; 
0.4  Mach,  sea  level;  0.7  Mach,  15  000  feet;  1.2  Mach,  15  000 
feet;  and  0.9  Mach,  sea  level.  A  small  transient  pitch-rate 
occurs,  but  this  quickly  decays  (see  Figure  4-17). 

The  shape  of  the  control  inputs  varies  somewhat  from 
flight  condition  to  flight  condition.  At  two  of  the  flight 
conditions  the  strake  transient  response  reaches  the  control 
input  limit.  At  the  other  flight  conditions,  however,  the 
steady-state  strake  response  is  the  limiting  factor.  At  1.2 
Mach,  15  000  feet,  the  steady-state  canard  and  strake  deflec- 
tions  are  -9  and  -30  degrees,  respectively  (see  Figure  4-18). 

Examining  the  pitch  and  flight-path  angle  plots  reveals 
that  the  minimum  theta  steady-state  errors  occur  at  0.4  Mach, 


sea 

level;  0.9  Mach,  50  000 

f  eet ; 

0.7 

Mach  , 

15  000 

feet;  1.2 

Mach 

,  15  000  feet;  and  0.9 

Mach  , 

sea 

level . 

At  1.2 

Mach,  15 

000 

feet,  the  steady -state 

the  ta 

and 

gamma 

are  .5 

and  -  .  8 

degrees,  respectively  (see  Figure  4-19). 
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fllPm.Q.  DECREES. DECREES/StT 


A  Regular  design,  which  requires  that  the  product  CB  have  full 
rank,  is  a  minimum  phase  application.  The  input  and  output 
matrices  have  full  rank,  permitting  a  Regular  design.  It  is 
stressed  that  constant  output  matrices  must  be  used  throughout 
the  flight  envelope.  Otherwise  the  pilot  has  to  flip  a  switch 
to  tell  the  aircraft's  flight  control  system  which  maneuver  he 
is  about  to  command.  It  is  noted  that,  with  the  advent  of 
automatic  digital  flight  control  systems,  (such  as  the  integ¬ 
rated  fire/flight  control  system)  maneuvers  are  now  possible 
that  the  pilot  would  otherwise  be  incapable  of  commanding 
because  of  his  reaction  time. 

A  description  is  provided  of  the  six  longitudinal 
maneuvers  which  the  aircraft  is  required  to  perform.  At  each 
flight  condition  three  controllers  are  designed.  One  control¬ 
ler  for  the  direct  climb,  direct  liftl,  and  direct  lift2 
maneuvers  is  used  with  minimal  degradation  of  performance. 
One  Controller  for  p i tch - p o i n t ing  is  designed.  One  controller 
for  vertical  translationl  and  vertical  translation2  is  de¬ 
signed.  These  consolidated  controllers  are  used  in  the  quest 
for  a  single  robust  controller  for  all  maneuvers  at  all  flight 
conditions.  An  acceptable  alternative  would  be  a  single  con¬ 
troller  at  each  flight  condition  with  g a i n - s ch edu 1 ing  used 
between  design  points. 

A  general  list  of  assumptions  is  presented  for  the 
development  of  aircraft  equations  of  motion.  A  ten  step 
general  outline  for  applying  Porter's  method  is  presented.  A 
discussion  of  design  parameters  is  presented.  An  outline  to 


82 


’■» -»  »f  V  «  E  »  V  T~*  V  >  ?*  TT?T  '.*'  .V1.*-1  '-*■ 


allow  the  design  engineer  to  immediately  start  using  the 
interactive  CAD  program  MULTI  is  given. 

Finally,  a  general  discussion  of  the  calcomp  plots  for  all 
forty  maneuvers  is  provided.  Specific  details  and  analysis 
are  furnished  at  a  representative  flight  condition,  1.2  Mach, 
15  000  feet.  Comparisons  between  performance  capabilities  at 
different  flight  conditions  for  each  maneuver  are  made. 

Although  allowable  control  input  rates  are  exceeded  for 
most  of  the  maneuvers,  it  is  shown  that  only  the  maximum 
commanded  inputs  must  be  reduced.  Retuning  of  the  controllers 
is  not  required. 
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CHAPTER  V 


LATERAL  CONTROLLER  DESIGNS 


Introduction 

This  chapter  presents  the  results  of  the  lateral 
controller  designs  for  the  five  flight  conditions  for  each 
of  two  maneuvers.  A  sample  flight  condition,  1.2  Mach,  15 
000  feet,  is  presented.  Although  each  flight  condition  is 
discussed,  supporting  state-space  matrices,  design 
parameters,  commanded  inputs,  controller  matrices,  figures 
of  merit,  and  plots  for  the  other  four  flight  conditions  are 
included  in  the  appendices. 

Since  the  X-29A  demonstrator  aircraft's  longitudinal  and 
lateral  equations  of  motion  are  decoupled,  separate  longitu¬ 
dinal  and  lateral  controllers  may  he  designed.  In  addition, 
the  controllers  may  be  implemented  on  separate  computer 
systems,  thus  increasing  processing  speed  and  improving 
software  reliability  and  maintenance. 

Grumman  developed  a  nonlinear  simulation  computer 
program  for  the  X-29A.  From  this  program,  a  simplified 
Tedueed-state  linearized  model  at  specified  flight 
conditions  is  used.  The  linear  approximation  is  valid,  but 
becomes  less  so  the  greater  the  perturbation  away  from  the 
specified  flight  condition.  The  lateral  states  used  are 
0,  roll  angle;^,  side  slip  angle;  p,  angular  roll-rate;  and 
r»  angular  yaw  rate. 

Differential  flapeTons  and  rudder  are  the  lateral 
control  inputs.  Note  that  control-configured  vehicle  (CCV) 


maneuvers  are  not  possible  since  onLy  two  control  inputs  are 
available.  CCV  maneuvers  use  two  control  inputs  to  cancel 
opposing  moments,  thus  providing  a  net  force  in  the  desired 
direction.  In  order  to  perform  lateral  CCV  maneuvers  for 
the  X-29A,  vertical  stabilizers  (canards)  would  have  to  be 
added  to  cancel  the  opposing  rudder  moment  about  the  center 
of  gravity.  This,  however,  would  disrupt  the  smooth  airflow 
beneath  the  aircraft  and  increase  drag. 

Porter's  method  requires  that  the  number  of  outputs  be 
equal  to  the  number  of  inputs.  Therefore  two  states,  roll 
angle  and  sideslip,  are  chosen  as  outputs  (see  Appendix  A 
for  other  design  approaches).  It  is  of  paramount  importance 
that  the  same  output  matrix  (C)  be  used.  If  not,  then  the 
pilot  has  to  notify  the  aircraft  in  advance  which  maneuver 
he  intends  to  perform.  In  addition,  different  output  ma¬ 
trices  require  additional  sensors  to  measure  the  extra  out¬ 
puts.  Redundancy  requirements  mean  that  additional  cost  and 
weight  are  added  to  the  aircraft,  and  reliability  is 
reduced . 

Description  of  Maneuvers 

As  in  longitudinal  commanded  inputs,  maximum  maneuvers 
are  commanded.  The  controller  is  designed  to  perform  two 
maneuvers,  beta-pointing  and  coordinated  turn.  The  maneuver 
times  and  commanded  steady-state  values  vary  from  flight 
condition  to  flight  condition. 

When  an  aircraft  lands  in  a  crosswind,  its  initial 
descent  aligns  the  fuselage  of  the  aircraft  parallel  with 
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the  centerline  on  the  runway.  The  wings  ere  straight  and 
level.  As  the  aircraft  starts  to  "flare",  the  pilot  turns 
the  nose  of  the  aircraft  into  the  wind  with  the  rudder, 
since  excessive  nonzero  roll  angle  is  undesirable.  The 
landing  gear  of  heavier  transport  aircraft  cannot  take  the 
strain  of  landing  canted,  nose  into  the  wind,  but  are  not  as 
easily  blown  sideways  on  the  runway.  A  more  agile  air¬ 
craft,  such  as  the  X-29A,  must  land  with  its  nose  into  the 
wind.  Beta-pointing  is  an  ideal  maneuver  for  this  situa¬ 
tion.  Sideslip  angle,  or  beta,  is  defined  as  the  angle 
between  the  velocity  vector  and  the  nose  of  the  aircraft  (in 
the  plane  formed  by  the  wings  and  fuselage).  Thus  beta¬ 
pointing  allows  the  aircraft  to  continue  traveling  down  the 
center  of  the  runway  while  the  nose  is  pointed  into  the 
wind.  Beta-pointing  is  also  useful  for  quick  minor  adjust¬ 
ments  in  situations,  such  as  flying  in  formation,  where  the 
pilot  does  not  wish  to  tufn.  Thus  beta-pointing  commands 
roll  angle  to  zero  while  commanding  beta  to  ramp  up  and  hold 
its  steady-state  value. 

In  a  coordinated  turn,  the  aircraft  is  commanded  to  a 
fixed  roll  angle,  while  sideslip  is  commanded  to  zero.  In 
performing  a  coordinated  turn,  it  is  undesirable  to  allow 
sideslip.  If  sideslip  occurs  during  a  turn,  the  nose  drops, 
and,  if  left  unchecked,  causes  the  aircraft  to  start  spiral¬ 
ing.  At  all  five  flight  conditions,  the  aircraft  is  com¬ 
manded  to  a  steady-state  roll  angle  of  88.815  degrees  while 
beta  is  commanded  to  zero. 
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Initially,  separate  controllers  are  designed  for  each 
maneuver  at  each  flight  condition.  However,  it  is  discov¬ 
ered  that  the  design  parameters  used  to  obtain  the  control¬ 
ler  for  a  coordinated  turn  at  one  flight  condition  provide 
excellent  results  in  obtaining  controllers  for  a  coordinated 
turn  at  other  flight  conditions.  Additionally,  it  is  found 
that  the  same  controller  obtained  for  a  coordinated  turn  at 
a  specified  flight  condition  can  be  used  for  beta-pointing 
as  well.  Next,  a  successful  search  is  made  for  one 
controller  that  will  produce  excellent  tracking  for  both 
maneuvers  at  all  flight  conditions.  Finally,  half  of  the 
gains  in  the  robust  controller  are  set  to  zero  with  only 
minimal  degradation  of  results. 

Since  the  aircraft  model  Grumman  used  for  the  X-29A  is 
not  available,  validation  of  the  model  is  not  performed. 
Data  set  #3  is  used  in  this  thesis  (see  Appendix  A  for  other 
design  approaches  and  data  sets  used). 

The  outputs  selected  are  two  states,  roll  angle  and 
sideslip  angle.  The  same  C  matrix  is  used  for  both  maneu¬ 
vers  at  all  flight  conditions. 

The  state  and  output  equations  for  flight  condition 
number  three  (1.2  Mach  at  15  000  feet)  are: 


0 

-  .4133 
-97.95 
18.68 


1 

.03774 

-7.959 

-.1450 


.03777 
-.9993 
2.658 
- . 4644 
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0  0 

r  1 

-.002243  .0004275 

1.169  .2437 

C 

.2339  -.08424 

L  J 

x*  Ax  +  Bu 


(5-1) 


Y-Cx-f  <j) 
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j  J  L° 


d> 
0 

p 

L  R 


(5-2) 


Clearly,  the  product  CB  Is  rank  deficient,  and  therefore 
an  Irregular  design  is  required.  Thus  the  design  problem  i 
non-minimum  phase  (zeros  in  the  right-half  plane).  A  mea¬ 
surement  matrix  (M)  is  used  to  produce  a  new  output: 


*  ■  jr  +  Mx  (5-3) 

1 

where  x  are  the  states  that  the  N  matrix  operates  on. 

1 

Controllability  and  observability  are  checked  using  the 
computer  program  ZERO.  The  initial  results  before  a 
measurement  matrix  is  picked  are: 

1.  No  decoupling  zeros  exist.  Thus  the  system  is 
controllable  and  observable. 

2.  The  number  of  transmission  zeros  is  equal  to  the 
number  of  states  minus  the  number  of  outputs.  The  two 
transmission  zeros  are  in  the  stable  region  (in  the  open 
left  half  Laplace  plane). 

The  addition  of  a  measurement  matrix  does  not  move  the 
transmission  zeros  to  the  unstable  region. 
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The  location  of  the  closed-loop  roots  are  given  by 
equation  (4-4).  Using  the  computer-aided-design  program 
TOTAL  allows  the  designer  to  easily  find  the  roots  of  G(s). 

Since  the  B  matrix  does  not  have  n-p  zero  rows,  where  n 
is  the  number  of  states  and  p  is  the  inputs,  the  B  matrix 


cannot  be  transformed  to  the 


form  by  simply 


rearrang¬ 


ing  rows.  However,  since  the  M  matrix  is  only  2x2,  it  is 
relatively  easy  to  pick  the  M  elements  by  trial  and  error. 
Thus  it  is  not  worth  the  effort  to  first  transform  the  B 
matrix  in  order  to  obtain  algorithms  to  pick  the  elements  of 
the  M  matrix.  A  few  steps  should  be  followed  when  selecting 

T 

an  M  matrix  when  the  system  is  not  in  (0  Ji  ]  form.  First, 

"*2 

write  the  system  such  that  the  first  1  rows  in  the  plant 
matrix  correspond  to  derivatives  of  the  states  that  the 
designer  wishes  the  N  matrix  to  operate  upon.  Secondly, 
initially  set  all  elements  to  some  nonzero  value  (0.25  is 
suggested).  Then,  keeping  in  mind  that  at  least  one  element 
in  each  column  must  be  nonzero,  zero  out  elements  one  at  a 
time  and  examine  the  effect.  The  goal  is  to  obtain  the 
sparsest  M  matrix  that  provides  the  best  results.  It  is 
significant  to  note  that  no  satisfactory  designs  with  the 

T 

system  written  in  [0  ]  form  are  obtained  in  this  thesis. 

2 

The  discussion  in  the  rest  of  this  section  pertains  to 


-\ 
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T 

selecting  en  M  metrix  for  systems  in  [0  ]  form. 

2 

The  N  matrix  must  be  added  so  that  the  matrix 

F  has  full  rank.  The  new  output  matrix  is: 

2 

F-[FF]-[C+HA  C+MA  ]  (5-4) 

12  1  11  ~2  12 

Where:  _C  *  [  _C  £  ] 

1  2 

Note  that  C  is  partitioned  such  that  C  is  pxp,  where  p  is 

2 

the  number  of  outputs  (p  must  equal  m) .  The  matrix  M  must  be 

chosen  so  that  F  has  full  rank.  Note  that  A  can  be 

2  12 

changed  by  merely  rewriting  the  order  of  the  state  equa¬ 
tions.  In  fact,  this  is  necessary  before  good  results  are 
obtained.  While  merely  changing  the  order  of  the  state 
equations  does  not  change  decoupling  or  transmission  zeros, 
it  does  change  the  selection  of  the  M  matrix.  The  M  matrix 
in  turn,  changes  output  decoupling  and  transmission  zeros. 

Next,  the  design  parameters  are  selected  by  following 
the  procedure  given  in  Chapter  IV.  Initial  attempts  did  not 
produce  good  results,  so  the  order  of  the  state  equations 
was  changed.  The  final  order  is  presented  in  equation  (5- 
1).  One  addition  is  made  to  the  Regular  design  process  for 
an  Irregular  design.  The  M  matrix  is  chosen  after  ALPHA, 
EPSILON,  and  SIGMA  elements  are  selected.  The  M  elements 
are  closely  coupled,  and  to  obtain  the  best  results  must  be 
chosen  simultaneously.  For  this  lateral  application  only 
two  M  elements  are  required,  so  selecting  the  best  set  of 
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elements  in  2-space  is  not  difficult.  However,  for  applica¬ 
tions  where  the  A  matrix  (of  size  (n-p)  x  p) 

12 

matrix  is  large,  and  the  resulting  M  matrix  is  not  sparse, 

this  can  be  more  difficult.  A  is  large,  for 

12 

example,  when  decoupling  of  the  longitudinal  and  lateral 
equations  of  motion  does  not  occur.  The  following 
guidelines  should  be  noted  when  choosing  M: 

1.  M  should  be  as  sparse  as  possible.  In  addition  to 
minimizing  design  complexity,  fewer  measurement  gains 
increase  processing  speed  and  improve  software  reliability. 

2.  H  must  be  chosen  such  that  F  B  has  full  rank  since 

2  2 

the  proportional  gain  matrix  is  computed  as  : 

-1 

K  -  (F  B  )  Y  (5-5) 

0  ~2~2  S 

3.  The  resulting  analog  transmission  zeros  (the 
sampling  period  is  not  taken  into  account)  must  be  in  the 
closed  (including  the  origin)  left  half  Laplace  plane.  The 
closed-loop  roots  approach  the  transmission  zeros  as  the 
sampling  frequency  is  increased  to  infinity.  If  the  analog 
transmission  zeros  are  in  right  half  plane,  it  may  be  possi¬ 
ble  to  keep  the  closed-loop  roots  in  the  stable  region  by 
selecting  a  long  enough  sample  period. 

4.  Decoupling  of  outputs  in  the  limit  as  the  sampling 
frequency  goes  to  infinity  is  guaranteed  only  if  the  product 

_1 

C  F  is  a  diagonal  matrix.  Note  further  that  the  diagonal 
2  2 
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terms  must  be  nonzero,  or  the  matrix  is  nonsingular. 


5.  Initially,  set  the  nonzero  H  elements  equal  to  0.25. 
Then  change  the  values  by  small  increments  and  observe  if 
the  response  is  improved.  If  M  has  many  nonzero  elements, 
then  the  elements  may  be  selected  one  at  a  time,  but 
superior  designs  are  likely  to  be  missed.  In  addition,  if  M 
elements  can  be  selected  to  have  the  same  value,  then  it  may 
be  possible  to  store  the  reduced  number  of  gains  in  regis¬ 
ters  or  "cache-memory"  .  Such  implementations  can  greatly 
improve  processing  speed. 

6.  The  H  elements  speed  up  response  and  improve 
tracking  as  the  elements  approach  zero.  Below  a  certain 
point,  though,  increasingly  smaller  improvements  in  tracking 
are  achieved  at  the  expense  of  a  degradation  in  the  control 
input  responses.  Judgement  must  be  used  in  choosing  the 
desired  tradeoff. 

Development  of  Minimum-Gain  Robust  Controller 

The  design  parameters  in  Table  5-1  are  obtained  by  using 
the  computer-aided-design  (CAD)  program  MULTI.  Once  the 
design  parameters  for  one  maneuver  are  selected,  these 
parameters  should  be  used  as  an  starting  set  of  parameters 
for  the  new  maneuver.  After  obtaining  the  optimum  set  of 
parameters  for  the  new  maneuver,  the  output  and  control 
input  responses  should  be  compared  to  the  responses  obtained 
when  using  the  initial  set  of  parameters.  This  is  done  to 
try  to  obtain  an  individual  robust  controller,  that  is,  one 
controller  that  provides  acceptable  response  for  all  maneu¬ 
vers  at  that  flight  condition.  With  this  milestone,  gain- 
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scheduling  of  the  individual  elements  of  the  K  ,  K  ,  and  M 

0  1 

elements  between  flight  conditions  is  possible. 

Once  this  is  achieved,  the  next  step  is  to  try  to  find 

one  robust  controller  that  provides  acceptable  responses  for 

all  maneuvers  at  all  flight  conditions.  In  general,  the 

fewer  the  gains,  the  better  the  design.  One  robust 

controller  greatly  eases  the  task  of  software 

implementation.  Even  finding  one  robust  gain  matrix,  K  , 

0 

K  ,  or  M  is  desirable.  Individual  gain  elements  may  be 
1 

adjusted  to  attain  this  goal.  The  individual  design  parame¬ 
ters  ALPHA,  EPSILON,  and  SIGMA  are  not  implemented  in  the 
digital  flight  control  system.  These  design  parameters  are 
used  only  to  provide  a  straightforward  method  of  obtaining 

the  best  set  of  gains,  K  and  K  . 

0  1 

Since  the  design  method  produces  a  robust  M  matrix,  only 

the  robust  K  and  K  matrices  have  to  be  determined.  Since 
0  1 

ALPHA  is  equal  to  one,  K  and  K  are  equal.  Thus  only  one 

0  1 

set  of  robust  gains  remains  to  be  determined.  Using  the 
average  of  the  K  elements  componentwise,  produces  good  re¬ 
sults  for  every  flight  condition  but  one,  0.9  Mach,  sea 
level.  It  is  noticed  that  componentwise,  the  K  elements 
at  0.9  Mach,  sea  level  are  the  smallest  of  the  five  flight 
conditions.  Thus  the  K  matrix  at  0.9  Mach,  sea  level,  is 
tried  as  a  universal  robust  controller.  Since  the  results 


appear  to  be  good  (see  Evaluation  of  Selected  Design  Parame 


ten)  ,  no  further  search  is  made.  It  is  quite  possible  that 
better  results  could  be  obtained  with  a  trimmed  robust 
contro 1 ler . 

Next,  a  minimum-gain  robust  controller  is  sought. 
Rather  than  guessing  which  gains  might  be  set  to  zero  with 
minimal  degradation  of  response,  it  is  useful  to  try  to 
obtain  some  insight.  The  rows  of  the  K  matrix  correspond  to 
the  outputs.  The  columns  of  the  K  matrix  correspond  to  the 
control  inputs.  The  lateral  K  matrix  is  thus: 


flaps 


rudder 


phi 

beta 


:] 


The  elements  whose  inputs  are  expected  to  have  minimal 
effect  upon  the  corresponding  outputs  are  set  to  zero.  The 
magnitude  of  the  element  has  nothing  to  do  with  its  likeli¬ 
hood  of  acceptable  degradation  of  response  when  set  to  zero. 
For  example,  the  magnitude  of  the  gain  corresponding  to 
(phi, rudder)  is  almost  thrice  the  value  of  the  gain 
corresponding  to  (phi, flaps).  An  interesting  sidelight  is 
that  the  gains  corresponding  to  (phi, flaps)  and 
(beta , rudder)  can  be  set  to  zero  for  beta-pointing  and 
coordinated  turn,  respectively.  This  is  not  done  for  the 
reasons  given  below. 

It  is  entirely  useless  to  obtain  zero  elements  unless  it 
is  obtained  for  a  robust  controller  for  all  maneuvers  at  all 
flight  conditions.  Otherwise,  the  memory  location  for  the 
element  must  still  be  reserved.  The  fact  the  memory  loca- 


tion  sometimes  has  a  zero  in  it  is  of  no  value. 


To  summarize,  a  minimum-gain  controller  requiring  a 
total  of  only  three  gains  is  obtained.  Two  gains  are 
required  for  the  K  matrices  corresponding  to  the  elements 
(phi, flaps)  and  (be ta , rudder )  .  The  two  gains  required  for 
the  measurement  matrix  M  have  the  same  value,  0.1,  so  only 
one  memory  location  is  required.  Therefore,  since  the  dis¬ 
tinct  gains  required  are  few,  they  may  be  stored  in  three 
registers . 

Evaluation  of  Selec  ted  Design  Parameters 

The  design  parameters  and  maximum  commanded  inputs 
selected  are  presented  in  Table  5-1.  The  measurement  matrix 
is  a  (2  x  2)  diagonal  matrix  whose  two  diagonal  elements  are 
chosen  to  be  0.1  for  all  lateral  controllers.  The  control¬ 
ler  matrices  resulting  from  the  selected  design  parameters 
are  presented  in  Table  5-2.  Figures  of  merit  are  provided 
in  Table  5-3.  Control  input  rates  are  given  in  Table  5-4. 
Calcomp  plots  of  outputs,  commanded  inputs,  control  inputs, 
and  states  are  provided  in  Figures  5-1  to  5-11.  Note  that 
only  nonzero  commanded  inputs  are  plotted.  Similar  informa¬ 
tion  is  provided  for  the  other  four  flight  conditions  in 
Appendices  C,  D,  E,  F,  G,  and  H.  Since  the  output  matrix  is 
constant  and  the  feedforward  matrix  is  zero,  only  the  plant 
and  input  matrices  for  the  other  flight  conditions  are 
different  in  Appendix  B. 

The  calcomp  plots  are  organized  in  the  following  manner: 
maneuver;  flight  condition;  outputs,  control  inputs,  and 
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states;  individual  robust  controllers,  full-gain  robust 
controller,  and  minimum-gain  robust  controller.  The  full- 
gain  controller  has  the  letter  "U"  (for  universal)  on  the 
time  axis  label.  The  minimum-gain  robust  controller  has  the 
letter  "Z"  (for  zeroed  out)  on  the  time  axis  label.  No 
plots  at  0.9  Mach,  sea  level,  are  given  for  the  full-gain 
robust  controller  (U)  since  the  individual  robust  controller 
at  this  flight  condition  is  selected  as  the  full-gain  robust 
controller  for  all  flight  conditions. 

Beta-pointing 

The  first  lateral  maneuver,  beta-pointing,  is  best 
performed  at  the  flight  condition  0.7  Mach,  15  000  feet.  A 
smoothed  ramp-up  of  15.5  degrees  of  sideslip  is  commanded 
while  roll-angle  is  commanded  to  zero.  The  largest  command¬ 
ed  inputs  (in  decreasing  order)  are  achieved  at  0.7  Mach,  15 
000  feet;  0.9  Mach,  sea  level;  0.4  Mach,  sea  level;  0.9 
Mach,  50  000  feet;  and  1.2  Mach,  15  000  feet.  The  outputs 
track  the  commanded  inputs  excellently,  with  no  overshoot  of 
the  beta  output.  In  fact,  the  commanded  beta  input  and 
output  are  indistinguishable  (see  Figure  5-1).  The  roll 
angle  undershoots  the  zero  axis,  providing  a  peak  transient 
roll  response  from  a  minimum  of  -0.10  degrees  at  0.4  Mach, 
sea  level  to  a  maximum  of  -0.8569  degrees  at  1.2  Mach,  15 
000  feet.  The  roll  response  then  asymptotically  approaches 
the  zero  axis,  with  a  maximum  roll  angle  of  -0.05  degrees  in 


four  seconds . 


TABLE  5-1 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS 


FOR  LATERAL  CONTROLLERS 


Flight  Condition  -  1.2  Mach,  15  000  feet 
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Note:  (IRC)  •  Individual  Robust  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(MGRC)-  Minimum  Gain  Robust  Controller 

Note:  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 
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TABLE  5-2 

LATERAL  CONTROLLER  MATRICES 
Flight  Condition  -  1.2  Mach,  15  000  feet 
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Note:  (IRC)  -  Individual  Robuat  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(MGRC)-  Minimum  Gain  Robust  Controller 

Note:  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0. 1,0.1)  for  all  lateral  controllers 
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LATERAL  FIGURES  OF  MERIT 
Flight  Condition  -  1.2  Mach,  15  000  feet 
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settling  tine  in  seconds 
peak  tine  in  seconds 

peak  value  (greatest  nagnitude),  angles  in  radians 


Note:  connanded  inputs  are  snoothed 
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TABLE  5-4 

LATERAL  CONTROL  INPUT  RATES 
Flight  Condition  -  1.2  Htch,  15  000  feet 
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♦Note!  maximum  differential  flaperon  rate  in  degrees/second 
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Note:  commanded  inputs  are  smoothed 

Note!  commanded  inputs  in  order  are  -  roll  angle 
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Note:  rates  presented  are  the  greatest  magnitude 
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The  full-gain  robust  controller  does  not  degrade  beta 
output  response,  except  for  minimal  transients  at  the  change 
in  slope  at  0.9  Mach,  50  000  feet,  and  1.2  Mach,  15  000 
feet.  At  0.9  Mach,  50  000  feet,  a  transient  dip  of  about 
0.4  degrees  occurs,  while  at  1.2  Mach,  15  000  feet  the  beta 
response  is  rounded  off  slightly  (see  Figure  5-2).  The  roll 
response  is  degraded  by  an  increase  in  the  transient  under¬ 
shoot.  The  peak  roll  undershoot  varies  from  -0.73  degrees 
at  0.4  Mach,  sea  level,  to  -2.02  degrees  at  1.2  Mach,  15  000 
f  eet . 

The  minimum-gain  robust  controller  provides  beta  output 
responses  that  are  indistinguishable  from  the  full-gain 
robust  controller.  The  transient  roll  response,  however,  is 
changed  from  an  undershoot  to  an  oscillation  about  the  zero 
axis.  The  oscillations  are  of  approximately  equal 
displacement  from  the  zero  axis.  The  peak  roll  displacement 
varies  from  a  minimum  of  -0.78  degrees  at  0.9  Mach,  sea 
level,  to  a  maximum  of  1.42  degrees  at  0.9  Mach,  50  000 
feet.  In  addition,  the  minimum-gain  robust  controller 
appears  to  appreciably  improve  roll  angle  transient-settling 
times.  The  roll  angle  transient  s e t t 1 ing- times  appear  to  be 
approximately  one-third  of  those  obtained  from  the  full-gain 
robust  controller.  Surprisingly,  the  roll  angle  transient 
se ttl ing- times  also  are  approximately  one-half  of  those 
obtained  by  the  individual  robust  controllers.  Thus,  if 
transient  peak  and  settling  time  are  the  criteria  used,  the 


minimum-gain  robust  controller  provides  superior  output 
response  to  the  full-gain  robust  controller. 

Next,  the  control  input  responses  are  evaluated.  An 
inherent  tradeoff  exists  between  output  response  and  control 
input  response.  If  the  designer  wishes  to  improve  control 
input  response,  the  value  of  the  measurement  elements  can  be 
increased  but  then  tracking  of  the  commanded  inputs  by  the 
outputs  is  degraded.  The  deflection  limits  of  the  flaps  are 
-19.5  and  13  degrees.  The  deflection  limits  of  the  rudder 
are  -30  and  30  degrees.  The  rudder  deflection  limits  an  in¬ 
crease  in  commanded  inputs.  Although  for  three  of  the 
flight  conditions  (0.4  Mach,  sea  level;  0.9  Mach,  50  000 
feet;  and  0.7  Mach,  15  000  feet)  the  maximum  rudder  deflec¬ 
tion  is  less  than  30  degrees,  no  appreciable  increase  in 
maximum  commanded  inputs  is  possible.  The  rudder  deflec¬ 
tions  seem  to  follow  differential  flaperon  deflections  with 
a  lag  of  about  0.1  seconds.  Every  time  a  small  change  in 
the  direction  of  flap  deflection  occurs,  the  rudder  deflec¬ 
tion  changes  direction  to  compensate  it.  Except  for  these 
transient  deflections,  the  rudder  follows  a  ramp-up  that 
provides  the  sideslip  ramp-up.  The  steady-state  flap  deflec¬ 
tion  is  approximately  two  degrees.  Mote  that  for  the  three 
flight  conditions  where  maximum  rudder  deflection  is  less 
than  the  maximum,  the  steady-state  roll  angle  is  positive, 
while  full  rudder  deflection  occurs  where  flap  deflection  is 
negative  (see  Figure  5-3). 
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The  full-gain  robust  controller  degrades  the  transient 
rudder  response.  In  addition,  the  transient  flap  response 
is  minimally  degraded  (see  Figure  5-3).  It  is  clear  that 
these  degradations  cause  an  increase  in  control  input  rates. 

The  minimum-gain  robust  controller  slightly  degrades  the 
flap  input  response  obtained  from  the  full-gain  robust  con¬ 
troller  (see  Figure  5-4).  The  best  measure  of  this  degrada¬ 
tion  is  again  the  control  input  rates. 

Note  that  the  overall  best  control  input  response  is 
obtained  from  1.2  Mach,  15  000  feet;  0.9  Mach,  sea  level; 
0.4  Mach,  sea  level;  0.7  Mach,  15  000  feet;  and  0.9  Mach,  50 
000  feet. 

The  roll  and  yav  rate  responses  have  smooth  transients 
(see  Figure  5-4).  The  steady-state  roll  rate  is  zero,  as 
expected.  Also,  as  expected,  a  negative  steady-state  yaw 
rate  is  required  to  provide  the  steady-state  sideslip.  The 
roll  rate  has  a  transient  overshoot,  while  the  yaw  rate  has 
a  transient  undershoot. 

The  full-gain  robust  controller  degrades  the  yaw  rate 
response  only  slightly,  but  causes  the  transient  roll  rate 
to  oscillate.  In  addition  the  transient  roll  rate  magnitude 
is  approximately  doubled,  but  this  is  of  little  consequence 
because  the  roll  rate  is  small.  Steady-state  responses  are 
unaffected  (see  Figure  5-5). 

The  minimum-gain  robust  controller  degrades  the 
transient  response  further  (see  Figure  5-5).  In  addition  to 
Increasing  the  magnitude  of  the  roll  rate  oscillations, 
oscillations  appear  in  the  yaw  rate  as  well.  Note,  however, 


that  the  largest  transient  roll  rate  is  nine  degrees/second 
at  0.7  Mach,  15  000  feet.  The  largest  transient  yaw  rate  is 
-22  d egr e e s / s econd  .  Se t tl ing- times  are  about  the  same  as 
for  the  full-gain  robust  controller,  but  are  three  times 
those  of  the  individual  robust  controllers. 

Again,  degradation  varies  somewhat  with  flight 
ccndition.  The  best  overall  roll  and  yaw  rate  responses  are 
provided  by  0.9  Mach,  sea  level;  1.2  Mach,  15  000  feet;  0.7 
Mach,  15  000  feet;  0.4  Mach,  sea  level;  and  0.9  Mach,  50  000 
feet. 

Maximum  permissable  control  input  rates  are  50 
degrees/second  for  the  differential  flaperons  and  105 
degrees/second  for  the  rudder.  The  maximum  negative  and 
positive  flap  rates  for  the  individual  robust  controllers 
are  -22  and  23  degrees/second,  respectively  at  0.9  Mach,  50 
000  feet.  The  maximum  negative  and  positive  rudder  rates 
for  the  individual  robust  controllers  are  -83  and  88 
degrees/second,  respectively,  at  0.9  Mach,  50  000  feet. 
Thus  the  maneuver  be ta - p o in t ing  with  individual  robust  con¬ 
trollers  produces  actual  control  input  rates  within  the 
allowable  limits. 

The  full-gain  robust  controller  has  maximum  negative  and 
positive  flap  rates  of  -22  and  15  degrees/second, 
respectively,  at  0.9  Mach,  50  000  feet.  The  maximum 
negative  and  positive  rudder  rates  are  -87  and  132 
d eg r e e s / s econd  at  0.7  Mach,  15  000  feet.  Thus  beta-pointing 
with  the  full-gain  robust  controller  requires  control  input 


rates  which  exceed  limits  for  two  flight  conditions. 

The  minimum-gain  robust  controller  has  maximum  negative 
and  positive  flap  rates  of  -32  and  31  degrees/second, 
respectively,  at  0.7  Mach,  15  000  feet.  The  maximum 
negative  and  positive  rudder  rates  are  -73  and  134 
degrees/second  at  0.4  Mach,  sea  level.  Allowable  rudder 
control  rates  are  also  exceeded  at  0.7  Mach,  15  000  feet. 

The  maximum  flap  and  rudder  rates  are  increased  in  three 
out  of  four  flight  conditions  (recall  that  the  0.9  Mach,  sea 
level,  individual  robust  controller  is  the  full-gain  robust 
controller)  when  comparing  the  individual  robust  controllers 
to  the  full-gain  robust  controller.  Therefore,  it  is  rea¬ 
sonable  to  conclude  that  in  general,  control  input  rates  are 
significantly  increased  when  comparing  performance  between 
the  individual  robust  controllers  and  the  full-gain  robust 
controller.  Practically,  the  maximum  commanded  input  must 
be  reduced  for  two  flight  conditions,  0.4  Mach,  sea  level, 
and  0.7  Mach,  15  000  feet.  The  controllers  do  not  have  to 
be  retuned.  A  reduction  in  the  commanded  steady-state  value 
of  107.  should  be  sufficient  to  bring  down  the  rudder  rates 
to  within  105  degrees/second.  Alternately,  the  maneuver 
time  must  be  lengthened  by  a  small  increment  in  order  to 
reduce  the  rudder  control  rates. 

Compared  to  the  full-gain  robust  controller,  the 
minimum-gain  robust  controller  slightly  increases  flap 
rates  in  three  of  five  flight  conditions.  Since  the  flap 
rates  are  well  within  the  50  degrees/second  limit,  this  is 


of  little  importance.  The  rudder  rates,  however,  are 
slightly  reduced  (two  degrees/second)  for  four  of  the  five 
flight  conditions.  One  flight  condition  has  an  increase  in 
rudder  rates  from  40  to  48  degrees/second.  Since  this  is 
much  less  than  the  105  degrees / second  limit,  it  is  not 
significant.  On  balance,  then,  for  beta-pointing,  minimum- 
gain  robust  controller  transient  response  is  superior  to  the 
full-gain  robust  controller  transient  response. 

Coordinated  turn 

The  second  lateral  maneuver,  coordinated  turn,  has  one 
flight  condition,  1.2  Mach,  15  000  feet,  where  the  maximum 
commanded  roll  angle  input  is  reduced  from  88.815  to  88.672 
degrees  for  the  minimum-gain  robust  controller.  This  slight 
degradation  is  compensated  by  tracking  greater  commanded 
inputs  for  shorter  maneuver  times  (not  presented).  The 
coordinated  turn  is  best  performed  (in  decreasing  order)  at 
0.7  Mach,  15  000  feet;  0.9  Mach,  sea  level;  1.2  Mach,  15  000 
feet;  0.4  Mach,  sea  level;  and  0.9  Mach,  50  000  feet.  The 
first  three  flight  conditions  command  the  roll  angle  ramp-up 
to  end  in  0.8  seconds.  The  outputs  track  the  commanded 
inputs  excellently,  with  little  overshoot  of  the  roll  out¬ 
put.  The  roll  peak  transient  ranges  from  a  minimum  of 
0.0000  degrees  at  0.7  Mach,  15  000  feet  to  a  maximum  of 
0.735  degrees  at  1.2  Mach,  15  000  feet.  In  fact,  the  roll 
output  and  commanded  input  are  indistinguishable  (see  Figure 
5-6).  The  sideslip  angle  undershoots  the  zero  axis,  pro¬ 
viding  a  peak  transient  sideslip  response  from  a  minimum  of 


-.01008  degrees  at  0.9  Mach,  50  000  feet  to  a  maximum  of 
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.01301  degrees  at  0.7  Mach,  IS  000  feet.  The  greatest 
steady-state  sideslip  angle  is  -.00131  degrees  at  0.4  Mach, 
sea  level.  Mo  steady-state  roll  angle  error  occurs. 

The  full-gain  robust  controller  slightly  degrades  roll 
output  response.  Again,  the  output  plots  are  indistinguish¬ 
able  (see  Figure  5-7).  The  maximum  roll  peak  transient  is 
1.95  degrees  at  1.2  Mach,  15  000  feet.  The  maximum  sideslip 
peak  transient  is  increased  to  0.21180  degrees  at  1.2  Mach, 
15  000  feet.  Mo  steady-state  roll  angle  error  occurs. 

The  minimum-gain  robust  controller  produces  roll  outputs 
that  improve  the  roll  peak  transient  (see  Figure  5-8).  The 
maximum  roll  peak  transient  is  1.337  degrees  at  1.2  Mach,  15 
000  feet.  The  maximum  sideslip  peak  transient  varies  from  a 
minimum  of  -0.62922  degrees  to  a  maximum  of  0.30930  degrees 
at  1.2  Mach,  15  000  feet.  Mo  steady-state  roll  angle  error 
occurs  . 

Next,  the  control  input  responses  are  evaluated.  An 
inherent  tradeoff  exists  between  output  response  and  control 
input  response.  An  improvement  in  control  input  response 
can  be  obtained  by  increasing  the  measurement  (M)  gains,  but 
this  degrades  output  response.  The  flap  deflection  limits 
are  -19.5  and  13  degrees.  The  rudder  deflection  limits  are 
-30  and  30  degrees.  The  actual  control  surface  deflections 
are  well  within  the  limits.  As  expected,  a  small  steady- 
state  rudder  deflection  of  about  -0.5  degrees  exists.  A 
steady-state  flap  deflection  of  about  -0.2  degrees  also 
exists  (see  Figure  5-8). 


The  full-gain  robust  controller  degrades  input  responses 
by  adding  greater  transient  deflections  and  oscillations. 
The  flap  peak  transient  deflection  increases  by  about  3.5 
degrees,  a  significant  amount.  These  transient  increases 
degrade  the  control  input  rate  as  well.  The  transient 
rudder  deflection  is  approximately  doubled.  Since,  the 
rudder  deflection  is  small  to  start  with,  this  is 
unimportant . 

The  minimum-gain  robust  controller  improves  input 
responses  at  0.4  Mach,  sea  level,  and  0.9  Mach,  50  000  feet 
by  decreasing  the  number  of  transient  oscillations  and  by 
reducing  the  transient  rudder  deflections.  At  0.7  Mach,  15 
000  feet  the  input  response  is  indistinguishable  from  that 
obtained  with  the  full-gain  robust  controller.  The  input 
response  obtained  at  0.9  Mach,  sea  level  is  close  to  the 
response  obtained  from  its  individual  (also  full-gain)  ro¬ 
bust  controller.  The  greatest  degradation  occurs  at  1.2 
Mach,  15  000  feet  (see  Figure  5-9),  where  the  rudder  peak 
transient  deflection  increases  by  about  4.5  degrees. 

The  roll  and  yaw  rate  responses  are  entirely  as 
expected.  The  roll  rate  ramps  up  to  a  maximum  value  that 
varies  from  a  minimum  of  37  degrees/second  at  0.9  Mach,  50 
000  feet,  to  a  maximum  of  110  degrees/second  at  0.7  Mach,  15 
000  feet.  After  staying  at  the  maximum  value,  the  roll  rate 
ramps-down  to  a  small  negative  steady-state  value  (-0.19 
degrees/second).  After  an  Initial  yaw  transient  peak  maxi¬ 
mum  of  8.03  degrees/second  at  0.7  Mach,  15  000  feet,  a 


steady-state  value  of  3.84  degrees/second  Is  established. 
The  is  necessary  to  continue  turning  the  aircraft. 


The  full-gain  robust  controller  does  not  degrade  the 
roll  rate  response.  The  yaw  rate  response  is  not 
significantly  changed.  Thus  the  full-gain  robust  controller 
does  not  degrade  the  roll  or  yaw  rate  responses  (see  Figure 
5-10)  . 

The  minimum-gain  robust  controller  also  does  not  degrade 


the  roll 

or  yaw  rate 

response . 

A  small 

osci 1 lation 

is 

introduced 

into  the 

yaw  rate 

transient , 

but  this 

is 

insignificant  (  see  Figure  5-11). 

Maximum  permissable  control  input  rates  are  50 
degrees/second  for  the  differential  flaperons  and  105 
degrees/second  for  the  rudder.  The  maximum  negative  and 
positive  flap  rates  for  the  individual  robust  controllers 
are  -53  and  67  degrees/second  at  1.2  Mach,  15  000  feet.  The 
maximum  negative  and  positive  rudder  rates  for  the 
individual  robust  controllers  are  -117  and  148 
degrees/second  at  1.2  Mach  ,  15  000  feet.  Thus  the  actual 
control  input  rates  for  a  coordinated  turn  using  individual 
robust  controllers  exceed  the  rate  limits  at  one  flight 
condition.  Since  88.82  degrees  of  roll  is  commanded  in  0.8 
seconds,  the  ramp-up  time  should  be  lengthened  to  0.9 
seconds  to  reduce  the  actual  rates  to  remain  within  limits. 

The  full-gain  robust  controller  has  maximum  negative  and 
positive  flap  rates  of  -49  and  65  degrees/second  at  1.2 
Mach,  15  000  feet.  The  maximum  negative  and  positive  rudder 


Figure  5-10 


DECREES/5EC0N0 
bo  oo  uo.co 


;c 


:c 


Figure  5-10 


rates  are  -82  and  157  degrees/second  at  1.2  Mach,  15  000 

feet.  While  the  flap  rates  are  slightly  reduced,  the  rudder 


rates  are  increased.  Thus  the  full-gain  robust  controller 
provides  a  slight  control  input  rate  degradation. 

The  minimum-gain  robust  controller  has  maximum  negative 
and  positive  flap  rates  of  -57  and  76  deg r ee s / s econd  at  1.2 
Mach,  15  000  feet.  The  maximum  negative  and  positive  rudder 
rates  are  -75  and  185  d egrees / s econd  at  1.2  Mach,  15  000 
feet.  Surprisingly,  this  is  the  only  flight  condition  to 
exceed  the  rudder  rate  limits.  Other  flight  conditions  have 
much  lower  rates,  with  minimum  negative  and  positive  rudder 
rates  of  -3  and  2  degrees/second  at  0.4  Mach,  sea  level. 
The  controller  does  not  have  to  be  retuned,  only  the  maximum 
commanded  input  value  must  be  reduced.  Except  for  the 
significant  degradation  of  rates  at  1.2  Mach,  15  000  feet, 
the  minimum-gain  controller  improves  or  does  not  degrade 
rates  compared  with  the  full-gain  controller. 


Summary 

Since  the  X-29A  demonstrator  aircraft's  longitudinal  and 
lateral  equations  of  motion  are  decoupled,  separate 
longitudinal  and  lateral  controllers  may  be  designed.  In 
addition,  the  controllers  may  be  implemented  on  separate 
computers,  thus  increasing  processing  speed  and  improving 
software  reliability  and  maintenance. 


The  outputs  selected  are  two  states,  roll  angle  and 


sideslip  angle.  The  same  output  (C)  matrix  is  used  for  both 
maneuvers  at  all  flight  conditions. 

As  in  longitudinal  commanded  inputs,  maximum  maneuvers 
are  commanded.  The  controller  is  designed  to  perform  two 
maneuvers,  b e ta -po in t ing  and  coordinated  turn.  The  maneuver 
times  and  commanded  steady-state  values  vary  from  flight 
condition  to  flight  condition. 

In  a  coordinated  turn,  the  aircraft  is  commanded  to  a 
fixed  roll  angle,  while  sideslip  is  commanded  to  zero.  In 
performing  a  coordinated  turn,  it  is  undesirable  to  allow 
sideslip.  If  sideslip  occurs  during  the  turn,  the  nose 
drops,  and  must  be  corrected  to  avoid  having  the  aircraft 
spiral.  At  all  five  flight  conditions,  the  aircraft  is 
commanded  to  a  steady-state  roll  angle  of  88.815  degrees 
while  beta  is  commanded  to  zero. 

Initially,  separate  controllers  are  designed  for  each 
maneuver  at  each  flight  condition.  However,  it  is 
discovered  that  the  design  parameters  used  to  obtain  the 
controller  for  a  coordinated  turn  at  one  flight  condition 
provide  excellent  results  in  obtaining  controllers  for  a 
coordinated  turn  at  the  other  flight  conditions. 
Additionally,  it  is  found  that  the  same  controller  obtained 
for  a  coordinated  turn  can  be  used  for  b e ta - po in t ing  at  the 
same  flight  condition.  Next,  a  successful  search  is  made 
for  one  controller  that  produces  excellent  tracking  for  both 
maneuvers  at  all  flight  conditions.  Finally,  half  of  the 


gains  in  the  robust  controller  are  set  to  zero  with  only 
minimal  degradation  of  results.  Thus  a  digital  minimum-gain 


controller  requiring  only  three  distinct  gains  is  achieved. 
Since  the  product  CB  is  rank  deficient,  ai  Irregular 


design  is  required.  Thus  the  design  problem  is  non-minimum 
phase  (zeros  in  the  right-half  Laplace  plane)  (Ref  3:7). 
An  M  matrix,  or  measurement  matrix,  is  used  to  provide 
measurements  of  the  derivatives  of  the  states  to  augment  the 
outputs. 

The  M  elements  are  closely  coupled  and  must  be  chosen 

simultaneously  to  obtain  superior  designs.  For  this  lateral 

application  only  two  M  elements  are  required,  and  selecting 

the  best  set  of  elements  in  2-space  is  not  difficult. 

However,  this  becomes  more  difficult  for  applications  where 

the  A  matrix  (of  size  (n-p)  x  p,  where  n=number  of  states 
12 

and  p=number  of  outputs)  is  large,  and  the  resulting  M 

matrix  is  not  sparse.  A  is  large,  for  example,  when  de- 

12 

coupling  of  the  longitudinal  and  lateral  equations  of  motion 
does  not  occur. 

Decreasing  the  measurement  gains  speeds  up  response  and 
improves  output  tracking  of  the  commanded  inputs.  Below  a 
certain  point,  though,  increasingly  smaller  improvements  in 
tracking  are  obtained  at  the  expense  of  degradation  of  con¬ 
trol  inputs. 

Since  the  design  method  produces  a  robust  M  matrix,  only 

the  robust  K  and  K  have  to  be  determined.  Since  ALPHA  is 
0  1 

selected  equal  to  one,  K  and  K  are  equal.  Thus  only  one 

0  1 
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set  of  robust  gains  must  be  determined.  It  is  found  that 
the  controller  which  has  the  smallest  component  gains,  0.9 
Mach,  sea  level,  produces  excellent  results  as  a  universal 
robust  controller  at  all  flight  conditions. 

Next,  a  minimum-gain  robust  controller  is  found.  The 
elements  whose  inputs  are  expected  to  have  minimal  effect 
upon  the  corresponding  outputs  are  set  to  zero.  The 
magnitude  of  the  element  has  nothing  to  do  with  its 
likelihood  of  acceptable  degradation  of  response  when  set  to 
zero.  Thus  the  elements  corresponding  to  rudder  affecting 
roll  and  differential  flaps  affecting  sideslip  are  set  to 
zero. 

Thus  a  minimum-gain  controller  requiring  a  total  of  only 
three  gains  is  obtained.  Two  gains  are  required  for  the  K 
matrices  corresponding  to  the  elements  (phi, flaps)  and 
(be ta , rudd er )  .  The  two  gains  required  for  the  measurement 
matrix  in  M  have  the  same  value,  0.1,  so  only  one  memory 
location  is  required.  Therefore,  since  only  three  distinct 
gains  are  required,  they  may  be  stored  in  registers. 

The  individual  robust  controllers  produce  excellent 
tracking  of  commanded  inputs  by  outputs.  The  full-gain 
robust  controller  minimally  degrades  tracking  by  increasing 
the  transients  at  two  flight  conditions.  The  settling  times 
are  lengthened  by  about  507. .  The  peak  times  for  three  of 
the  four  flight  conditions  (the  full-gain  robust  controller 
is  the  individual  robust  controller  at  0.9  Mach,  sea  level) 
are  the  same.  The  minimum-gain  robust  controller  reduces 


this  transient  peak  degradation  by  more  than  507.  at  1.2 


Mach,  15  000  feet.  The  peak  and  settling  times  are  slightly 
degraded  further  for  most  flight  conditions.  In  addition, 
the  plots  reveal  that  transient  settling  times  are  markedly 
improved.  The  transient  settling  times  are  half  that  of  the 
individual  robust  controllers  and  one  third  of  those  from 
the  full-gain  robust  controllers. 

The  control  inputs  are  degraded  by  the  full-gain  robust 
controller  but  less  by  the  minimum-gain  robust  controller. 
It  is  clear  that  this  causes  an  increase  in  control  input 
rates. 

For  beta-pointing,  the  individual  robust  controllers 
produce  actual  control  rates  within  allowable  limits. 
Unfortunately,  for  beta-pointing,  the  full-gain  robust 
controller  produces  a  degradation  of  control  input  rates. 
This  requires  that  the  maximum  commanded  steady-state  input 
be  reduced  at  two  flight  conditions.  A  reduction  of  no  more 
than  107.  of  the  commanded  steady-state  value  should  be 
sufficient  to  reduce  the  rudder  rates  to  within  the  limit  of 
105  degrees / s econd .  Alternately,  the  maneuver  time  must  be 
slightly  increased.  It  is  important  to  note  that  the 
controller  does  not  have  to  be  retuned.  The  minimum-gain 
Tobust  controller  ameliorates  the  degradation  of  rudder 
rates  slightly. 

On  balance  the  minimum-gain  robust  controller  provides 
superior  transient  response  to  that  obtained  from  the  full- 
gain  robust  controller.  Clearly,  since  the  degradation  is 
slight,  a  minimum-gain  robust  controller  offers  advantages 
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of  implementation  that  outweigh  the  degradation. 

The  second  lateral  maneuver,  coordinated  turn,  commands 
a  roll  angle  of  88.8  degrees  at  all  flight  conditions.  For 
the  individual  robust  controllers,  the  greatest  transient 
sideslip  is  0.013  degrees.  The  full-gain  robust  controller 
degrades  the  transient  by  increasing  the  peak  transient  to 
0.21  degrees.  The  minimum-gain  robust  controller  increases 
the  magnitude  of  the  peak  sideslip  transient  to  0.63 
degrees.  In  general,  transient  roll  overshoot  is  small  for 
the  individual  robust  controller  (maximum  0.735  degrees), 
full-gain  robust  controller  (1.95  degrees),  and  minimum-gain 
robust  controller  (1.337  degrees). 

The  full-gain  robust  controller  degrades  input  responses 
by  adding  greater  transient  deflections  and  oscillations. 
The  minimum-gain  robust  controller  improves  input  responses. 

The  actual  control  input  rates  for  a  coordinated  turn 
using  individual  robust  controllers  exceed  the  rate  limits 
at  one  flight  condition.  The  can  be  remedied  by  increasing 
the  maneuver  ramp-up  time  from  0.8  seconds  to  0.9  seconds. 
The  full  and  minimum  gain  robust  controllers  degrade  the 
rudder  rates  but  these  can  also  be  reduced  to  acceptable 
limits  by  increasing  the  ramp-up  time  or  decreasing  the  roll 
angle . 

The  transient  peak  times  are  reduced  by  about  157.  by  the 
full-gain  robust  controller.  This  is  balanced  by  about  a 
15%  increase  in  settling  times.  The  minimum-gain  robust 
controller  produces  transient  peak  times  equal  to  those  of 
the  individual  robust  controllers  but  the  settling  times  are 


lengthened  by  about  207.. 

In  general,  very  slight  degradation  occurs  when  com¬ 
paring  the  performance  of  minimum  and  full  gain  robust 
controllers  to  the  individual  robust  controllers.  Also,  the 
difference  in  performance  between  the  full  and  minimum  gain 
robust  controllers  is  minimal.  Clearly,  the  reduction  in 
performance  is  not  enough  to  warrant  using  individual  robust 
controllers  at  each  flight  condition.  The  superior  imple¬ 
mentation  advantages  of  a  digital  flight  controller  re¬ 
quiring  only  three  gains  are  considerable. 

An  error  was  found  in  MULTI  that  affected  Irregular 
designs  not  in 

creased  the  transient  and  steady-state  errors  obtained. 


£ 

B 


form.  The  "bug"  significantly  in- 


Subsequent  computer  simulations  with  the  corrected  program 
show  marked  improvement  over  the  excellent  results  presented 
here . 
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CHAPTER  VI 


CONCLUSIONS  and  RECOMMENDATIONS 

General  Comments 

Porter's  method  offers  the  advantage  that  the  analog 
theoretical  equations  have  straightforward  equivalent  digi¬ 
tal  equations  with  no  loss  of  generality.  As  yet,  the 
location  of  the  transmission  zeros  gives  no  a  priori  infor¬ 
mation  as  to  the  actual  performance  of  the  system.  Only 
general  statements  about  performance  in  the  limit  as  the 
sampling  frequency  goes  to  infinity  can  be  made. 

Whenever  the  system  model  represents  a  Regular  system, 
the  designer  has  the  option  of  using  an  Irregular  design 
approach  instead.  A  Regular  design  often  gives  better  per¬ 
formance  in  the  theoretical  limit  because  the  closed-loop 
poles  corresponding  to  the  transmission  zeros  become  unob¬ 
servable.  However,  heuristlcally  it  seems  that  Irregular 
designs  offer  better  actual  performance.  In  addition  to  the 
advantage  in  the  theoretical  limit,  the  Regular  design  has 
the  advantage  of  requiring  fewer  gains. 

Porter's  method  offers  a  quick,  fairly  straightforward 
design.  Once  the  designer  is  familiar  with  the  aircraft  and 
Porter's  method,  many  approaches  can  be  tried  in  a  short 
time  period.  Furthermore,  once  good  results  at  a  particular 
flight  condition  are  obtained,  the  same  design  parameters 
can  often  be  used  as  good  initial  values  for  other  flight 


conditions.  For  example,  one  set  of  lateral  design  parame¬ 
ters  produced  individual  robust  controllers  at  all  five 
flight  conditions. 

The  design  techniques  currently  implemented  in  MULTI  do 
not  allow  the  designer  to  have  fewer  outputs  than  inputs. 
For  a  control-configured-vehicle  (CCV)  with  many  control 
surfaces  and  hence  outputs,  it  may  be  difficult  to  determine 
what  each  commanded  input  should  be  (see  Appendix  A). 

Another  limitation  of  Porter's  method  that  should  be 
noted  is  the  requirement  that  the  number  of  outputs  be  less 
than  the  number  of  states  when  using  an  Irregular  design. 
MULTI  does  allow  the  user  to  have  as  many  outputs  as  states 

T 

when  the  system  is  not  in  10  JJ  ]  form.  This  is  a  possible 

2 

error  and  should  be  examined  (see  Appendix  A). 

Also,  MULTI  does  not  currently  have  the  capability  to 
model  the  effects  of  noise  upon  outputs.  Consider  that  the 
measured  value  of  an  output  is  equal  to  the  actual  value  of 
the  output  added  to  some  noise,  whose  statistics  are  un¬ 
known  (see  Figure  6-1): 

I  *  Cx  +  v  (6-1) 
where  Y  are  the  outputs,  C  is  the  output  matrix,  x  are  the 
states,  and  v  is  the  noise.  As  the  sampling  period 
approaches  zero,  the  gain  which  is  equal  to  the  sampling 
frequency  approaches  infinity.  For  high-gain  systems,  the 
noise  is  Increased  as  well.  Thus  the  noise  can  saturate  the 


sensors.  As  sample  periods  decrease  with  the  improvement  of 
digital  flight  control  systems,  the  noise  saturation  problem 


becomes  even  worse. 


An  inherent  tradeoff  exists  for  condi¬ 


tionally  stable  systems.  As  sampling  frequencies  increase 
to  obtain  better  precision,  attention  must  be  paid  to  sta¬ 
bility  concerns.  Fortunately,  it  would  not  be  difficult  to 
modify  MULTI  to  accept  a  stochastic  model. 

An  assumption  of  white  Gaussian  noise  can  be  made  for  v. 
In  general,  sensors  do  produce  noise  that  exhibits  the 
characteristics  of  whiteness  and  normality.  Since  Gaussian 
random  variables  are  completely  determined  by  their  first 
two  moments,  only  the  mean  and  covariance  must  be  deter¬ 
mined.  It  is  not  unreasonable  to  expect  that  the  noise 
might  have  a  bias  level  and  be  time-correlated  when  the 
sensors  are  first  turned  on.  Whiteness  refers  to  the  pro¬ 
perty  that  the  noises  are  not  time-correlated;  the  as¬ 
sumption  of  small  perturbations  in  steady  level  flight  is 
consistent  with  whiteness. 

Since  sensors  can  be  corrected  for  biases,  a  zero  mean 
assumption  of  the  sensor  noises  is  a  good  one.  To  simplify 
the  analysis,  assume  that  the  non-diagonal  terms  of  the 
covariance  matrix  are  zero.  This  is  equivalent  to 
independence  of  sensor  noises  from  one  another.  In  order  to 
determine  the  diagonal  values  of  the  covariance  matrix,  a 
determination  or  assumption  must  be  made  as  to  the  quality 
of  the  sensors  used.  In  fact,  a  sensitivity  analysis  of  the 
quality  of  sensors  required  can  be  made.  If  very  noisy 
sensors  do  not  degrade  stability  and  tracking,  then  cheaper 
sensors  for  the  digital  flight  control  system  may  be  used. 

By  the  central  limit  theorem,  a  Gaussian  random  variable 


can  be  well-modeled  by  the  sum  of  random  variables  of  uni¬ 
form  distribution,  which  are  easy  to  generate.  Random 
number  generators  could  be  used  to  generate  the  random 
variables.  This  white  Gaussian  noise  should  then  be  added 
to  the  newly  generated  states  returned  from  the  subroutine 
ODE  (Ordinary  Differential  Equations).  This  sum,  the  output 
of  the  noisy  sensor,  is  then  fed  back  to  the  summer  and 
differenced  with  the  commanded  input. 

This  is  only  correct,  though,  for  the  assumption  of  a 
perfect  truth  model  for  the  actual  system.  In  addition  to 
modeling  sensor  noise,  uncertainties  in  the  state  equation 
can  be  modeled  as  well  (see  Figure  6-2): 

e 

x  ■  Ax  +  Bu  +  Gw  (6-2) 
where  x  are  the  states,  A  is  the  plant  matrix,  B  is  the 
input  matrix,  u  are  the  inputs,  G  is  the  noise-input  matrix, 
and  w  are  white  Gaussian  noises.  Assuming  that  the  driving 
noises  in  different  state  equations  are  independent,  then 
the  covariance  of  w  is  a  diagonal  matrix.  Information  on 
wind  gusts  must  then  be  obtained  to  determine  the  covariance 
of  w.  Since  wind  gusts  are  not  time-correlated,  whiteness 
is  assured.  In  addition,  w  and  v  are  assumed  to  be  indepen¬ 
dent.  This,  however,  may  not  be  an  adequate  model  since  it 
is  expected  that  wind  gusts  also  increase  the  sensor  noise. 
Sensitivity  analysis  of  w  and  v  can  then  be  performed  to 
determine  under  what  conditions  significant  degradation  of 
the  performance  of  digital  flight  control  system  occurs. 

After  the  nev  states  are  generated  by  ODE,  and  multi- 
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plied  by  the  plant  matrix  A,  the  noise  is  added  and  the  sum 
fed  back  (see  Figure  6-2). 

Note  that  acceleration  outputs  often  require  a 
feedforward  (D)  matrix.  However,  the  state  equations  can  be 
rewritten  to  imbed  the  D  matrix  into  the  output  (C) 
matrix.  Since  Porter's  method  was  not  developed  for  systems 
with  a  D  matrix,  it  is  not  clear  whether  good  designs  under 
these  conditions  can  consistently  be  achieved.  If  the  out¬ 
put  is  the  derivative  of  one  of  the  states,  then  a  transmis¬ 
sion  zero  is  present  at  the  origin  in  the  Laplace  domain 
(equivalently,  is  situated  on  the  unit  circle  in  the  z 
domain  for  the  discrete  case),  indicating  conditional  sta¬ 
bility  (as  the  sampling  period  approaches  zero  for  the 
discrete  case).  This  poses  some  serious  limitations  since 
the  measurement  of  angles  is  difficult.  The  pilot  actually 
commands  rates  with  the  control  stick  and  feels  accel¬ 
erations  with  his  body. 

Cone lus ions 

After  an  introduction  and  description  of  the  X-29A,  a 
set  of  fifteen  longitudinal  controllers  is  developed  for  six 
maneuvers  at  five  flight  conditions.  Next,  Porter's  method 
is  applied  to  the  lateral  state-space  equations.  First, 
individual  lateral  robust  controllers  are  developed  for  each 
flight  condition.  Next,  a  single  full-gain  lateral  robust 
controller  which  achieves  good  performance  at  all  flight 
conditions  is  chosen.  Finally,  a  minimum-gain  lateral  ro¬ 
bust  controller  is  developed. 

The  straightforwardness  of  Porter's  method  is  due 


primarily  to  the  structure  of  decoupling  of  the  outputs  in 
the  limit  as  the  sampling  frequency  goes  to  infinity  for 


Regular  systems  and  for  Irregular  systems  with  a  C  F  di- 

2  2 

agonal  matrix.  Actual  decoupling  is  reasonably  good  not 
only  for  the  above-mentioned  systems  but  also  for  Irregular 

_1 

systems  where  the  C  F  matrix  is  only  diagonally  dominant. 

2  2 

As  the  order  of  the  system  increases,  the  possibility  of 
in t e r co up  1 ing  increases  as  (p-1)!,  where  p  is  the  number  of 
outputs  (assuming  the  restriction  that  the  number  of  out¬ 
puts,  p,  must  equal  the  number  of  inputs,  m).  Clearly,  such 
design  methods  can  be  extremely  complex. 

The  X-29A  demonstrator  aircraft  integrates  several 
different  advanced  technologies  of  aircraft  design.  In 
Chapter  III  forward-swept  wings,  c 1 os e - coup  led  canards, 
s tr ake- f lap s  ,  variable  camber,  fly-by-wire,  and  relaxed 
static  stability  are  discussed.  Although  knowledge  of  the 
aircraft  design  is  not  needed  to  apply  Porter's  method, 
insight  into  how  the  aircraft  "should"  behave  can  be 
helpful. 

Chapter  IV  presents  the  results  of  a  Regular  design 
application  to  develop  longitudinal  controllers  for  the  X- 
29A.  It  is  stressed  that  the  output  matrix  (C)  must  be  the 
same  for  all  maneuvers  throughout  the  flight  envelope.  The 
Regular  design,  which  requires  that  the  product  CB  be  non¬ 
singular,  is  also  a  minimum  phase  application  (transmission 
zeros  in  the  left-half  Laplace  plane).  If  it  is  non-minimum 
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phase,  then  the  system  is  conditionally  stable. 


Since  the 


closed-loop  poles  are  driven  to  the  transmission  zeros  as 
the  sampling  frequency  is  increased,  there  is  a  lower  bound 
on  the  sampling  period. 

A  description  of  the  six  longitudinal  maneuvers  designed 
for  is  provided.  At  each  flight  condition  three  controllers 
are  designed.  One  controller  for  the  direct  climb,  direct 
liftl,  and  direct  lift2  maneuvers  is  used  with  minimal 
degradation  of  performance.  One  controller  for  pitch¬ 
pointing  is  designed.  One  controller  for  vertical  transla- 
tionl  and  vertical  translation2  is  designed.  These  consoli¬ 
dated  controllers  are  used  in  the  quest  for  a  single  robust 
controller  for  all  maneuvers  at  all  flight  conditions.  An 
acceptable,  less  desirable  alternative  to  a  robust  control¬ 
ler  would  be  a  single  controller  at  each  flight  condition 
with  gain-scheduling  used  between  design  points. 

A  ten  step  general  outline  for  applying  Porter's  method 
is  presented.  A  discussion  of  design  parameters  follows. 
An  outline  to  allow  the  design  engineer  to  immediately  start 
using  the  interactive  CAD  program  MULTI  is  given. 

When  the  allowable  control  input  rates  are  exceeded  as 
occurred  for  some  of  the  longitudinal  maneuvers,  the  maximum 
commanded  inputs  must  be  reduced  or  the  maneuver  time  must 
be  lengthened.  Retuning  of  the  controllers  is  not  required. 

Since  the  X-29A  demonstrator  aircraft’s  longitudinal  and 
lateral  equations  of  motion  are  decoupled,  separate 
longitudinal  and  lateral  controllers  may  be  designed.  In 
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addition,  the  controllers  may  be  implemented  on  separate 
computers,  thus  increasing  processing  speed  and  improving 
software  reliability  and  maintenance. 

As  in  longitudinal  commanded  inputs,  maximum  lateral 
maneuvers  are  commanded.  Two  maneuvers,  b e ta - po i n t ing  and 
coordinated  turn,  are  designed  for.  The  maneuver  times  and 
commanded  steady-state  values  vary  from  flight  condition  to 
flight  condition. 

Initially,  separate  lateral  controllers  are  designed  for 
each  maneuver  at  each  flight  condition.  However,  it  is 
discovered  that  the  design  parameters  used  to  obtain  the 
controller  for  a  coordinated  turn  at  one  flight  condition 
provide  excellent  results  in  obtaining  controllers  for  a 
coordinated  turn  at  the  other  flight  conditions. 
Additionally,  it  is  found  that  the  same  controller  obtained 
for  a  coordinated  turn  can  be  used  for  be ta - po in t ing  at  the 
same  flight  condition.  Next,  a  successful  search  is  made 
for  one  controller  that  will  produce  excellent  tracking  for 
both  maneuvers  at  all  flight  conditions.  Finally,  half  of 
the  gains  in  the  robust  controller  are  set  to  zero  with  only 
minimal  degradation  of  results.  Thus  a  digital  minimum-gain 
controller  requiring  only  three  distinct  gains  is  achieved. 

Since  the  product  CB  is  rank  deficient,  an  Irregular 
design  is  required.  Thus  the  design  problem  is  non-minimum 
phase  (zeros  in  the  right-half  Laplace  plane).  An  M  matrix, 
or  measurement  matrix,  is  used  to  provide  measurements  of 
the  derivatives  of  specified  states  in  order  to  form  a  new 
output  matrix.  The  M  elements  are  closely  coupled,  and  must 


be  chosen  simultaneously  to  obtain  superior  designs.  The 

choice  of  which  states  are  fed  back  is  determined  by  the 

order  in  which  the  state  equations  are  written.  Changing 

the  order  changes  the  A  matrix;  the  measurement  matrix  (M) 
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operates  on  A  .  More  specifically,  the  first  (n-p) 

12 

derivatives  of  the  states  are  operated  on  by  M. 

It  is  found  that  the  controller  which  has  the  smallest 
component  gains,  0.9  Mach,  sea  level,  produces  excellent 
results  as  a  universal  robust  controller.  Next,  a  minimum- 
gain  robust  controller  is  found  by  setting  to  zero  those 
elements  whose  inputs  are  expected  to  have  minimal  effect 
upon  the  corresponding  outputs. 

The  individual  robust  controllers  produce  excellent 
tracking  of  commanded  inputs  by  outputs.  For  be ta - po in t ing , 
the  control  inputs  are  degraded  by  the  full-gain  robust 
controller  but  less  by  the  minimum-gain  robust  controller. 
It  is  clear  that  this  causes  an  increase  in  control  input 
rates.  In  order  to  restrain  the  control  input  rates  the 
maximum  commanded  steady-state  inputs  must  be  reduced  at  two 
flight  conditions.  Alternately,  the  maneuver  time  must  be 
slightly  increased.  It  is  important  to  note  that  again,  the 
controller  does  not  need  to  be  retuned. 

On  balance  the  minimum-gain  robust  controller  provides 
superior  transient  response  to  that  obtained  from  the  full- 
gain  robust  controller.  Clearly,  since  the  degradation  is 
slight,  a  minimum-gain  robust  controller  offers  advantages 


of  Implementation  that  outweigh  the  degradation. 


For  the  coordinated  turn,  slight  degradation  occurs  when 
comparing  the  performance  of  minimum  and  full  gain  robust 
controllers  to  the  individual  robust  controllers.  Also,  the 
difference  in  performance  between  the  full  and  minimum  gain 
robust  controllers  is  minimal.  Clearly,  the  reduction  in 
performance  is  not  enough  to  warrant  using  individual  robust 
controllers  at  each  flight  condition.  The  superior  imple¬ 
mentation  advantages  of  a  digital  flight  controller  re¬ 
quiring  only  three  gains  are  considerable.  Indirect 
control  of  rates  and  accelerations  is  possible  but 
unsa  tis  fying . 

Recommendations 

First,  suggestions  for  improvements  in  MULTI's 
capabilities  and  enhancements  are  provided.  Next,  a  series 
of  design  algorithms  that  significantly  reduce  the  trial- 
and-error  required  with  Porter's  method  are  presented.  Fi¬ 
nally,  suggestions  for  the  direction  of  future  digital 
flight  control  designs  for  the  X-29A  are  given. 

Sugges  tions  for  improvements 

MULTI,  as  it  is  presently  constituted,  is  not  written 
with  structured  programming.  In  addition  to  adding 
substantial  enhancement  of  MULTI's  capabilities,  this 
year's  group  of  theses  have  achieved  an  order  of  magnitude 
reduction  in  the  central  processing  units  (CPUs)  required  to 
run  a  simulation.  Future  improvements  to  MULTI  should 
include  "timing  runs"  in  order  to  discover  what  part  of  the 
code  is  used  most  often.  This  code  should  then  be  rewritten 


to  reflect  maximum  computational  efficiency.  In  addition, 
"blending",  or  the  use  of  assembly  language  to  speed  up  the 
parts  of  the  code  used  most  often  should  be  considered. 
These  improvements  in  efficiency  would  permit  the  additional 
CPU  time  required  for  the  addition  of  algorithms  to 
eliminate  the  t r i a  1 -and- e r r o r  presently  necessary.  To  in¬ 
crease  maintainability,  comment  statements  should  be  added 
to  MULTI.  At  present,  MULTI  is  not  well  documented. 

The  capability  to  command  a  series  of  maneuvers  should 
be  added  to  MULTI.  This  capability  allows  the  determination 
of  the  number  of  flight  condition  design  points  required  to 
adequately  design  for  the  entire  flight  envelope.  For  exam¬ 
ple,  after  the  direct  climb  maneuver  at  1.2  Mach,  15  000 
feet,  the  aircraft  has  changed  its  flight  condition  to  1.6 
Mach,  21  000  feet.  The  simulation  parameters  at  the  end  of 
the  simulation  run  are  then  used  as  the  starting  parameters 
for  a  new  maneuver.  The  results  are  then  compared  with 
those  obtained  by  using  linearized  state-space  matrices  at 
1.6  Mach,  21  000  feet  that  are  provided  by  Grumman's  nonli¬ 
near  simulation  program.  If  the  difference  in  simulations 
is  slight,  then  the  assumption  of  small  perturbations  is  not 
violated.  If  the  difference  is  great,  then  perhaps  another 
intermediate  flight  condition  between  1.2  Mach,  15  000  feet 
and  1.6  Mach,  21  000  feet  should  be  used. 

After  a  robust  con tro 1 1 e r ( s )  is  attained,  a  dynamic 


range  analysis  of  the  controller  and  measurement  matrices 
should  be  performed.  If  the  initial  heuristic  objective  of 


obtaining  controller  gains  with  magnitudes  from  0.001  to 
1000  is  achieved,  then  the  implementation  of  the  software 
can  be  done  in  single  precision.  This  increases  processing 
speed  and  addresses  the  problem  of  finite  wordlength  imple- 
men  ta  ti ons  . 

Enough  applications  of  Porter's  method  have  been  made 
that  an  analysis  of  digital  transmission  zeros  should  be 
made.  Perhaps  some  insight  into  heuris tically  determined 
good  values  of  transmission  zeros  would  assist  the  designer 
in  choosing  the  measurement  matrix  in  Irregular  designs. 
This  is  analogous  to  choosing  the  closed-loop  eigenvalues  of 
the  system.  A  priori  information  might  also  help  in  the 
selection  of  design  parameters  and  controller  gains. 

In  addition,  perhaps  an  analog  to  the  frequency-domain 
Neil-Smith  analysis  of  handling  qualities  can  be  made  with 
transmission  zeros  and  closed-loop  poles. 

The  program  ZERO  should  be  rewritten  and  incorporated 
into  MULTI.  This  would  enable  one  of  proposed  algorithms  to 
accept  only  a  measurement  matrix  that  moved  the  transmission 
zeros  into  the  unconditionally  stable  region  (within  the 
unit  circle  in  the  z-domain). 

Since  MULTI  already  calculates  the  control  input  rates, 
it  would  not  be  difficult  to  add  rate  limiting  to  MULTI 
simulations.  Assuming  that  the  commanded  maneuver  is 
reasonable,  this  would  prevent  the  designer  from  pursuing 
designs  that  require  excessive  aircraft  hinge  moments. 

A  major  improvement  to  Porter's  method  would  be  to  allow 
the  number  of  control  inputs  to  be  greater  than  the  outputs. 


If  another  independent  control  input,  such  as  the 
symmetrical  flaperons  on  the  X-29A,  is  added,  MULTI  requires 
that  another  output  be  added.  It  is  possible  to  "fool"  the 
design  method  by  commanding  flight  path  angle  when  pitch 

rate  (and  thus  pitch  angle)  and  angle-of-attack  is  already 

« 

commanded,  but  this  can  complicate  the  design.  Adding 
unneeded  outputs  increases  the  number  of  sigma  elements  and 
for  Irregular  designs  increases  the  order  of  the  measurement 
matrix. 

A  simple,  useful  addition  would  be  the  addition  of  a 
parsing  routine  to  allow  the  user  to  Interactively  enter 
several  commands  on  one  input  line  at  once. 

Design  a  Igor i thms 

The  following  is  a  series  of  algorithms  that  mimic  the 
design  procedure.  These  algorithms,  if  implemented  in 
MULTI,  would  permit  a  significant  reduction  in  the  trial- 
and-error  presently  required.  Each  option  stops  when  one  of 
two  conditions  is  met.  If  the  performance  specifications 
are  met,  then  MULTI  returns  the  answer.  If,  however,  the 
specification  is  not  met  the  designer  has  the  option  of 
using  the  best  answer  obtained,  or  giving  MULTI  more  CPU 
time.  If  the  user  does  not  specify  how  much  computer  time 
should  be  used  in  selecting  an  option,  then  a  default  value 
is  used. 

An  Implicit  assumption  in  these  algorithms  is  that  the 
model  of  small  linear  perturbations  is  accurate.  Since  the 


system  is  linear,  the  designer  can  determine  the  direction 


from  the  initial  value  that  a  design  parameter  must  be 
adjusted  to  provide  better  performance.  For  example,  if 
increasing  the  value  of  ALPHA  causes  a  degradation  in  per¬ 
formance,  then  increasing  ALPHA  further  will  only  make  mat¬ 
ters  worse.  Many  good  designs  are  possible  for  the  set  of 
unspecified  design  parameters.  However,  for  a  single  design 
parameter  there  is  a  single  best  value  (for  values  already 
assigned  to  other  design  parameters).  Thus,  part  of  the 
design  problem  is:  in  what  order  should  the  design  parame¬ 
ters  be  chosen?  Note  that  the  order  in  which  the  design 
parameters  are  chosen  will  limit  the  amount  of  further 
tuning  that  is  possible.  In  practice,  the  designer  will 
retune  each  parameter  until  significant  additional  improve¬ 
ments  in  performance  are  unobtainable. 

It  is  unknown  if  a  single  design  superior  to  others 
exists,  or  how  to  find  such  a  design.  More  likely,  many 
good  designs  exist  that  each  exhibit  superior  qualities  in 
some  respects  and  deficiencies  in  others.  Also,  diminishing 
marginal  returns  force  the  designer  to  decide  at  what  point 
additional  design  effort  is  worth  the  improvement  in 
performance . 

It  is  suggested  that  the  following  proposed  options  have 
a  prefix  of  two. 

Option  211-  This  option  selects  the  Initial  values  of 
ALPHA,  EPSILON,  and  SIGMA  based  upon  the  algorithm  presented 
in  Chapter  IV  (see  suggested  design  procedure  for  using 
MULTI ,  step  3)  . 

Option  212-  Given  a  range  of  values  for  ALPHA  and 


EPSILON,  find  the  best  set  of  ALPHA  and  EPSILON  (e.g.  ALPHA 


between  zero  and  1.5,  EPSILON  between  zero  and  four).  This 
option,  used  in  the  retuning  phase  of  design,  solves  the 
problem  in  2-space.  Note  that  the  other  design  parameters 
are  held  constant.  Thus,  no  more  than  four  simulations  must 
be  run  to  determine  the  direction  that  (ALPHA , EPSILON)  must 
be  varied. 

Option  213-  (For  Irregular  designs)  Regenerate  the 

state-space  model  in  an  equivalent  form,  with  the  rows  in 

different  order.  This  allows  a  different  set  of  states  to 

be  fed  back,  a  different  A  matrix,  and  hence  a  different  M 
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(measurement)  matrix.  The  designer  would  input  the  initial 

state  order.  Then  he  would  specify  which  states  (up  to 

maximum  of  n-p ,  where  n  is  the  number  of  states  and  p  is  the 

number  of  outputs)  he  wanted  fed  back. 

Option  214-  (For  Irregular  designs)  Given  that  the 

states  to  be  fed  back  are  already  selected,  find  ail 

acceptable  M  matrices.  Note  that  the  actual  gains  are  not 

found,  Just  the  location  of  the  nonzero  elements.  The 

default  is  to  list  the  sparsest  M  matrix. 

Note  that  F  B  must  have  full  rank.  Since  B  must  have 
2  2  2 

full  rank  (or  the  control  inputs  are  redundant)  this 

requirement  can  probably  be  met  by  choosing  F  to  have  full 

2 

rank.  It  is  possible,  though,  for  F  and  B  to  have  full 

2  2 

rank,  and  yet  their  product  could  be  rank  deficient.  Thus 


the  nonzero  M  elements  should  be  chosen  such  that  F  has 

2 

nonzero  entries  along  either,  but  not  both,  diagonals  (this 

ensures  decoupling  in  the  limit).  If  a  diagonal  F  matrix 

2 


is  not  possible,  then  M  is  chosen  to  have  the  fewest  ele¬ 
ments.  This  can  be  done  in  following  manner: 

1.  Set  each  M  element  to  an  initial  value  between 
zero  and  one  such  that  no  M  element  in  the  same  column  has 
the  same  value. 

2.  Set  one  M  element  to  zero. 

3 .  Fo  rm  F 

2 

4.  Try  to  invert  F 

2 

5.  If  successful,  keep  that  element  set  to  zero. 
If  not,  reset  the  M  element  to  its  initial  value. 

6.  Repeat  steps  2-5  until  all  M  elements  have  been 


tried . 

Option  215-  (For  Irregular  designs)  Given  that  the 

states  to  be  fed  back  are  already  selected,  print  the 

transformation  (T)  matrix  that  changes  the  B  matrix  into 

J)  J  form.  This  allows  the  designer  to  pick  the  M  elements 
B 


using  existing  algorithms. 

Option  216-  (For  Irregular  designs)  Given  a  selected  M 
matrix,  find  the  best  values  of  the  M  matrix.  This  is  done 
as  f o 1 lows : 

1 . 


Set  all  nonzero  M  elements  to  their  initial 


values  selected  la  Option  214. 


2.  Next,  find  the  solution  in  N-space,  where  N  is 
the  number  of  nonzero  M  elements.  This  is  done  by  .n  ex¬ 
haustive  search  in  which  each  element  is  increased  and 
decreased  simultaneously  with  the  other  elements  and  the 

best  improvement  in  performance  determined.  Note  that  each 

N 

iteration  will  require  2  simulations  to  be  run.  Where  N  is 
large,  MULTI  may  be  run  in  batch  or  the  M  elements  can  be 
selected  one  row  at  a  time. 

It  is  reasonable  to  assume  that  the  design  which  best 
tracks  (as  determined  by  the  figures  of  merit)  the  maximum 
commanded  input  will  also  best  track  a  less  demanding  maneu¬ 
ver.  The  less  demanding  maneuver  commands  a  fraction  of  the 
maximum  steady-state  value.  In  addition,  the  length  of  the 
maneuver  is  reduced.  In  this  thesis  all  maneuvers  and 
simulation  times  are  16  seconds  long.  A  maneuver  that  is 
only  three  seconds,  for  example,  reduces  the  amount  of  CPU 
time  required  to  evaluate  each  design. 

Option  217-  For  a  given  set  of  TREND(I),  TPEND(I), 
and  TZERO(I),  determine  find  the  maximum  VSS(I).  TREND(I), 
TPEND(I),  TZERO(I),  and  VSS(I)  are  values  used  in  Option  22 
in  MULTI  to  ramp  up  a  commanded  value  to  a  steady-state 
value  and  then  ramp  it  back  to  zero.  TREND(I)  is  the  time 
the  ramp  ends  ( s teady- s ta te  is  reached).  TPEND(I)  is  the 
time  the  input  leaves  steady-state.  TZERO(I)  is  the  time 
the  input  returns  to  zero.  This  option  finds  the  maximum 
commanded  maneuver  within  control  input  limit  and  rate  re¬ 


strictions. 


Option  218-  For  a  final  set  of  individual  robust 
controllers  with  a  constant  output  (C)  matrix,  find  a 
universal  robust  controller  of  best  performance.  This  is 
done  as  follows  : 

1.  First  try  the  controller  with  the  most  gains  of 
minimum  magnitude.  If  the  performance  is  within  specified 
limits,  stop.  Otherwise  proceed  to  step  2. 

2.  Next  try  individual  robust  controllers  at  each 
flight  condition.  After  selecting  the  individual  robust 
controller  with  the  best  Performance,  check  to  see  if  the 
performance  is  acceptable.  If  acceptable,  stop.  If  not, 
proceed  to  step  3. 

3.  Form  a  composite  robust  controller  composed  of 
the  lowest  gains,  component-wise,  from  each  of  the 
individual  robust  controllers.  If  performance  is 
acceptable,  stop.  If  not,  abandon  the  search  for  a 
universal  robust  controller. 

Option  219-  For  a  given  universal  robust  controller, 
find  a  minimum-gain  robust  controller.  For  given 
performance  specifications,  determine  which  controller  gains 
can  be  set  to  zero.  This  is  done  by  setting  the  controller 
elements  to  zero  one  at  a  time  and  determining  if  the 
performance  specifications  are  met. 

Option  221-  This  option  is  used  to  obtain  a  measure  of 
the  robustness  of  the  controller  matrices.  This  is  done  as 
f o 1 lows : 

1.  Set  the  commanded  steady-state  value  for  a  given 


maneuver  to  one-half  Its  maximum  value. 


C„- 


2.  Vary  controller  gains  one  at  a  time  until 
control  input  limits  are  reached  or  performance  degrades 
unacceptably.  This  gives  an  indication  of  how  robust  each 
gain  is  in  the  design. 

3.  Next,  vary  the  controller  gains  in  sets,  but 
only  up  to  half  the  gain  variation  found  in  step  2.  This 
gives  an  overall  robustness  measure. 

Option  222-  For  a  given  set  of  performance 
specifications,  determine  how  much  sensor  noise  can  be 
tolerated  by  varying  the  mean  and  variance  (assume  white 
Gaussian  noise).  This  gives  an  indication  of  how  good  the 
sensors  mus  t  be  . 

Option  223-  For  a  given  set  of  performance 
specifications,  determine  how  much  uncertainty  in  the  states 
can  be  tolerated  (such  as  that  caused  by  wind  gusts). 
Again,  vary  the  mean  and  variance  (assume  white  Gaussian 
noise)  . 

Note  that  Option  28  can  be  used  to  set  performance 
specifications.  Also,  MULTI  must  be  modified  to  enable  it 
to  simultaneously  store  all  controller  and  measurement 
matrices  from  each  flight  condition.  In  addition,  MULTI 
must  be  modified  to  simultaneously  store  all  maneuvers. 
Note  that  the  search  for  a  universal  robust  controller, 
search  for  a  minimum-gain  robust  controller,  and  tests  for 
robustness  must  be  repeated  for  each  maneuver. 

Future  designs  for  the  X  -  2  9 A 

In  general,  it  is  better  to  command  rates  and 


accelerations  than  angles.  The  pilot  commands  rates  with 
the  control  stick  and  can  feel  accelerations.  Angles  do  not 
have  such  a  physical  interpretation.  In  addition,  angles 
are  often  difficult  to  measure.  For  Irregular  designs  where 
the  derivative  of  a  state  is  fed  back,  noise  can  be  a  se¬ 
rious  difficulty.  Some  angles,  such  as  ang 1 e - o f - a t ta ck  and 
sideslip,  are  very  difficult  to  measure,  hence  their  deriva¬ 
tives  are  noisy.  It  is  with  these  considerations  in  mind 
that  the  following  suggestions  are  made. 

Two  longitudinal  controllers  should  be  designed  for  the 

X-29A.  One  would  provide  minimum  drag  by  incorporating 

variable  camber.  This  controller  would  have  only  two 

independent  control  inputs  and  thus  would  command  two 

outputs.  The  first  independent  control  input,  thrust,  would 

permit  control  of  u,  change  in  forward  velocity.  The 

second  independent  control  input  would  consist  of  the 

canards,  symmetrical  flaperons,  and  strake  tied  together. 

In  practice,  Grumman's  present  design  feeds  a  signal  from 

the  canards  to  the  symmetrical  flaperons  and  then  the  strake 

deflection  is  determined  by  a  computer  "table  look-up"  of 

the  symmetrical  flaperons'  deflection.  This  design  models 

the  airflow.  The  second  output  chosen  could  be  either  q, 

pitch  rate,  which  is  the  output  chosen  for  Grumman's  analog 

backup  design,  or  N  . 

z 


A  second  flight  control  system,  which  would  be  activated 
by  the  pilot's  control  stick,  is  required  for  maximum 
maneuverability.  This  design  would  treat  the  canards, 


symmetrical  flaperons,  and  strakes  as  independent  control 
inputs.  Since  four  independent  control  inputs  are 

available,  four  outputs  must  be  commanded.  Since  a  reduced- 
order  state-space  model  of  four  states  is  used,  all  four 
states  could  be  commanded  (see  Appendix  A).  The  addition  of 
pitch  angle  to  the  outputs  is  done  only  to  meet  a 
requirement  of  Porter's  method,  that  is,  the  number  of 
independent  control  inputs  must  equal  the  number  of  outputs. 
Since  a  primary-lift  control  surface,  symmetrical  flaperons, 
is  being  added  to  the  aircraft  model,  a  substantial  increase 
in  performance  over  the  results  presented  in  this  thesis  is 
expec  ted . 

For  the  lateral  mode,  an  Irregular  design  using  p,  roll 
rate,  and  r,  yaw  rate,  should  be  made  (see  Appendix  A).  In 
addition,  an  Irregular  design  with  N  and  p  as  outputs 

y 

should  also  be  attempted.  In  both  cases,  high  priority 
should  be  place  on  not  choosing  the  derivative  of  sideslip 
as  one  of  the  states  to  be  fed  back. 

The  results  of  this  thesis  may  be  useful  in  determining 
the  values  to  which  p  and  r  should  be  commanded  for  the 
first  design.  For  the  second  lateral  design  attempt, 
for  be ta - po i n t i ng  ,  N  should  be  commanded  as  a  pulse  while  p 

y 

should  be  commanded  to  zero.  For  a  coordinated  turn  N 

y 

should  be  commanded  to  zero  while  p  is  commanded  as  a  pulse. 
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Appendix  A 

Other  Design  Approaches 


This  appendix  outlines  some  design  approaches  that  were 
tried  and  found  to  be  unsatisfactory.  First,  a  discussion 
of  longitudinal  design  approaches  is  presented.  Next,  a 
discussion  of  lateral  design  approaches  is  given.  Lastly,  a 
comparison  in  performance  between  two  different  models  of 
the  X-29A  is  made.  The  results  of  the  first  model,  "data 
set2",  are  not  given.  This  thesis  presents  only  the  results 
obtained  using  "data  set3"  supplied  by  Grumman  Aerospace 
Corporation.  For  a  list  of  future  suggested  designs,  see 
Recommendations  in  Chapter  VI. 

Long i tud ina  1  Designs 

The  first  longitudinal  design  attempt  was  applied  to  a  A 
x  A  system  model.  That  is,  four  control  inputs  and  four 
outputs  were  incorporated  into  the  model.  The  control  in¬ 
puts  available  were  horizontal  canard,  flaperons,  strakes, 
and  thrust.  The  outputs  chosen  were  pitch  rate,  change  in 
forward  velocity,  ang 1 e- o f -a t tack  ,  and  normal  acceleration. 
This  is  the  same  as  the  model  chosen  for  this  thesis,  except 
that  one  additional  control  input  (flaperons)  and  output 
(normal  acceleration)  were  added. 

The  acceleration  output  added  a  D  (feedforward)  matrix 
to  the  system  model.  Consequently,  an  Irregular  design  was 
required.  Although  no  decoupling  zeros  were  present,  the 
program  ZERO  output  one  transmission  zero  in  the  right-half 
Laplace  domain.  That  was  an  error  in  the  program  ZERO, 


since  the  number  of  transmission  zeros  should  be  equal  to 
the  number  of  states  minus  the  number  of  outputs.  Since  the 
reduced-order  state-space  model  uses  only  four  states  (pitch 
angle,  pitch  rate,  change  in  total  velocity,  and  angle-of- 
attack) ,  the  number  of  transmission  zeros  should  have  been 
zero . 

_1 

If  the  system  model  is  not  in  [zero'|B  ]  form,  then  the 

2 

measurement  (M)  matrix  must  be  chosen  by  guessing,  subject 
to  the  constraint  that  the  resulting  product  FB  matrix  must 
have  full  rank.  Mote  that  since  the  number  of  states  is 
equal  to  the  number  of  outputs,  Porter's  method  collapses 

_1 

for  an  Irregular  design  applied  to  a  system  in  [zero! B  ] 

2 

form.  Porter's  method  would  require  the  M  matrix  to  have 
dimension  4x0.  MULTI,  however,  did  allow  design  attempts 

_1 

to  be  made  (maybe  because  the  system  was  not  in  [zerojB  ] 

2 

form).  Results  were  poor. 

The  second  approach  tried  was  to  imbed  first  order 
actuators  into  the  system  to  accomplish  two  objectives: 

_1 

first,  to  transform  the  system  into  [zerojB  ]  form;  and 
second,  to  increase  the  number  of  states  to  eight  while 
keeping  the  number  of  outputs  constant  at  four.  Although 
algorithms  were  now  available  for  picking  the  location  of 
nonzero  elements  in  the  M  matrix,  the  algorithms  did  not 
assist  in  picking  the  value  of  the  nonzero  elements.  Since 


the  M  matrix  was  now  a  nonsparse  4x4  matrix,  the 


s 

assumption  of  linear  behavior  requires  that  2  simulations 
(where  s  is  the  number  of  nonzero  elements  in  the  M  matrix) 
be  performed  in  order  to  determine  in  what  direction  each 
element  must  be  changed  (increased  or  decreased)  in  order  to 
improve  performance  in  the  best  possible  manner.  Since  s>8, 
guessing  the  value  of  the  M  elements  was  all  that  could  be 
practically  accomplished.  Although  a  stable  system  could  be 
achieved,  the  most  modest  maneuver  would  cause  increasing 
oscillations  in  the  control  inputs. 

Next,  a  pseudo-short-period  approximation  was  made  to 
the  system.  The  term  "pseudo"  is  used  because  the  short- 
period  roots  were  on  the  real  axis.  The  short-period 

approximation,  which  assumed  that  the  derivative  of  the 
change  in  forward  velocity  is  zero,  also  allowed  the  removal 
of  the  theta  state.  Thus  an  Irregular  design  was  applied 
to  a  5x5  system  with  three  control  inputs  and  outputs. 
The  states  were  pitch  rate,  angle-of-attack ,  canard 

deflection,  strake  deflection,  and  flaperon  deflection.  The 
control  inputs  were  the  canard  actuator,  strake  actuator, 
and  flaperon  actuator.  The  outputs  were  pitch  rate,  angle- 
of-attack,  and  normal  acceleration.  Unfortunately,  for  the 
same  reasons  as  given  before,  results  were  unacceptable 
because  of  growing  oscillations  in  the  control  inputs. 
Lateral  Designs 

Regular  designs  were  achieved  for  the  following  flight 
conditions  (the  others  were  not  attempted):  0.4  Mach,  sea 


ciLUUUL'd.Si.VL'iViWiUUK: 


level;  0.9  Mach,  50  000  feet;  and  1.2  Mach,  15  000  feet. 
The  outputs  were  roll  rate  and  yaw  rate.  These  outputs  were 
chosen  because  rates  are  easier  to  measure  than  angles. 
Unfortunately,  the  best  designs  achieved  could  not  remove  a 
large  transient  roll  angle  (e.  g.,  -4  degrees  at  0.9  Mach, 
50  000  feet)  in  the  beta-pointing  maneuver.  Similarly,  a 
large  transient  sideslip  angle  occurred  in  the  coordinated 
turn  (e.  g.,  4  degrees  at  0.4  Mach,  sea  level  with  roll 
angle  indirectly  commanded  to  90  degrees). 

Da ta  set2 

Data  set2  was  designed  for  at  three  flight  conditions: 
0.4  Mach,  sea  level;  0.9  Mach,  50  000  feet;  and  1.2  Mach,  15 
000  feet.  For  most  maneuvers,  the  maximum  commanded  input 
for  data  set2  was  407.  greater  than  the  maximum  commanded 
input  achieved  for  data  set3.  This  applies  for  both 
longitudinal  and  lateral  maneuvers.  Data  set3  was  used  in 
this  thesis  because  it  was  judged  by  Grumman  to  be  a  better 


model  of  the  X-29A  demonstrator  aircraft. 
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TABLE  D-l 


LONGITUDINAL  STATE-SPACE  MATRICES 
Flight  Condition  -  0.4  Mach,  sea  level 


IT  IT 

IT  A  (PLANT  MATRIX)  IT 

IT . - . V 

IT 

IT  0  0  0  1  IT 

IT  IT 

IT  -32.17  -.01692  -1.89  -0.3742  IT 

IT  IT 

IT  .000001876  -.0003094  -.  9991  0.9880  IT 

5T  IT 

IT  0  .0005282  10.37  -0.4340  IT 

V  IT 

IT . IT 

IT  IT 

IT  ^  (CONTROL  INPUT  MATRIX)  IT 

IT . IT 

IT  IT 

IT  0  0  0  5T 

H  IT 

IT  -.  1179  -.002723  37.82  IT 

IT  IT 

IT  -.0009614  -.0007841  -.00592  IT 

IT  IT 

IT  .0926  -.  03481  -.2680  IT 

IT  IT 

V  . IT 

^T  IT 

IT  C  (OUTPUT  MATRIX)  5T 

IT . IT 

IT  IT 

IT  0  0  0  1  IT 

IT  0  1  0  0  IT 

IT  0  0  1  0  IT 

IT _  IT 


Note:  the  states  are  (listed  in  order)  -  pitch  angle,  change 

in  forward  velocity,  ang 1 e - o f -a t tack  ,  and  pitch  rate. 

Note:  the  control  inputs  are  (listed  in  order)  -  horizontal 

canards,  s trake- f laperons ,  and  thrust. 

Note:  the  outputs  are  (listed  in  order)  -  pitch  rate,  change 

in  forward  velocity,  and  ang  1  e - o f - a 1 1 a ck  . 
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TABLE  D-2 


LONGITUDINAL  STATE-SPACE  MATRICES 
Flight  Condition  -  0.9  Mach,  50  000  feet 


IT 

IT 

IT 

A  (PLANT 

MATRIX) 

5T 

IT . 

-  -5T 

IT 

5T 

IT 

0 

0 

0 

1 

5T 

IT 

5T 

IT 

-32.02 

-  .01536 

- 

17.76 

-  .1665 

5T 

IT 

IT 

IT 

0 

- . 00008461 

- 

.  3942 

.  9978 

5T 

IT 

5T 

IT 

0 

- . 0009140 

7.143 

-  .1584 

IT 

IT 

IT 

IT . 

-  -5T 

IT 

5T 

IT 

B 

(CONTROL  INPUT  MATRIX) 

5T 

IT . 

-  -5T 

IT 

IT 

IT 

0 

0 

0 

5T 

IT 

5T 

IT 

-  .  1124 

- 

.01618 

6.443 

5T 

IT 

5T 

51 

-.0002668 

.0002067 

-  .  0006428 

IT 

5T 

5T 

5T 

.05860 

- 

.01878 

-  .  04575 

5T 

5T 

IT 

IT . 

-  -5T 

5T 

5T 

IT 

C  (OUTPUT 

MATRIX) 

IT 

5T . 

- -5T 

5T 

5T 

5T 

0 

0 

0 

1 

5T 

5T 

0 

1 

0 

0 

IT 

5T 

0 

0 

1 

0 

5T 

IT 

5T 

Note: 

the  states  are  (listed  in  order)  -  pitch 
in  forward  velocity,  angle-of-attack,  and 

angle, 

pitch 

change 

rate. 

Note: 

the  control  inputs  are  (listed  in  order) 
canards,  s trake- f laperons  ,  and  thrust. 

-  horizontal 

Note: 

the  outputs  are  (listed  in  order)  -  pitch 
in  forward  velocity,  and  angle-of-attack. 

rate, 

change 
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TABLE  D -  3 


LONGITUDINAL  STATE-SPACE  MATRICES 


Flight  Condition  -  0.7  Mach,  15  000  feet 


IT  IT 

IT  A  (PLANT  MATRIX)  IT 

IT . IT 

IT  IT 

IT  0  0  0  1  IT 

IT  IT 

IT  -32.13  -.01548  .  5001  -.3035  IT 

IT  51 

IT  -.000001132  -.0001218  -1.152  .  9917  IT 

IT  IT 

IT  .  0  .0002201  20.17  -.4928  IT 

IT  IT 

IT . - . IT 

^  IT 

V  B  (CONTROL  INPUT  MATRIX)  IT 

IT . . . - . IT 

IT  IT 

^T  0  0  0  5T 

TT  IT 

TT  -.  1877  -.005027  27.25  IT 

nr  IT 

IT  -.001136  -.0007193  -.001819  IT 

IT  IT 

IT  .  1708  -.05466  -.  1930  IT 

IT  IT 

IT . IT 

IT  5T 

5T  C  (OUTPUT  MATRIX)  5T 

* . IT 

IT  IT 

IT  0  0  0  1  IT 

IT  0  1  0  0  IT 

IT  0  0  1  0  IT 

IT _  IT 


Note: 

the  states  are  (listed  in  order)  -  pitch 
in  forward  velocity,  angle-of-attack,  and 

angle  , 
pitch 

change 

rate. 

Note: 

the  control  inputs  are  (listed  in  order) 
canards,  s t rake - f 1  a p er ons  ,  and  thrust. 

-  hor izonta  1 

Note: 

the  outputs  are  (listed  in  order)  -  pitch 
in  forward  velocity,  and  angle-of-attack. 

rate, 

change 
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TABLE  D-4 


IT 

it 

n 

IT 

IT 

IT 

IT 

IT 

IT 

tt 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

V 
IT 
IT 
IT 
IT 
IT 
TT 
IT 

V 
IT 
IT 
IT 
IT 
IT 
IT 
5T 


LONGITUDINAL  STATE-SPACE  MATRICES 
FLight  Condition  -  0.9  Mach,  sea  level 


A  (PLANT  MATRIX) 


32.17 


0637 


4 . 044 


-  .  5743 


0 

-.0002263  -3.542 

.  9856 

0 

.0002669  72. 

18 

1 . 189 

£ 

(CONTROL  INPUT  MATRIX) 

0 

0 

0 

- .6987 

-  .03225 

32.98 

- .003985 

-  .0009334 

- .001305 

.  6060 

- .1008 

- .2334 

C  (OUTPUT  MATRIX) 

0 

0  0 

1 

0 

1  0 

0 

0 

0  1 

0 

Note: 

the  states 

are  (listed  in  order)  - 

pitch  angle, 

change 

in  forward 

velocity,  ang  1  e - o f - a t ta ck  ,  and  pitch 

rate. 

Note: 


the  control  inputs  are  (listed  in  order)  -  horizontal 
canards,  s trake - f lap erons ,  and  thrust. 


Note: 


the  outputs  are  (listed  in  order)  -  pitch  rate,  change 
in  forward  velocity,  and  ang 1 e - o f - a t ta ck . 


IT 

IT 

IT 

5T 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 


TABLE  D-6 


LATERAL  STATE-SPACE  MATRICES 
Flight  Condition  -  0.9  Mach,  50  000  feet 


A  (PLANT  MATRIX) 


u 

03661 

0 

0 


0 

1 

-  .07105 

. 08658 

-9.59 

- . 8622 

3.306 

- .03038 

(CONTROL 

INPUT  MATRIX) 

0 

0 

0003309 

.0002964 

.  4374 

.1144 

02538 

-  .  03312 

. 08690 
-  .  9962 
.4814 
.01828 


C  (OUTPUT  MATRIX) 


the  states  are  (listed  in  order)  -  roll  angle,  sideslip 
angle,  roll  rate,  and  yaw  rate 

the  control  inputs  are  (listed  in  order)  -  differential 
flaperons  and  rudder. 


the  outputs  are  (listed  in  order) 
sideslip  angle. 


roll  angle  and 


v'.-'  ,v.v 

L-  • 


TABLE  D-7 

LATERAL  STATE-SPACE  MATRICES 
Flight  Condition  -  0.7  Mach,  15  000  feet 


A  (PLANT  MATRIX) 


0 

04335 

0 

0 


0 

.  2229 
■19.48 
7.912 


1 

.  04938 

-  .  2  9  7  6 

-  .  09287 


.  04944 

-  .  9988 
1.008 

-  .  09  788 


B  (CONTROL  INPUT  MATRIX) 


.0008422 

1.497 

.09646 


. 0009084 
.  2680 
-  .  08652 


C  (OUTPUT  MATRIX) 


mr. 

Note: 

the 

states  a 

re  (listed 

i  n 

order) 

roll 

angle,  sideslip 

.’a; 

ang 

le  ,  roll 

rate,  and  y 

aw 

rate 

>.V 

Note: 

th  e 

contro  1 

inputs  are 

(  li 

s  t e d  in 

order) 

-  differential 

V.V 

f  la 

perons  an 

d  rudder . 

\W0 

v*. 

Note: 

the 

sid 

outputs 
eslip  ang 

are  (listed 
1  e  . 

i  n 

order) 

-  roll 

angle  and 

V\  v'.v’JW'  WV-V  v"  v\  v'  vVO.  *'  Cv'.V  ’  ■ 


TABLE  D-8 


LATERAL  STATE-SPACE  MATRICES 
Flight  Condition  -  0.9  Mach,  sea  level 


51 

51 

51- 

51 

5T 

5T 

5T 

5T 

51 

51 

51 

5T 

5T 

5T 

51 

5T 

5T 

51 

5t 

5T 

51 

5T 

5T 

5T 

5T 

51 

51 

51 

5T 

5T 

51 

5T 

51 

51 


A  (PLANT  MATRIX) 


0 

03199 

0 

0 


0 

.4704 

74.83 

18.93 


1 

.03976 
-8.95 
-  .  2587 


B  (CONTROL  INPUT  MATRIX) 


•  .  002488 
2  .  789 
.1802 


.001397 
.6026 
-  .1831 


C  (OUTPUT  MATRIX) 


.03979 
-  .  9992 
2 .094 
-  .  1653 


1  0 

0  1 

0 

0 

0 

0 

Note: 

the  states  are  (listed  in  order)  -  roll 
angle,  roll  rate,  and  yaw  rate 

angle,  sideslip 

Note: 

the  control  inputs  are 
flaperons  and  rudder. 

(listed  in  order) 

-  differential 

Note: 

the  outputs  are  (listed 
sideslip  angle. 

in  order)  -  roll 

angle  and 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

5T 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

51 

5T 
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TABLE  E - 1 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS  FOR  LONGITUDINAL  CONTROLLER 
Flight  Condition  -  0.4  Mach,  sea  level 


Maneuver  IT  Command  Vector  v_  IT  ALPHA  5T  EPSILON  5T  SIGMA  T 

IT  *  ^T  IT  IT  (Diagonal  values)  5 


Direct  IT  1.2,  .4, 2, 3. 2 

Climb  IT  1.2,135,16,16 

IT  0 , 0 , 0 , 0 


Direct 


Direct 


. 2 ,  .  27  ,  .  5  , 
0  ,  0 , 0 , 0 
0 , 0 , 0 , 0 


1.6,  .37,2.4,3.2  IT 
0 , 0 , 0 , 0  IT 

0  ,  0 , 0 , 0  IT 


Pitch- 

V 

.8,  .036,1.2,1.6 

Po inting 

IT 

o 

o 

o 

o 

IT 

.8, .  036,16,16 

Vertical 

IT 

V 

o 

o 

o 

o 

T  rans 1  a  - 

IT 

0 , 0 , 0 , 0 

t  i  onl 

IT 

IT 

1.8, .03,3.2,4.5 

Vertical 

TT 

IT 

0 , 0 , 0 , 0 

Trans  la- 

IT 

0 , 0 , 0 , 0 

tion2 

5T 

5T 

IT 

0.8,-024,16,16 

. 2  ,  .  9  ,  .  016 


*Note:  Each  pulse  entry  m  v  has  four  parts: 

1)  The  time  (seconds)  the  input  reached  steady-state, 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady -state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


H|b-A140  982  MULTIVARIABLE  DIGITAL  FLIGHT  CONTROL  DESIGN  OF  THE 
^  X-29ACU)  AIR  FORCE  INST  OF  TECH  HRIGHT-PATTERSON  AFB  OH 
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TABLE  E-2 

DESIGN  PARAMETERS  AND  COMMANDED  INPUTS  FOR  LONGITUDINAL  CONTROLLER 
Flight  Condition  -  0.9  Mach,  50  000  feet 


IV 

Maneuver 

TV 

Command  Vector  v 

TV 

ALPHA 

TV  EPSILON 

TV 

SIGMA 

TV 

IT 

TV 

•k 

TV 

TV 

TV 

(Diagonal  values) 

TV 

IV 

TV 

TV 

TV 

TV 

TV 

IV 

Direct 

TV 

1.2,  .25,2,3.2 

TV 

TV 

TV 

TV 

IV 

Climb 

TV 

1.2,50,16,16 

TV 

1.1 

TV  0.9 

TV 

1.75,  .5,0 

TV 

IV 

TV 

0 , 0 , 0 , 0 

TV 

TV 

TV 

TV 

IV 

TV 

TV 

TV 

TV 

TV 

IV 

TV 

TV 

TV 

TV 

TV 

IT 

Direct 

TV 

TV 

TV 

TV 

TV 

IV 

Lif  tl 

TV 

0,0, 0,0 

TV 

TV 

TV 

TV 

IV 

TV 

0 , 0 , 0 , 0 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

IV 

TV 

"tv 

TV 

TV 

TV 

TV 

Direct 

TV 

1.6,  .22,2.4,3.2 

TV 

TV 

TV 

TV 

IV 

Lift2 

TV 

0,0, 0,0 

TV 

TV 

TV 

TV 

IV 

TV 

0 , 0 , 0 , 0 

TV 

TV 

TV 

TV 

IV 

TV 

TV 

TV 

TV 

TV 

IT 

TV 

TV 

TV 

TV 

TV 

IV 

Pitch- 

TV 

.8, .012,1.2,1.6 

TV 

TV 

TV 

TV 

TV 

Pointing 

TV 

0,0, 0,0 

TV 

TV 

TV 

TV 

TV 

TV 

.8,  .012,16,16 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

"TV 

TV 

TV 

TV 

TV 

Vertical 

TV 

0,0, 0,0 

TV 

TV 

TV 

TV 

TV 

Trans  la- 

TV 

0,0, 0,0 

TV 

TV 

TV 

TV 

TV 

tionl 

TV 

1.8,  .012,3.2,4.5 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

1.1 

TV  0.9 

TV 

2,1.5,  .05 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

Vertical 

TV 

0 ,0 ,0 , 0 

TV 

TV 

TV 

TV 

TV 

Trans  la- 

TV 

0 , 0 , 0 , 0 

TV 

TV 

TV 

TV 

TV 

tion2 

TV 

0.8,  .0075,16,16 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

TV 

*Note:  Each  pulse  entry  in  has  four  parts: 


1)  The  time  (seconds)  the  input  reached  steady-state, 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


TABLE  E-3 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS  FOR  LONGITUDINAL  CONTROLLER 
Flight  Condition  -  0.7  Mach,  15  000  feet 


ir 

Maneuver 

IT 

Command  Vector  v 

IT 

ALPHA 

5T  EPSILON 

IT 

SIGMA 

51 

5T 

IT 

* 

IT 

5T 

51 

(Diagonal  values) 

5T 

IT 

IT 

IT 

IT 

IT 

5T 

5T 

Direct 

IT 

1.2,. 32, 2, 3. 2 

IT 

IT 

5T 

5T 

IT 

Climb 

IT 

1.2,165,16,16 

IT 

1 

IT  0.9 

5T 

2. 2, 2,0 

5T 

IT 

IT 

0,0, 0,0 

IT 

5T 

51 

5T 

IT 

IT 

IT 

IT 

5T 

5T 

IT 

"IT 

"it 

5T 

5T 

5T 

IT 

Direct 

IT 

•  2  }  .30,  #  3 | .7 

IT 

IT 

5T 

5T 

IT 

Lif  tl 

IT 

0 ,0,0,0 

IT 

IT 

5T 

5T 

IT 

IT 

0 , 0 , 0 , 0 

IT 

IT 

5T 

5T 

IT 

IT 

IT 

5T 

5T 

5T 

5T 

IT 

IT 

5T 

5T 

5T 

IT 

Direct 

IT 

1.6,  .29,2.4,3.2 

IT 

IT 

5T 

5T 

5T 

Lif  t2 

IT 

0 ,0 ,0 , 0 

IT 

IT 

5T 

5T 

IT 

IT 

0 , 0 , 0 , 0 

IT 

IT 

5T 

5T 

IT 

IT 

IT 

IT 

5T 

5T 

IT 

5T 

IT 

IT 

5T 

5T 

IT 

Pitch- 

IT 

.8,  .032,1.2,1.6 

IT 

IT 

51 

5T 

IT 

Pointing 

IT 

0,0, 0,0 

IT 

.9 

IT  0.85 

5T 

2.65,2, .08 

5T 

IT 

IT 

.8,  .032,16,16 

5T 

IT 

5T 

5T 

IT 

IT 

IT 

IT 

5T 

5T 

IT 

IT 

IT 

5T 

5T 

5T 

IT 

Vertical 

IT 

0,0, 0,0 

•iT 

IT 

5T 

51 

IT 

Trans  la- 

IT 

0 ,0 ,0 ,0 

5T 

IT 

5T 

5T 

IT 

tionl 

IT 

1.8, .032,3. 2,4.5 

IT 

IT 

5T 

5T 

IT 

5T 

IT 

IT 

5T 

2.55,2,  .0158 

5T 

IT 

IT 

IT 

5T 

5T 

5T 

IT 

Vertical 

IT 

0,0, 0,0 

IT 

IT 

5T 

5T 

IT 

Trans  la- 

IT 

0 ,0 , 0 , 0 

IT 

IT 

5T 

5T 

IT 

tion2 

5T 

0.8,  .032,16,16 

IT 

IT 

5T 

5T 

IT 

IT 

IT 

IT 

5T 

5T 

IT 

IT 

IT 

IT 

5T 

5T 

♦Notes  Each  pulse  entry  in  v  has  four  parts: 


1)  The  time  (seconds)  the  input  reached  s teady- s ta te , 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


vt.  1  vw  v 


I'i'i'U-.U' 


yT7rrT*7wy?>r*‘  ^  w  *.  *.  *?  v  w  ■» 
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i 

i 
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TABLE  E-4 

DESIGN  PARAMETERS  AND  COMMANDED  INPUTS  FOR  LONGITUDINAL  CONTROLLER 
Flight  Condition  -  0.9  Mach,  sea  level 


it 

Maneuver 

IT 

Command  Vector  v 

IT 

ALPHA 

IT  EPSILON 

IT 

SIGMA 

IT 

ir 

IT 

* 

5T 

IT 

IT 

(Diagonal  values) 

IT 

IT 

IT 

IT 

IT 

5T 

IT 

IT 

Direct 

IT 

1.2,  .4, 2, 3. 2 

IT 

IT 

IT 

IT 

IT 

Climb 

IT 

1.2,325,16,16 

IT 

1 

IT  0.9 

IT 

2. 2, 2,0 

IT 

IT 

IT 

0,0, 0,0 

IT 

IT 

IT 

IT 

* 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

5T 

IT 

IT 

IT 

Direct 

IT 

•  2p  •  2  8  )  • 5  )  *7 

IT 

IT 

IT 

IT 

IT 

Lif  tl 

IT 

0,0, 0,0 

IT 

5T 

IT 

IT 

IT 

IT 

0 ,0 ,0 , 0 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Direc  t 

IT 

1.6, .33,2.4,3.2 

IT 

IT 

IT 

IT 

IT 

Lif  t2 

IT 

0,0, 0,0 

IT 

IT 

IT 

IT 

IT 

IT 

0 ,0,0,0 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Pitch- 

IT 

.8,  .015,1.2,1.6 

IT 

IT 

IT 

IT 

IT 

Pointing 

IT 

0 ,0,0,0 

IT 

IT 

IT 

2.2,2,  .4 

IT 

IT 

IT 

.8,  .015,16,16 

IT 

IT 

IT 

IT 

IT 

IT 

5T 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Vertical 

IT 

0,0, 0,0 

IT 

IT 

IT 

IT 

IT 

Trans  la- 

IT 

0,0, 0,0 

51 

IT 

IT 

IT 

IT 

tionl 

IT 

1.8, .02,3.2,4.5 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

0.4 

IT  3.6 

IT 

1.4, 1.2,  .016 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Vertical 

IT 

0,0, 0,0 

IT 

IT 

IT 

IT 

IT 

Trans  la- 

IT 

0,0, 0,0 

5T 

5T 

IT 

IT 

IT 

tion2 

IT 

0.8, .015,16,16 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

*Note:  Each  pulse  entry  in  has  four  parts: 


1)  The  time  (seconds)  the  input  reached  s teady- s ta te , 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 
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TABLE  E- 5 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS 
FOR  LATERAL  CONTROLLERS 

I 

Flight  Condition  -  0.4  Mach,  sea  level 


IT 

Maneuver 

IT 

Command  Vector  v_ 

IT  ALPHA 

IT 

EPSILON 

IT 

SIGMA 

IT 

IT 

IT 

IT 

IT 

IT 

(Diagonal  values) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

if 

Beta- 

"if 

IT 

"it" 

"it" 

"it 

IT 

Pointing 

IT 

IT 

IT 

IT 

IT 

H 

(IRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

if 

"it 

IT 

IT 

IT 

IT 

IT 

Beta- 

IT 

0 , 0 , 0 , 0 

IT 

IT 

IT 

IT 

IT 

Pointing 

IT 

1.2,  .24,16,16 

IT 

IT 

IT 

IT 

IT 

(FGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

if 

Beta- 

"it 

IT 

IT 

IT 

IT 

IT 

Pointing 

IT 

IT 

IT 

IT 

IT 

IT 

(MGRC ) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT  1.0 

IT 

1.0 

IT 

1.0, 1.0 

IT 

if 

"it" 

IT 

IT 

IT 

IT 

IT 

Coordina-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(IRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

1.3,1.55,16,16 

IT 

IT 

IT 

IT 

if 

"it 

0 , 0 , 0 , 0 

IT 

IT 

IT 

IT 

IT 

Coordina 

-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(FGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

if 

"it 

IT 

IT 

IT 

IT 

IT 

Coordina' 

-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(MGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Note:  (IRC)  -  Individual  Robust  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(MGRC)-  Minimum  Gain  Robust  Controller 

Note:  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 
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TABLE  E-6 

DESIGN  PARAMETERS  AND  COMMANDED  INPUTS 
FOR  LATERAL  CONTROLLERS 
Flight  Condition  -  0.9  Mach,  50  000  feet 


it 

Maneuver 

IT 

Command  Vector  v 

IT 

ALPHA 

IT 

EPSILON 

IT 

SIGMA 

V 

IT 

IT 

IT 

IT 

IT 

(Diagonal  values) 

IT 

V 

IT 

IT 

IT 

IT 

IT 

it 

Beta- 

IT 

IT 

IT 

IT 

"it 

n 

Pointing 

IT 

IT 

IT 

IT 

IT 

ir 

(IRC) 

IT 

IT 

IT 

IT 

IT 

ir 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

n 

Beta- 

IT 

0,0, 0,0 

IT 

IT 

IT 

IT 

IT 

Pointing 

IT 

0.8,0.1,16,16 

IT 

IT 

IT 

IT 

IT 

(FGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Beta 

"it 

IT 

IT 

IT 

IT 

IT 

Pointing 

IT 

IT 

IT 

IT 

IT 

IT 

(MGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

1.0 

IT 

1.0 

IT 

1.0, 1.0 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Coordina' 

-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(IRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

2.9,1.55,16,16 

IT 

IT 

IT 

IT 

IT 

IT 

0 ,0,0,0 

IT 

IT 

IT 

IT 

IT 

Coordina 

-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(FGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Coordina' 

-IT 

IT 

IT 

IT 

IT 

IT 

ted  Turn 

IT 

IT 

IT 

IT 

IT 

IT 

(MGRC) 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

IT 

Note:  (IRC)  -  Individual  Robust  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(MGRC)-  Minimum  Gain  Robust  Controller 

Note:  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 
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TABLE  E-7 


DESIGN  PARAMETERS  AND  COMMANDED  INPUTS 
FOR  LATERAL  CONTROLLERS 
Flight  Condition  -  0.7  Mach,  15  000  feet 


Maneuver  ^  Comand  Vector  v  ^  ALPHA  I"  EPSILON 


SIGMA 


I 

_ I 

Beta- 
Pointing  I 
(IRC) 


I  (Diagonal  values) 


I 


Beta-  I  0,0, 0,0 
Pointing  .  0.8,0.27,16,16 
(PGRC)  I 

beta-  . 

Pointing  I 
(M3RC) 


Coordina- I 
ted  Turn 
(IRC)  | 

_  0,8,1.55,16,16 

I  0,0, 0,0 
Coordina-1 
ted  Turn  . 

(PGRC)  I 


I 

1 


Coordina¬ 
ted  Turn  I 
(SCRC)  1 
l 

_ I 


l 


Note  1  (IRC)  -  Individual  Robust  Controller 
(PGRC)-  Full  Gain  Robust  Controller 
(ICRC)-  Minimum  Gain  Robust  Controller 

Note*  The  measurement  (M)  matrix  is  a  (2x2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 


TABLE  E-8 


DESIGN  PARAMETERS  AND  CONMANDED  INPUTS 
FOR  LATERAL  CONTROLLERS 
Flight  Condition  -  0.9  Ffech,  sea  level 


Note i  (IRC)  -  Individual  Robust  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(NDRC)-  Minimum  Gain  Robust  Controller 

Note i  The  measurement  (M)  matrix  is  a  (2x2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 


TABLE  F-l 


LONGITUDINAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.4  Mach,  sea  level 


Jfeneuver  ^  Conmand  Vector  v  1  ALPHA 

'  *  t 


T 


K. 


Direct 

1.2, .4, 2, 3. 2  “ 

I 

0.9 

Climb 

1  1.2,135,16,16 

1 

(  0,0, 0,0 

I 

1 

Direct 

I  .2, .27, .5, *7 

1 

1 

Liftl 

0,0,0, 0 

1 

.1454E+02 

.7336E-03 

0 

1  0,0, 0,0 

|  -.1816E+02 

-.1455E-00 

0 

.4401E-01 

.1915E-01 

0 

Direct 

Lift2 

*  1.6,. 37,2.4,3.2 

1  0,0, 0,0 

0,0, 0,0 

'1 

1 

1 

1 

.1 _ 

.1 - 

I 

Pitch- 

Pointing 


|  .8, .036,1.2,1.6 

1  0,0, 0,0 
|  .8, .036,16,16 


Vertical 

Transla¬ 

tion! 


T 


I  0, 0,0,0 

0,0, 0,0 

I  1.8, .03,3*2,4.5 

- 1 - 

Vertical  I  0,0, 0,0 
Transla-  0,0, 0,0 
tion2  1  0.8,. 024, 16, 16 


0*5 


.1 


.1464E+2  .1504E-2  -.8665EI 
- . 1828E2  - . 2982E+0  - . 2283E2 
.4432E-1  .3925E-1  -.2866E-1 


♦Note «  Each  pulse  entry  in  v  has  four  parts  i 


1)  The  time  (seconds)  the  input  reached  steady-state, 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-sta^e  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


*«Notei  KQ  =  ALPHA  *  ^ 


TABLE  P-2 


LONGITUDINAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.9  Mach,  50  000  feet 


♦Note  i  Each  pulse  entry  in  v  has  four  parts  i 


1)  The  time  (seconds)  the  input  reached  steady-state, 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  tine  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 

♦♦Note :  K-  *  ALPHA  ♦  K. 

0  1 


LONGITUDINAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.7  Mach,  15  000  feet 


- 7* 

*7 - 

Maneuver 

Conmand  Vector  v  .  ALPHA 

1 

1  0 

1 

*  * 

1 

1 

1 

Direct  . 
Climb  1 

1.2,02,2,3.2  1  “ 

1.2,165,16,16  ,  1 

1 

| 

1 

1 

0,0, 0,0 

1 

1 

Direct 

1 

•2, .30, *5. *7  1 

1 

.  .77UE1  .1408E-1 

0 

Liftl  l 

0,0, 0,0  1 

1  -.1231E2  -.1894 

0 

1 

0,0, 0,0  , 

|  .5085E-1  .6612E-1 

0 

Direct  5 
Lift2  1 

.  _  _  1 

1.6, .29,2.4,3.2  , 

0,0, 0,0 

I 

1 

1 

1 

0,0, 0,0  1 

1 

t 

1 

Pitch-  1 

.8, .032,1.2,1.6 

1  .7895E1  .1197E-1 

-.1812E2 

Pointing  | 

0,0, 0,0  1  .9 

1  -.126012  -.1610 

-.5613E2 

l 

.8,. 032, 16, 16  1 

j 

,  .5206E-1  .5620E-1 

>.1351 

Vertical  i 

1 

0,0, 0,0  1 

1 

| 

Transla- 

0,0, 0,0  > 

1 

tionl  ? 

I 

1.8, .032,3.2,4.5  , 

1 

,  .7597E1  .1197E-1 

-.3578E1 

1 

1  -.1213E2  -.1610 

-.1109E2 

Vertical  j 

0,0, 0,0  I 

1  .5009E-1  .5620E-1 

-.2669E-1 

Transla- 

0,0, 0,0  1 

| 

tion2  ! 

i 

1 

0.8, .032, 16, 16 

1 

1 

1 

♦Note  i  Each  pulse  entry  in  v  has  four  parts  i 


1)  The  time  (seconds)  the  input  reached  steady-state, 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 


to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 


♦♦Notei  K-  ■  ALPHA  *  K. 


TABLE  P-4 


$ 

5 

6 

l 


LONGITUDINAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.9  Mach,  sea  level 


Maneuver  ^  Comnand  Vector  v  ^  ALPHA  ^ 


Direct  .  1.2, .4,2, 3. 2 

Climb  1  1.2,325,16,16 

1  0,0, 0,0 


Direct  .2,. 28,. 5.. 7 

Liftl  '  0, 0,0,0 

I  0,0, 0,0 


•1913E1  .4871E-2  0 
-.8215E1  -.9711B-1  0 
.3250E-1  .5459 E-l  0 


Direct  1  1.6, .33,2.4,3.2 

Lift2  1  0,0, 0,0 

•  0,0, 0,0 


r>J 


Pitch- 

Pointing 


.8, .015,1.2,1.6 
0,0, 0,0 
.8,. 015, 16, 16 


.1913E1  .4871E-2  -.3760E2 

-.8215E1  -.9711E-1  -.2237E3 
.3250E-1  .5459E-1  -.1015E1 


Vertical 

0,0, 0,0  1 

1 

Transla- 

1  0,0, 0,0  1 

1 

tionl 

|  1.8, .02,3.2,4.5  . 

1  0.4 

1 

1  .1948E1  .4676E-2  -.2407E1 

-.8364E1  -.9323E-1  -.1432E2 

Vertical 

f  0,0, 0,0  • 

1  .33IOE-I  .5240E-1  -.6499E-1 

Transla- 

|  0,0, 0,0 

l 

tion2 

,  0.8, .015, 16, 16  1 

1  1 

1  1 

1 

1 

♦Note :  Each  pulse  entry  in  v  has  four  part  si 


1)  The  time  (seconds)  the  input  reached  steady-state 

2)  steady-state  value  (angles  are  in  radians), 

3)  the  time  the  input  leaves  steady-state  and  starts 
to  return  to  zero, 

4)  and  the  time  the  input  returns  to  zero. 

♦♦Note 1  K_  =  ALPHA  *  K. 


wa: 


TABLE  P-7 

LATERAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.7  Mach,  15  000  feet 


Maneuver  Conmand  Vector  v 
t 


Beta- 

Pointing 

(IRC) 


Coordina¬ 
ted  Turn 
(IRC) 


Beta- 

Pointing 

(PCRC) 


Coordina¬ 
ted  Turn 
(PGRC) 


Beta- 

Pointing 

(r.DRC) 


Coordina¬ 
ted  Turn 
(r.DRC) 


0,0, 0,0 

0.8,0.27,16,16 


0.8,1.55.16,16 

0,0, 0,0 


0,0, 0,0 

0.8,0.27,16,16 


0.8,1.55.16,16 

0,0, 0,0 


0,0, 0,0 

0.8,0.27,16,16 


0.8,1.55.16,16 

0,0, 0,0 


ALPHA 


6.3*14  -15.42 

1.792  87.78 


3.282  -8.871 

1.384  41.67 


3.282 

0 


0 

41.67 


♦Note »  K.  =  ALPHA  *  K, 

0  1 


Note:  (IRC)  -  Individual  Robust  Controller 
(PGRC)-  Full  Gain  Robust  Controller 
(TDRC)-  Minimum  Gain  Robust  Controller 

Note*  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 


TABLE  P-8 

LATERAL  CONTROLLER  MATRICES 
Flight  Condition  -  0.9  Mach,  sea  level 


Maneuver  Command  Vector  v 


Pointing  I  0,0,0, 0 
(IRC)  r  0.8,0.27,16,16 


Coordina- 1  0.8,1.55.16,16 
ted  Turn  (  0,0, 0,0 


Beta-  |  0,0, 0,0 
Pointing  0.8,0.27,16,16 
(PGRC)  1 


Coordina-*  0.8,1.55,16,16 
ted  Turn  |  0,0, 0,0 


.282  - 
1.384  41.67 


Pointing  0.8,0.27,16,16 
(M3RC)  I 


Coordina- 1  0.8,1.55,16,16 
ted  Turn  .  0,0, 0,0 


3.282 

0 


♦Note 1  Ka  =  ALPHA  *  K. 

0  1 


Note 1  (IRC)  -  Individual  Robust  Controller 
(FGRC)-  Full  Gain  Robust  Controller 
(rvERC)-  Minimum  Gain  Robust  Controller 

Notei  The  measurement  (M)  matrix  is  a  (2  x  2)  diagonal  matrix 

with  diagonal  entries  (0.1, 0.1)  for  all  lateral  controllers 
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TABLE  G-l 


LONGITUDINAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.4  Mach,  sea  level 


♦Note :  settling  time  in  seconds 
**Notei  peak  time  in  seconds 

♦♦♦Note t  peak  value  (greatest  magnitude),  angles  in  radians 
Notet  coirmanded  inputs  are  smoothed 

Notet  coirmanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of-attack 


J. 

V 


v 

>• 

a 

1 


TABLE  G-2 

LOfCITUDINAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.9  Mach,  50  000  feet 


Maneuver  l  Command  Vector  v 


Direct 

Liftl 


Vertical 

Transla- 

tion2 


.2, .19, .5. .7 
0,0, 0,0 
0,0, 0,0 


15.92  1  .35 

15.92  I  .7 
15.92  ,  .7 


1.6,. 22, 2.4, 3. 2  15.92  .  2.275 
0,0, 0,0  ,  15.92  2.80 
0,0, 0,0  15.92  I  2.8 


.2681 

50.30 

.3226 


.  .1932 
-.2042 
•  .0857 


.2379 

-.6226 

.2975 


Pitch- 

1 

1 

.8, .012,1.2,1.6 

1  15.92 

1.05 

.0129 

Pointing 

1 

0,0, 0,0 

15.92 

1.40 

1  -.0082 

1 

.8, .012, 16, 16 

|  7.00 

1.40 

.0141 

Vertical 

T 

0,0, 0,0 

1 

15.92 

2.275 

1  .0006 

Transla¬ 

1 

0,0, 0,0 

•  15.92 

2.275 

1  -.0024 

tion! 

1 

1.8, .012,3.2,4.5 

15.92 

1 

3.15 

.0120 

0, 0,0,0  15.92  1  1.4 

0,0, 0,0  ,  15.92  |  1.575 

0.8,. 0075, 16, 16  1.925  2.80 


♦Note*  settling  time  in  seconds 
**Notei  peak  time  in  seconds 

♦♦♦Note «  peak  value  (greatest  magnitude),  angles  in  radians 
Note*  commanded  inputs  are  smoothed 

Note*  commanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of -attack 


V f*i rfi  **i  1T1  *‘m *n -W mlT n  f  i  ‘V*^ i  * \  *‘  **  .**«  **i 


TABLE  G-3 

LONGITUDINAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.7  Mach,  15  000  feet 


♦Notei  settling  time  in  seconds 
♦♦Note «  peak  time  in  seconds 

♦♦♦Note «  peak  value  (greatest  magnitude),  angles  in  radians 


Note  i  comnanded  inputs  are  smoothed 
Note«  conmanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of -attack 
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TABLE  G-4 

LONGITUDINAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.9  Mach,  sea  level 


Maneuver 

T 

Command  Vector  v 

r-  ,T 

t  ** 

1  fj  *** 

l 

-l- 

S  1 

!  > 

P 

1  P 

1 

Direct 

1.2, .4, 2,j. 2 

1.5 

Climb 

i 

1.2,325,16,16 

1  1.35 1 

1.5 

1  326.6 

1 

0,0, 0,0 

6.6  | 

1.8 

|  .0993 

Direct 

.2,. 28,. 5,. 7 

1 

l  15.00 

.45 

1  .3232 

Liftl 

| 

0,0, 0,0 

15.00  1 

.6 

1  -.0507 

0,0, 0,0 

1  15.00 

.6 

'  .0677 

Direct 

1 

1.6,. 33, 2.4, 3.2 

,  15.00  1 

1.95 

1  .3683 

Lift2 

1 

0,0, 0,0 

15.00  j 

7.65 

1  -.1288 

1 

0,0, 0,0 

1  15.OO 

2.25 

1  .0999 

Pitch- 

J 

.8, .015,1.2,1.6 

1  15.00  * 

0.9 

1  .0231 

Pointing 

1 

0,0, 0,0 

i  15.00  , 

1.35 

-.0196 

I 

1 

.8, .015, 16, 16 

'  3.75 
,  1 

15.00 

1  .0150 

Vertical 

1 

Of  Of  Of  0 

,  15.00  1 

1  15.00  , 

1.95 

.0024 

Transla- 

0,0, 0,0 

2.4 

1  -.0003 

tionl 

1 

1.8,. 02, 3.2, 4.5 

15.00  1 

3.45 

|  .0154 

Vertical 

0,0, 0,0 

'  15.00  . 

1  15.00  1 

1.05 

1  .0023 

Transla- 

1 

0,0, 0,0 

1.65 

I  -.0002 

.0150 

tion2 

1 

0.8, .02, 16, 16 

|  7.05  | 

15.00 

.1 _ L 

*Note«  settling  time  in  seconds 
**Notet  peak  time  in  seconds 

***Note :  peak  value  (greatest  magnitude),  angles  in  radians 


Notet  commanded  inputs  are  smoothed 

Notei  commanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of -attack 
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TABLE  G-5 

LATERAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.4  Mach,  sea  level 


♦Note:  settling  time  in  seconds 
♦♦Note:  peak  time  in  seconds 

♦♦♦Note i  peak  value  (greatest  magnitude) ,  angles  in  radians 

Note:  cornnonded  inputs  are  smoothed 
Note i  coirmanded  inputs  in  order  are  -  roll  angle,  and 
sideslip  angle 
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TABLE  G-6 


I 

i 


LATERAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.9  Mach,  50  000  feet 


♦Note  i  settling  tiro  in  seconds 
♦♦Notes  peak  tiro  in  seconds 

♦♦♦Notes  peak  value  (greatest  magnitude),  angles  in  radians 

Notes  commanded  inputs  are  smoothed 
Notes  cosnnanded  inputs  in  order  are  -  roll  angle,  and 
sideslip  angle 


TABLE  G-7 

LATERAL  FIGURES  OF  MERIT 


Flight  Condition  -  0.7  Mach,  15  000  feet 


Maneuver  ^  Command  Vector  v  t_  *  t 
.  “  I  S  ,  ? 


M  *** 
? 


Beta- 

Pointing 

(IRC) 

1 

1 

1  15.00 
i1-2 

1  0.9 
15.00 

1- 

1  -.0082 
(  .2722 

Beta- 

! 

0,0, 0,0 

1  15.00 

1  0.9 

1  -.0172 

Pointing 

|  .8, .27, 16, 16 

1.65 

1  15.00 

,  .2723 

(FGRC) 

1 

• 

1 

1 

1 

-l - 

Coondina- 1 

0.8,1.55.16,16 

1 

1  u°5  1 

1 15 

1.5500 

ted  Turn 

0,0, 0,0 

1 9.9 

0.3 

1 

|  .0011 

(FGRC)  1 

1 

1 

Coordina¬ 
ted  Turn  ' 

1  1.05 

1 

1  l-2 

1  1.5504 

1  10.05 

1  .45 

.0011 

1 

(IDRC)  | 

1 

1 

1 

1 

1 

1 

♦Note  i  settling  time  in  seconds 
♦♦Note i  peak  time  in  seconds 

♦♦♦Notes  peak  value  (greatest  magnitude),  angles  in  radians 

Notes  corranded  inputs  are  smoothed 
Notes  contended  inputs  in  order  are  -  roll  angle,  and 
sideslip  angle 


a -4  **•/>  i  n.1  * 


•V;'l 


m 


v 


tti:UAActoM:W»tM^L'a:ai>toL’>ata 


TABLE  G-8 

LATERAL  FIGURES  OF  MERIT 
Flight  Condition  -  0.9  Mach,  sea  level 


Maneuver  Coirmand  Vector  v 


Pointing  1 

1  15.00  0.9 

(IRC)  , 

|  1.8  |  15.00 

•.0310 

.2738 


Pointing 

(FGRC) 


vy w| v 


.8,. 27, 16, 16 


Beta- 

Pointing 


I  15.00  *  0.9 
*  -  *  -  /> 


-.0137 

.2738 


Coordina-  ■ 

1  1.05 

1.35 

ted  Turn 
(IRC)  l 

i 

,  6.6 

1 

1 

1  0.6 

1 

Coordina- | 

0.8,1.55.16,16 

1 

ted  Turn 

0,0, 0,0 

1 

(FGRC)  1 

1 

I 

1 

Coordina- 1 

1 

1.05 

1 

.  1-35 
1  .30 

ted  Turn  , 

1  10.95 

*Hotei  settling  time  in  seconds 
**Note»  peak  time  in  seconds 

***Notei  peak  value  (greatest  magnitude),  angles  in  radians 


Note  i  comnanded  inputs  are  smoothed 
Notet  commanded  inputs  in  order  are  -  roll  angle,  and 
sideslip  angle 


Aw 
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TABLE  H-0 


Note  (+)  Canard  -  trailing  edge  down 

(+)  Right  Flaperon  -  trailing  edge  up 
(+)  Left  Flaperon  -  trailing  edge  down 
(+)  Strake-flap  -  trailing  edge  down 
(+)  Rudder  -  trailing  edge  left 

Information  on  thrust  dead-zone,  time-lags,  rates,  and 
maximum  output  at  various  flight  conditions  is  avail¬ 
able  from  the  Grumnan  reduced-state  simulation  program 
Also,  the  limits  on  the  inboard  and  outboard  flaperons 
have  been  changed  and  are  different  for  the  digital 
and  analog  reversion  mode. 


TABLE  H-l 


LONGITUDINAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.4  Mach,  sea  level 


Maneuver 

T 

Command  Vector  v  *  * 

• 

r 

- f - 

** 

0 

c 

W**  * 

1  c 

l 

.L 

1 

1 

K 

1 

$s 

1 

1 

6t 

Direct 

Climb 

1 

1 

1.2, .4, 2, 3. 2  . 
1.2,135,16,16  1 
0,0, 0,0  1 

46 

! 

I 

-191 

l 

1 

21 

Direct 

Liftl 

1 

1 

1 

•  2 ,  •  27 ,  •  5 » •  7  1 

0,0, 0,0  1 

0,0, 0,0  1 

143 

I 

| 

176 

1 

1 

0.5 

Direct 

Lift2 

T 

I 

1 

( 

1.6,07,2.4,3.2 

0,0, 0,0  1 

0,0, 0,0  1 

59 

1 

1 

-74 

T 

1 

_ 1 _ 

.31 

Pitch- 

Pointing 

i 

1 

-  u 

.8, .036,1.2,1.6  1 
0,0, 0,0  I 

.8, .036,16,16  1 

-17 

"T" 

1 

1 

-39 

1 

1 

-.05 

Vertical 

Transla- 

tionl 

r 

1 

1 

0,0, 0,0  1 

0,0, 0,0  1 

1.8, .03,3.2,4.5  | 

12 

I 

1 

27 

1 

1 

1 

.04 

Vertical 

Transla- 

tion2 

I 

1 

-L 

I 

0,0, 0,0  l 

0,0, 0,0 

0.8,. 027, 16, 16  1 

\ 

1 _ 

-16 

1 

1 

1 

1 

-37 

1 

1 

1 

_ 1 _ 

-.05 

•Note*  maximum  canard  rate  (magnitude)  in  degrees/second 
••Note*  maximum  strake  rate  (magnitude)  in  degrees/second 
•••Note*  maximum  thrust  rate  (magnitude)  in  hundreds  of 
lbs  ./second 

Note*  commanded  inputs  are  smoothed 

Note*  comnanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of-attack 


TABLE  H-2 


LONGITUDINAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.9  Ifech,  50  000  feet 


♦Note i  maximum  canard  rate  (magnitude)  in  degrees/second 
♦♦Note*  maximum  strake  rate  (magnitude)  in  degrees/second 
***Notet  maximum  thrust  rate  (magnitude)  in  hundreds  of 
lbs ./second 

Note:  commanded  inputs  are  smoothed 

Note*  commanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of -attack 
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TABLE  H-3 

LONGITUDINAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.7  Mach,  15  000  feet 


*Notei  maximum  canard  rate  (magnitude)  in  degrees/second 
**Notet  maximum  strake  rate  (magnitude)  in  degrees/second 
***Notej  maximum  thrust  rate  (magnitude)  in  hundreds  of 
lbs  ./second 

Note i  commanded  inputs  are  smoothed 
Notej  commanded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of-attack 
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TABLE  H-4 


Direct 


LONGITUDINAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.9  Mach,  sea  level 


Maneuver  Comnand  Vector  v  |  *  • 

& 


1.2,325,16,16 

0,0, 0,0 


"f  **  * 

,  £ 


.2, .28, .5, .7 
0,0, 0,0 
0,0, 0,0 


Vertical  0,0, 0,0 

Transla-  1  0,0, 0,0 
tionl  1.8, .02,3.2,4.5 

1 

1 

4  1 

i 

,  1 

16  .09 

!  1 

Vertical  1  0,0, 0,0 

Transla-  0,0, 0,0 

tion2  |  0.8, .02,16,16 

1 

-4  1 

1  1 

1  -17  ,  -.09 

1 

_ 1 

1  1 

1  _ 1  _ 

♦Note :  maximum  canard  rate  (magnitude)  in  degrees/second 
♦♦Note s  maximum  strake  rate  (magnitude)  in  degrees/second 
♦♦♦Notei  maximum  thrust  rate  (magnitude)  in  hundreds  of 
lbs ./second 

Notei  commanded  inputs  are  smoothed 
Notei  corananded  inputs  in  order  are  -  pitch  rate,  change 
in  velocity,  and  angle-of -attack 
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V  ■.-  -■-  -.-  v-  ^•,.+jr.ir.  «-.«-.r-.i \r,rmr.i  -r*~:  r^t..  , 


TABLE  H-5 

LATERAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.4  Mach,  sea  level 


Notei  commanded  inputs  are  smoothed 

Note j  commanded  inputs  in  order  are  -  roll  angle,  and  side¬ 
slip  angle 

Note*  rates  presented  are  the  greatest  magnitude 
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TABLE  H-6 

LATERAL  CONTROL  INPUT  RATES 
Flight  Condition  -  0.9  Mach,  $0  000  feet 


I 


8 


JOJ 


Pointing 
IRC 


Pointing  .8,.l,l6,l6 
(RJRC)  | 


* 

*«  . 

i 

23 

88 

-22 

95 

23 

93 

Beta-  I 
Pointing  , 
■DRC 


Coordina- 1 
ted  Turn  , 
(IRC)  I 


Coordina-  2.9.1.55.16,16 
ted  Turn  I  0,0, 0,0 
(PGRC) 

Coordina- 1 
ted  Turn 
(rCRC)  I 


•Note  i  maximum  differential  flaperon  rate  in  degrees/second 
**Notei  maximum  rudder  rate  in  degrees/second 

Note i  commanded  inputs  are  smoothed 

Note j  commanded  inputs  in  order  are  -  roll  angle,  and  side¬ 
slip  angle 

Note i  rates  presented  are  the  greatest  magnitude 
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TABLE  H-7 

LATERAL  CONTROL  INPUT  RATES 


FHght  Condition  -  0.7  Mach,  15  000  feet 


Note  i  commanded  inputs  are  smoothed 

Note*  commanded  inputs  in  order  are  -  roll  angle,  and  side¬ 
slip  angle 

Notei  rates  presented  are  the  greatest  magnitude 


TABLE  H-8 


LATERAL  CONTROL  INPUT  RATES 
FHght  Condition  -  0.9  Mach,  sea  level 


Note «  commanded  inputs  are  smoothed 

Note i  coitranded  inputs  in  order  are  -  roll  angle,  and  side¬ 
slip  angle 

Note i  rates  presented  are  the  greatest  magnitude 
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